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I. INTRODUCTION 


The overall goal of this Investigation has been to study the feasibi- 
lity of modeling the magnetic fields produced by certain electrical 
currents flowing In the earth's Ionosphere-magnetosphere system. Vector 
magnetic field measurements from the near-polar orbiting Magsat satellite 
contain, In addition to the main geomagnetic field and crustal anomaly 
fields, contributions that arise from these external currents. In 
fulfilling the ultimate goals of the Magsat project. It Is desirable that 
the external current effects be Identified In the observations and sub- 
sequently separated from the Internal field. The objective of this 
Investigative effort has been to determine the capability of a modeling 
procedure to facilitate the separation of these external and Internal com- 
ponents. 

The approach of this feasibility study was to develop forward modeling 
procedures through which the magnetic effects of model currents may be 
derived. It Is Intended to enable the modeling of, separately, the 
equatorial electrojet, Sq currents, and the effects due to auroral zone 
and polar cap currents Including the high latitude Ionosphere-magnetosphere 
coupling currents. Candidate current systems were devised and resulting 
'typical' magnetic field signatures calculated for comparison with Magsat 
observations. 

The effort has successfully demonstrated that an efficient, com- 
putationally economical, yet accurate, model of the magnetic effects of 
distributed space currents could be developed. 
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II. PROJECT EVOLUTION 


A. Historical Perspectives 

Prior to the Magnetic Fields Satellite (MAGSAT), previous satellite 
measurements of near-earth vector magnetic fields had been severely 
limited, both In number and In quality. The Navy TRIAD satellite carried a 
tri-axial n.agnetometer into a near earth-polar orbit at 800 km. A number 
of Investigators have used the TRIAD data In studying the signatures of 
high latitude field-aligned currents (^rmrtrong and Zmuda, 1970, 1973; 

Zmuda and Armstrong, 1974; lijima and Potemra, 1976a, b, 1978; Suglura and 
Potemra, 1976). Due to a relatively poor knowledge of the spacecraft atti- 
tude, however, TRIAD data has not been of significance for detailed 
modeling of the earth's main magnetic field. 

At the time of Magsat, magnetic field observations made from the 
ISIS-2 satellite were also being successfully employed to study high lati- 
tude perturbations associated with the Birkeland currents (Klumpar et al , 
1976; Burrows et al . , 1976; McDlarmid et al . , 1977, 1978). Deviations on 
the order or 1000 nT In the component transverse to the main field are not 
uncommon In the Birkeland current regions. The ISIS data is unique by vir- 
tue of the fact that direct comparisons of magnetic perturbations with 
charged particle and plasma observations are being made simultaneously from 
a single sctellite (Klumpar et al . , 1976). Results of these comparisons 
illustrate that the magnetic signatures attributed to Birkeland currents 
are not always unique to locations where auroral particle precipitation is 
present (Klumpar, 1976, 1979). 

The first near-earth satellite measurements of the earth's global 
magnetic field usable for providing large scale field models were those of 
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the Polar Orbiting Geophysical Observatory (POGO) series. These satellites 
measured only the scalar magnetic field primarly during the last half of 
the 1960's over an altitude range from 400 km to 1510 km. Several magnetic 
field models based upon the POGO data had been published. Langel (1974a, b, 
1975a, b) carried out a detailed analysis of the POGO scalar magnetic field 
data at high latitudes. His analysis shows that there are extended regions 
over which the field Is enhanced or depressed with respect to the model 
magnetic field. These variations can only partially be Interpreted In 
terms of horizontal Ionospheric currents. Magnetospheric current systems 
undoubtedly produce an Important contribution to these field variations, 
and would be no less Important In the MAGSAT data. 

The MAGSAT mission would represent the first time U.S. scientists 
would have access to detailed, accurate vector magnetic fields from a near- 
earth orbiting satellite. 

B. Proposed Investigation 

In response to NASA's Announcement of Opportunity from the Office of 
Space and Terrestrial Applications (A.O. No. OSTA ’8-1) dated September 1, 
1<J’9 for "Data Use Investigations for the Magnetic Fields Satellite 
(MAGSAT) Mission" a pr.nposal was submitted from the University of Texas at 
Dallas in January, 1979. The proposed study, entitled "Investigation of 
the Effects of External Current Systems on the MAGSA’’ Data Utilizing Grid 
Cell Modeling Techniques" under the direction of Dr. David M. klumpar, of 
UT-Dallas as Principal Invest! gator , and Drs. Jerry l. kisabeth (’Jniv. of 
Alberta) and Walter J. Heikkila (UT-Dallas) as co-investigators, was to be 
performed during the period from 1 September, 19’9 through 31 August, 1981. 
This original proposal had the following broad objective: 
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To apply a modeling procedure to the vector MAGSAT data In order to 
separate the terrestrial component from that due to extraterrestrial 
sources. The proposed study would contribute to the MAGSAT program 
objectives In two ways: 

1) Provide removal of those contributions to the measured field 
that are undesirable for studies of the main core field and 
localized crustal variations. 

11) Provide detailed accurate vector measurements of the field due 
to Ionospheric and magnetospheric currents for the purpose of 
studying these current systems. 

The proposed study was an ambitious one, Involving 3.6 scientific man- 
years of effort at a p-oposed cost of S196,000 In January, 1979 dollars. 

As originally proposed this effort had two major parts. The first part of 
this proposed Investigation was to model the effects of the various exter- 
nal current systems (Including Induced currents In the earth) at obser- 
vation positions relevant to the MAGSAT satellite orbit In order to 
ascertain whether or not the fields due to external sources may be removed 
from the MAGSAT data with any degree of accuracy. The second part of that 
proposed Investigation was the application of the resulting modeling tech- 
niques to the MAGSAT data by taking into account the state (temporal and 
spatial variations) of the external current systems at the time the data 
was recorded. The Investigation would be a natural extension of magnetic 
field modeling techniques under development for a number of years at the 
University of Alberta, Edmonton, Canada, and the Texas A. and M. University 
at College Station, Texas (see Kisabeth, 1975, Kisabeth and Rostoker, 

1977). Those modeling techniques have been successfully used in the 
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modeHng of magnetic field measurements from many sources Including ground- 
based magnetometer arrays In northern Canada and In Scandinavia (Oldenburg, 
1976; Kisabeth and Rostoker, 1977; Hughes and Rostoker, 1979; Bannister and 
Gough, 1978), high latitude total field variations (aB) from the OGO 2, 4, 
and 6 satellites, and vecto»* field component measurements from polar 
orbiting satellites ISIS-2 and TRIAD. 

In Its original form, the concept of the Investigation Included 
extensive detailed modeling as well as direct application to the reduction 
of MAGSAT data. 

The first part of the proposed research would deal primarily with 
the development of modeling programs necessary to predict magnetic field 
perturbations at MAGSAT orbital altitudes due to external current systems 
(and Induced current:- In the earth). These external current systems may be 
subdivided Into the following categories for study: 

1. Sq current system and the equatorial electrojet. 

11. Auroral zone and polar cap current systems along 
with lonospherlc-magnetospherlc coupling currents 
(e.g.. Region 1 and 2 flelc. aligned currents, 
polar electrojets, currents associated with magne- 
tospherlc convection). 

111. Ring current (asymmetries In the ring current could 
also be Included In category (ii)). 

1v. Bow shock, magnetosheath, magnetopause and magne- 
totail currents. 

The effects of induction In the conducting earth due to time 
variations of the external systems could be treated by assuming the earth 
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to be Infinitely conducting below a given depth. This mc'l.-i has worked 
especially well for substonn current systems (see Kisabeth, 1975). A more 
realistic conductivity structure would be extremely difficult to Include 
with sources as complex as the current systems given above (Mareschal and 
Kisabeth, 1977). 

After suitable computer programs had been developed for current 
systems In each of the categories above, a detailed study was proposed to 
determine the magnetic signatures of each current system along the MA6SAT 
orbit for varying degrees of magnetic activity. 

The second part of this proposed research was the application of the 
results of the modeling study to MAGSAT data reduction, the main goal being 
to see If contributions from the various external current systems can be 
systematically Identified and thus removed from the magnetometer data. The 
degree to which this can be done would provide valuable Information con- 
cerning the accuracy of the reduced data set which Is to be used to study 
the main field and crustal anomalies. 

The UT-Dallas proposal was provisionally accepted on July 5, 1979, for 
negotiation as one of the investigations to be performed using data from 
the MAGSAT mission. 

C. The Descoping Phase 

During the seven months following provisional acceptance, the 
UT-Dallas proposal proceeded through the so-called negotiation stage, 
during which a draft statement of work was prepared and circulated. Also, 
during this 1-year period following proposal submission one of the proposed 
co-investigators (J. L. Kisabeth) assumed a position with Industry outside 
of the space science research stream. This necessitated his withdrawal 
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from the project as a direct active participant In Its development. 

Because of his acknowledged expertise In the field and the recognized value 
of his advice we sought and received from his new employer consent to uti- 
lize his expertise as an Informal consultant If only on a very limited 
basis. 

Shortly thereafter UT-Dallas was Informed that there had been severe 
cutbacks In Investigatory funding and that the UT-Dallas proposed effort 
was being reduced to that of a “feasibility study". By the end of January, 
1980, a new statement of work had been drafted by NASA and accepted by the 
P.I. It represented only a small fraction of the level of effort that was 
deemed necessary for a full treatment ov the Important problem of external 
cur'-ent effects on low-altitude magnetic measurements. Still the delays 
continued. Active negotiation of the contract was held In limbo between 
February and August of 1980. During these seven very frustrating months 
conflicting reports were received that only some or even none of the Magsat 
Investigations might receive funding. Finally In September of 1980; two 
years after the Initial f^agsat Investigations Announcement of Opportuni ty; 
eleven months after the MAGSAT launch; and 3 months after satellite reentry 
the UT-Dallas contract was funded. The final contract called for 0.7 
scientific man-years of effort In comparison to the proposed 3.6 man-years 
and was funded at $64,000, a reduction by more than two- thirds of that Ini- 
tially proposed. 

0. Final Statement of Work 

The Statement of Work consists of “he objective to be satisfied by the 
effort, the approach to be used In satisfying the objective, and the 
tasks regulredf to satisfy It. For Magsat Investigation V-013, the 
Statement cf 'Work Is as follows: 



Objective ; 

The objective of this effort Is to perform a feasibility study of a 
modeling procedure to Magsat data to separate the terrestrial com- 
ponent from that due to extraterrestrial sources. 

Approach ; 

The Investigator shall perform a feasibility study of modeling 
programs necessary to predict magnetic field perturbations at Magsat 
orbital altitudes due to external current systems and Induced currents 
In the Earth. These external current systems will ce subdivided Into 
the following categories for study; 

a. Sq current systems and the equatorial electrojet 

b. Auroral zone and polar cap current systems along with 
lorospheric-magnetospheric coupling curents (e.g., Region 1 aid 2 
field-aligned currents, polar electrojects, currents associated with 
magnetospheric convection) 

c. Ring current (asymmetries In the ring current will also be 
Included In category b. 

d. Bow shock, magnetoshf.ath, magnetopause and magnetotall currents. 
The effects of Induction In the conducting Earth due to time 

variations of the external systems shall be treated by assuming the 
Earth to be infinitely conducting below a given depth. 

Because the software developed under this contract may be Important to 
the ability of other Magsat Investigators In separating external 
fields from anomaly fields, all results. Including software and docu- 
mentation, are to be made available (pre-publication) to otner 
Investigators via the Magsat Project Scientist. This Insure that the 
cooperative scientific effort described In Task d (below) will be 
fulfilled. 

Tasks ; 

The following casks «;hall be performed by the Investigator in 
fulfillment of the above objectives; 
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a. Obtain and render operable, existing software developed by Dr. 
Kisabeth. 

b. Extend the capability of software, as required, to enable 
modeling of: 

(1) The Sq current system and equatorial electrojet, and 

(2) Auroral zone and polar cap current systems along with 
ionospheric-magnetospheric coupling currents. 

c. Utilize the software to derive "typical" magnetic signatures of 
these current systems for local times pertinent to the Magsat orbit. 

d. Compare these signatures to signatures derived from Magsat data. 
This should be carried out cooperatively with other Magsat investi- 
gators (e.g., the U.S. Geological Survey, Dr. Pocemra, Or. Burrows). 

e. Prepare and submit to NASA periodic progress reports and a 
detailed final report documenting the results of this investigation. 
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III. ACCOMPLISHMENTS (RESULTS) 


The overall accomplishment of this Investigation has been the develop- 
ment of a new modeling technique that accurately determines the vector 
magnetic field that would arise from electrical currents distributed 
throughout the near-earth space environment. The modeling procedure has 
been applied specifically to the terrestrial high latitude Ionospheric 
currents and to the Ionosphere-magnetosphere coupling currents to ascertain 
their contributions to the magnetic fields measured by high latitude 
ground-based magnetometer arrays and to the field perturbations observed 
from low altitude polar orbiting satellites. The principles of the 
modeling technique are general and may readily be applied to other distri- 
buted current systems found In the Ionosphere and In the magnetosphere. 

The mod'*'' was used to derive typical magnetic s1gn;tures of these high 
latitude currents and was applied specifically to Magsat orbits for '‘Irect 
comparison with empirically determined magnetic perturbations derived from 
the Magsat data. To support these comparisons considerable effort was also 
expended under this contract to develop software for the reduction of 
Magsat data co produce the empirical magnetic perturbations. The following 
paragraphs describe, In detail, the accomplishments under the two c:te- 
gorles of Magsat data reduction and Field Modeling. 

A. Magsat Data Reduction 

Although the primary goal of this contract »<as to det^elop field 
modeling techniques for the near-earth magnetic field arising from external 
currents, such development cafinot successfully be carried out without con- 
comitant study and analysis of the actual Magsat Data. The following 
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paragraphs describe the Magsat data reduction efforts carried out at 
UT-Dallas under this contract. 

The Initial Magsat data reduction program was developed to read the 
Magsat Chronicle format data tapes on a U.T.D. POP 11/45 computer. The 
capability to read the Chronicle tapes and printout either the orbital data 
alone or both orbital and magnetic field-values from both scalar and vector 
magnetometers for any specified time period contained on the source tape 
was the first step to accessing the Magsat data. The program computes 
geodetic longitude and latitude, and altitude of the spacecraft and outputs 
this Information along with Inertial and magnetic coordinate positions. 

The magnetic field observations during each second are scanned with maxi- 
mum, minimum, and average values for each scalar head and each vector com- 
ponent being printed at one second Intervals. This software package formed 
the basis of additional data reduction software that stores selected por- 
tions of Magsat data on magnetic disk. Additional data reduction software 
was then developed to access the stored data, subtract a spherical harmonic 
model field, and plot each of the three resultant magnetic component 
deviations and the scalar field difference on a high resolution Interactive 
vector graphics terminal. 

The graphics capability f^r visualizing the Magsat ouservatlons as 
deviations with respect to a spherical harmonic main field model Is 
Illustrated In the accompanying figures. Figure 1 shows from top to bottom 
tne deviations 1n the North, East, Vertical components, and the scalar 
magnitude deviation with respect to the MGST(6/80) spherical harmonic field 
model * the Magsat orbit over the northern hemisphere from 10:51 - 11:21 
UT on November ''3. 1979. The polar dial on the right hand side of the 
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figure depicts the location of the satellite in an Invariant 
Latitude/Magnetic Local Time coordinate system for reference to the data 
plot. Perturbations in all three vector components, but particularly in 
thr E-W component, are evident as the satellite passes thorough the 
Birkeland current system near the dusk and dawn locations of the auroral 
oval centered near 60“ invariant latitude. This plot shows the large-scale 
features over a 30-minute portion of the orbit with a plotted resolution of 
two-seconds along the satellite orbit. In order to examine details and 
^iner scale perturbations, the graphics routine developed at UTD has the 
capability to "home in" on a selected time interval and plot the obser- 
vations at higher resolution. This capability is illustrated in figure 2 
where now only a five-minute portion of the data of figure 1 are presented 
over that segment of the orbit that crosses the dawn sector of the 
Birkeland curent system. In this plot the resolution has been increased by 
a factor of two to one-second along the orbit. A number of small-scale 
features are now visible on this plot that were unresolved in Figure 1. 

It is sometimes useful to view the magnetic perturbations as they would 
appear in an alternate coordinate system. Figure 3 illustrates the effect 
of a coordinate transformation to a system in which the field vector is 
resolved into components along (S) and perpendicul ar (D) to the earth-sun 
line. In this coordinate system the polar cap top hat in the sunward- 
antisunward component of the magnetic field, often attributed to magne- 
tospheric convection, stands out quite clearly as a positive perturbation 
in the S component at invariant latitudes above 69“ . Note in this figure 
that the 0 component perturbation at these high latitudes essentially 
zero, in sharp contrast to figure 1 in which the polar top hat contributed 
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about equally to the N and E components. 

The desire to visualize, with respect to the auroral oval, the Magsat 
derived magnetic perturbations along the orbit path has prompted us to 
develop as a first step, an orbit plotting program. Ultimately this 
program will be able to display horizontal magnetic perturbation vectors at 
regular intervals along the satellite orbit. The entire sequence of Magsat 
orbits in magnetic local time, invariant latitude coordinates over northern 
latitudes (^50°) is shown in Figure 4 for November 4, 1979. Note that as 
universal time advances through the day beginning in the upper left plot 
the orbital path of the satellite progresses across the polar ionosphere 
from the dayside of the polar cap to the nightside. Those orbits occurring 
during the first few hours of the day pass through or close to the so- 
called magnetospheric cusp where there are unique currents. At later uni- 
versal times (near 1300 UT) the satellite cuts lie along or near the 
0600-1800 magnetic meridian where the high latitude current pattern differs 
from that near the cusp. Later the orbit plane passes on the nightside of 
the 0600-1800 meridian and then returns to the dayside before midnight. 

Thus due to the location of the Magsat orbit and to the offset of the 
Magnetic pole relative to the geographic pole there is a systematic and 
repeatable 24-hour periodicity in the location of the satellite orbit with 
respect to the auroral and polar currents. 

It follows that there will also be a systematic and repeatable 24-hour 
pattern in the magnetic perturbation signatures derived from the Magsat 
data. To demonstrate this repeatability we have selected a sequence of six 
orbits over the northern high latitudes on each of three different days. 
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Figure 5 Illustrates the remarkable similarity between the Magsat orbital 
paths for similar universal times on each of the three days. The sunward 
components of the relative magnetic perturbations for each of the 18 passes 
(two missing) are shown In Figure 6 with the same relative positions as In 
the previous figure. The day to day similarity of the perturbations on 
each of the three days during Identical universal time periods Is apparent 
as Is the consistent change In the character of the per-turbatlon signature 
with Increasing universal time during each day. 

During the early hours of the day (02-08 UT) the sunward component pertur- 
bations appear to be highly structured with a tendency for a positive per- 
turbation over the center of the polar cap. At later times as the 
satellite orbit approaches the dawn-dusk magnetic meridian a nearly 
constant positive top hat develops with steep negative perturbations on 
either side of the polar cap. The dawn-dusk component of the magnetic 
field Is shown for the same orbits In Figure 7. Again there Is a systema- 
tic variation with universal time. Near the 0600-1800 meridian the pertur- 
bation signature Is negative on the afternoon side of the oolar cap and 
positive on the morning side. At later times both sides of the polar oval 
display negative D-component perturbations. 

In addition to addressing the repeatabll Ity and overall constancy of 
the polar external current systems, the significance of these periodic and 
repeatable patterns lies In their potential effect on crustal anomaly stu- 
dies. An example will Illustrate the problem. Consider a particular 
region of the earth's surface centered at some Intermediate geograph c 
latitude, say at N. The magnetic field attributed to crustal anomalies 
within this region Is determined f^om magnetic perturbations deduced f-om 
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Magsat observations taken when the satellite passes over the surface loca- 
tion In question. All Magsat passes over the location In question will 
occur at one or the other of two particular universal times, namely when 
the earth rotates through the orbital plane of the satellite. Thus each 
point on the ground has magnetic values associated with It that are taken 
at only two specific and unique universal times. Since at each particular 
universal time there Is a consistent and unique magnetic perturbation asso- 
ciated with the external currents, the magnetic signature attributed to 
each surface location will have a consistent component due ".o the external 
high latitude currents as well as to any crustal anomaly that might be pre- 
sent. Thus great caution must be exercised when Inferring magnetic anoma- 
lies at middle and high latitudes. Simply removing temporal variations or 
picking magnetically quiet days will not remove the consistent and 
repeatable external current effects discussed here. 

F lel d Model Ing 

It was originally envisioned that this research project to model the 
magnetic fields of certain distributed currents In space would utilize the 
techniques developed by J. L. Kisabeth and would extend the capabilities of 
that method as required. During the early stages of the project that tech- 
nique was evaluated as to Its applicability to the problem at hand with due 
consideration of the available resources. It was found that, although the 
technique Is a sound one. It nevertheless carried with It a number of 
constraints which limited Its scope and made It less desl’-able than It 
first appeared. One severe limitation contained In that modeling technique 
was the existence of artificial upper and lower latitude cutoffs for the 
Input current system. Currents existing high In the polar cap could not be 
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conveniently handled by the Kisabeth model. As It was our desire to model 
the fields due to distributed currents with little or no artificial 
constraints on their geometrical distribution, this limitation was con- 
sidered to be undesirable. Correcting this situation would have required 
considerable modification to the existing software. A second limitation 
lay In the resolution within which currents could be specified. The 
Kisabeth grid cell model existing In 1980 could handle currents specified 
In a grid cell having the dimensions of only 2° of latitude by 15° of 
longitude. It was rnn«<.ifred desirable to work with higher resolution than 
this since real ionospheric current variations on a smaller spatial scale 
than 2° by 15° were expected. This limitation was not considered to be 
Insurmountable. It would have been possible to expand Che number of grid 
cells to gain higher resolution, however the price would have entailed an 
Inordinate Increase In the required computational resources. Already ti.e 
existing programs required a large computer and a relatively large com- 
puting budget. A computer run costing several hundred dollars would have 
been necessary each time a set of kernals was needed to be generated. 
Because of the exploratory nature of the proposed investigation our limited 
computer budget would have been quickly depleted. 

In order to bypass the limitations discussed above we began to consider 
alternative techniques. We were driven by the deslf'e to develop an effi- 
cient but accurate technique that would not seriously restrict the form or 
magnitude of the Input current distribution. Furthermore, any new method 
developed would have to be capable of operating on a small computer such 
our POP 11/45 computer, for which no time-base operating charges are 
assessed. As various alternatives were reviewed. It became clear that a 
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relatively straight forward approach utilizing the most basic laws of 
magnetism and vector mathematics might be best suited to the available 
resources. 

The magnetic field computation technique thus developed "s based upoA 
the additive properties of vector fields. In general, the field 'ector at 
a point In space Is the vector sum of the vector components arising from 
all of the elemental field sources In the universe. In the present r.*:. . 
the magnetic field at a point Is computed by summing the contributions of 
all of the assumed current^ that exist everywhere In space. The assumed 
current distribution Is modeled by decomposing the actual current distribu- 
tion Into an arbitrary number of finite length cu-rent elements. The tech- 
nique Itself relies upon the use of an analytical expresrion for the 
magnetic field of a straight current carrying filament having an extended 
and smoothly varying cross-sectional cur-ent density. The cross-sectional 
current density profile looks somewhat like a square wave pulse with 
rounded corners. The us^ of such a platykurtic distribution ha.* been found 
to eliminate discontinuities that exist in a square wave representation and 
allows for easy calculation of the vector magnetic field at any point in 
the world space. 

The total current distribution to be calculated Is represented by an 
arbitrary number of these finite length current elements. Typically 
several hundred such current elements are used to represent the horizontal 
and field-aligned current distribution over the high latitude Ionosphere. 

By a suitable summation of the field at eacn point due to the contributions 
from all current elements, the magnetic field may be calculated anywhere, 
such as on the earth's surface or along a satellite orbit. The resultant 
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magnetic perturbations for each vector-field component are displayed on a 
high resolution vector graphics terminal by means of a computer program 
designed to allow the operator to interactively modify the model parame- 
ters. The initial development of this model was restricted to a hypotheti- 
cal satellite orbit at 90° inclination in the dawn dusk meridian plane and 
the initial computations included only the sunward component of the pertur- 
bation field at various constant altitudes. A sample plot of the output 
from this initial model is shown in figure 8. The main plot on this figure 
shows the sunward directed component of the magnetic field calculated at 94 
separate observation points as a function of latitude along a hypothetical 
orbit at a constant altitude of 500 km. The input current system is a 
"classical" large-scale Birkeland sheet current model with downward 
directed currents in the high latitude postmidnight and the low latitude 
pre-midnight portions and upward directed currents in the high latitude 
pre-midnight and low latitude postmidnight sectors. This current system is 
represented computationally in this example by 324 linear current elements 
as described above. The field-aligned currents are closed by N-S currents 
in the ionosphere at 110 km altitude. No E-W ionospheric currents have 
been included. The center of the current sheets are located in the figure 
by the vertical lines at 70.5® and 77.5® latitude. The latitudinal 
distribution of current intensity is plotted along the top of the main 
panel. The clock dial in the lower right-hand quadrant depicts the 
satel 1 i te orbi t. 

Following this initial development the capabilities of the model were 
extended considerably to allow the computation of all three vector com- 
ponents of the magnetic field arising from an assumed current system. The 
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restriction on possible satellite orbits was also almost entirely elimi- 
nated allowing the magnetic field components to be computed for virtually 
any satellite orbit over a range of Inclinations and altitudes and having 
an arbitrary angle between the orbital plane and the earth-sun line. 

Figures 9 through 12 show sample plots of the model magnetic field output 
for a satellite at 50Q km altitude with various orbital Inclination angles. 
The main plots on these figures show three orthogonal components of the 
magnetic field each calculated at 94 separate observation points as a func- 
tion of latitude. The solid curve Is the sunward component, the dot-dash 
curve Is the dawn-dusk component and the dashed curve the vertical com- 
ponent. The Input current system Is a "classical" large-scale Birkeland 
sheet current model with downward directed currents In the high latitude 
postmidnight and the low latitude pre-midnight portions and upward directed 
currents In the high latitude pre-midnight and low latitude postmidnight 
sectors. This current system is represented computationally In this 
example by 324 linear current elements as described above. The field- 
aligned currents are closed by t-S currents in the ionosphere at 110 km 
altitude. No E-W ionospheric currents have been included. The centers of 
the current sheets are located in the figures by the vertical lines at 61“ 
and 70“ invariant latitude. The latitudinal distribution of current inten- 
sity is plotted along the top of the main panel. The clock dial in the 
lower right-hand quadrant again depicts the satellite orbit. 

Two additional extensions of the magnetic field modeling software were 
then made that significantly extended the ability of the model to handle 
realistic current systems and display the results more realistically. The 
capability for ionospheric current closure in the east-west direction was 

27 


Finure 


5UNWRRD B-MELD OF BIRKELRND CURRENTS 



HLTITUOC 500KH 
INCLINRTION 0 
CUnnENT .400 HMP5 


laioo 


' \ '• • i : 

v \\\'' ''!// 

>' ''i' 

/ / ' • » \ \ '• 

* • • * 


06100 


OOlOO 


*n 5. 
•D© 

o 

50 « 

o ^ 

5 S 

p 


CLl 20.0 CL2 29.0 CL3 10.0 FR .000 lU 10.0 10 20.0 DSf .UOO 




Figure 10 


5UNWRRD B-FICLD OF BIRKELRND CURRENTS 



I2i00 


RLTITUOE SOOKH 
INCLINRTION -ID 
CURHCNT .HOO HHP5 



o p 

7, - 

; > 
X f- 

<c T 
t- 

nr * 

•"» . 


CLl 20.0 CL2 29.0 Ct3 10. 0 fn .000 lU 10. 0 ID 20.0 OST .900 


OOiOO 


Figure 11 


SUNWRRD B-FIELD OF BIRKELRND CURRENTS 



1B<00 


I2t00 
• < 




RLTITUDC 500RH 
INCL1NRTION -13 

cunnENT .uoo hhps 


IBiOO 


06 1 00 


''/V; . I . 


OOtOO 


CLl 20.0 CL2 29.0 Ct3 10. 0 FR .000 lU 10. 0 10 20.0 05T .000 


Ti 5 

TJ S 
O Z 
O s 

73 r- 

•o 

^ tS 





built Into the model. Furthermore, the resulting magnetic field pertur- 
bations could additionally be computed and plotted In an N,E, V coordinate 
system. 

Figure 13 Is a schematic representation 1n the magnetic local time- 
magnetic latitude coordinate system of the linear current elements used to 
model, to zeroth order, the naturally occurring field-aligned currents 
above the high latitude Ionosphere. Each circle with Its appropriate 
cross-hatching represents the location, current Intensity, and current flow 
direction for a linear current element. The superposition of all these 
current elements approximates a "classical" large-scale Birkeland sheet 
current model with downward directed currents In the high latitude postmid- 
night and the low latitude pre-midnight portions and upward directed 
currents In the high latitude pre-midnight and low latitude post-midnight 
sectors. For this particular Instance low-level distributed Inward field- 
aligned currents exist between 0800 and 1600 hours on the dayside and simi- 
lar outward directed currents appear on the nightside. These distributed 
currents are necessary to maintain continuity of the horizontal Ionospheric 
closure currents shown In Figure 14. In this figure the arrows represent 
the direction and relative magnitude of the Hall and Pedersen Ionospheric 
closure currents. For the current system Illustrated here the majority of 
the Birkeland current closure occurs 1n the N-S direction with the eastward 
and westward closure currents becoming proportionally stronger near the 
dusk and dawn sectors respectively and decreasing to zero near noon and 
midnight. 

Figures 15 and 16 depict the magnetic perturbations that would be 
observed by magnetometers on a satellite along two hypothetical orbits that 
pass through the modeled current system at an altitude of 450 km. In each 
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figure the upper panel shows the latitudinal profile of the three vector 
components of B in an X,Y,Z coordinate system while the center panel 
depicts the field components In the more conventional N,E,V coordinate 
system. Shown In the bottom panel of Figures 15 and 16 are the field- 
aligned current densities encountered at each point along the satellite 
orbit. The clock dial In the upper right quadrant depicts the satellite 
orbit In a Magnetic Local-Time and Invariant Latitude coordinate system. 

In Figure 15 the satellite has an orbital Inclination of 90° and Its orbi- 
tal plane Is contained In the dawn-dusk meridian. As expected the major 
perturbation In the magnetic field appears In the east-west component a*d 
has Its greatest gradient co-loca«ted with the local field-aligned current. 
In Figure 16 a slightly different orbit has been chosen with an orbital 
Inclination of 96°. Now the satellite passes slightly to the dayside of 
the dawn-dusk meridian. Comparison with Figure 15 reveals that the east- 
west magnetic component Is virtually unchanged whereas a substantial N-S 
component has now developed In the magnetic field. This kind of comparison 
Illustrates the strong effect that relatively small displacements In the 
location of the measuring point can have upon the vector maqne^l* field. 

At this stage the model had reached a level of development whereby 
extensive testing of Its predictive capabilities could begin. Pursuant to 
that end, we took the Initial steps to conduct comparative modeling of an 
agreed-upon Input Ionospheric and field-aligned current system In coopera- 
tion with the National Research Council of Canada group. A second test of 
one portion of our new modeling package was carried out using a horizontal 
Ionospheric current system published by Akasofu et al (1981). These 
researchers have deduced a model of the total Ionospheric current distribu- 



t1on based upon five minute averages of the magnetic fle’d measured on the 
ground with the Alaska Meridian Chain of Magnetometers. Using a likeness 
of their Ionospheric current distribution from Figure 4 of the referenced 
paper, we have calculated the three-dimensicnal field-aligned current 
distribution around the polar 1onosphe**e that would be requlreo to maintain 
continuity of the total current system. Our results are shown In Figure 
17. This figure Is a color-coded representation of the direction and 
magnitu.'e of the resulting field-aligned current densities over the polar 
region. The color hues of the red-orange-yellow portion of spectrum repre- 
sent varying Intensities of downward directed currents and those of green 
and blue represent various Intensities of upward directed currents as shown 
In the bar scale*" In the lower right hand corner. The gross features of 
this plot closely resemble those deduced by Akasofu et al . and shown In 
their Figure 5. This pattern 1r. also generally consistent with the overall 
empirical distribution of Region 1 and Region 2 field-aligned currents 
deduced from satellite data. Our model has also been used to p'*-oduce the 
magnetic field perturbations that would be observed from a setelllte 
passing through this current system, as well as those that would be 
observed on the ground below. Those latter distributions should compare 
with the Input magnetic field measurements used by Akasofu. 

A further test of our modeling capability Is, at the time of this 
writing, being undertaken. It serves as an example of our ability to model 
complex Input current distributions. This test utilizes as Input to the 
model a complex Ionospheric current distribution deduced by Y. Kamide of 
Kyoto Sangyo University In Kyoto, Japan from an extensive set of ground- 
based magnetic field observations. The ionospheric current distribution 
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used as Input to the model Is shown in Figure 18. This figure illustrates 
the degree of complexity with which our model is capable of operating. 

With these highly spatial varying horizontal currents and the requirement 
that the currents be continuous, the field-aligned current distribution 
required for closure is shown in Figure 19. Combining these ionospheric 
and field-aligned currents together we model the magnetic perturbations 
that would be observed at a satellite crossing over the current system as 
shown in Figure 20. The satellite orbit is shown ir. the polar dial at the 
right side of the figure. 

One further capability of our modeling procedure is illustrated in 
Figure 21. Owing to the high computational efficiency and the com- 
putational organization of the modeling technique we are not limited to 
calculating the magnetic perturbations at just a few points or just along a 
particular satellite orbit. The model easily calculates the vector com- 
ponent perturbations everywhere. This figure shows in a three-dimensional 
perspective drawing the relative amplitude of one component of the magnetic 
perturbations at equally spaced grid points everywhere on the surface of a 
spherical cap at 500 km altitude over the north polar regions down to 50“ 
latitude. 

C. Comparisons of Model with Magsat Perturbation Signatures 

Our modeling capability now allows us to make direct comparisons bet- 
ween model predictions and actual Magsat perturbations by additional soft- 
ware that allows us to calculate along an an actual Magsat orbit the 
magnetic perturbations that would be seen at the satellite for an assumed 
current distribution. This gives us the capability for direct comparison 
between measured and predicted perturbation fields. By successive iteration 
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of the input current system the current distribution that yields the best 
fit between the measured and predicted magnetic perturbations can be deter- 
mined. 

An illustration of th? effec" of large-scale ionospheric and field- 
aligned (Birkeland) Currents on tne Magsat data is shown in Figure 22. The 
left side of the figure represents model magnetic perturbations calculated 
from a relatively simple distribution of north-south and east-west 
ionospheric currents that are fed by classical Region 1 and Region 2 
Birkeland currents. The right hand side shows the actual data from Magsat 
as it passed over the northern polar regions between 10:45 and 11:06 UT on 
December 4, 1979. The large scale features in all three compoents (North, 
East, and Vertical) measured at Magsat are reproduced in the model calcula- 
tions. Note that the relative current strengths have not been optimized in 
the model so that the scales cannot be directly compared. 
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IV. SOFTWARE DESCRIPTIONS 


A. MAGSAT Chronicle Tape Manipulation Programs 

Software described In this section was developed for or modified for 
use In reducing the Magsat data at the University of Texas at Dallas and 
graphically displaying the reduced data for further analysis. 

1 . Tape Reading Programs 

MSMDI - Magsat Monitor to Disk 

Intended by the original author Ron Cook as a monitor to enable 
several user selected modes of operation within the data analysis 
subroutines, and now serves only to tell the user the present date and 
time and then to call MSDI. This routine was written by Ron Cook and 
modified by Dale Greer, both of UT-Dallas. 

MSDI - Magsat Data to Disk 

This program Is of major Importance to Magsat data reduction. It 
reformats Magsat data and stores It In a disk file for subsequent 
plotting. 

In addition to selecting a time window, the user must select either 
the "AL" or the "OR" mode. In the "AL" mode, the disk file Is called 
MAGSAT.DAT and will contain all orbital observations within tr,e time 
window followed by all the vector observations within the time window. 
In the "OR" mode, the disk file will be called SATORB.DAT and will 
contain orbital observations only. 

Finally, In the "AL" mode, the user must select either regular den- 
sity, whereby each disk file vector record will contain tf’e average of 
16 vector observations, or hexadecuple density whereby each vector 
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observation 1s written to disk. 


Each record contains five four byte words. Integers written to 

disk are double precision. 

Development of this program was made more difficult than necessary 

by Irregularities In the tape format and by Incomplete documentation 

of same. The following Is a list of surprises and Irritations. 

(1) There Is a four minute overlap at both ends of each orbital 
record, but there Is no mention of this In the CSC manual. 

(2) Orbital records begin on odd numbered hours plus 56 minutes. This 
makes day boundaries somewhat difficult to cross, especially when 
the manual says nothing about It. 

(3) Some vector records begin at negative times, no docu":6ntat1on, of 
course. 

(4) Some vector records are types 5, 6, 7 while others are types 
8,9,10. The manual mentions 5,6,7 only. 

(5) Some expected vecto'" records are missing. This Is mentioned In 
the manual but It would be better to fill In missing data with 
dummy data Indicating that these data are missing. Orbital Inter- 
polation goes haywire when the observation time Increment jumps 
from a small value to up to one complete revolution. 

(6) Some scalar records are missing Independently of the presence of 
vector records. This Is also mentioned In the manual but It Is 
bothersome to work around. 

This software was authored by Ron Cook and was extensively modified by 
Dale Greer, both of UT-Dallas. 

RACTW - Request and Convert Time Window 
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Called by MSDI this program requests the time window In year-month 
day, hour-minute-second and converts It to modified Julian day and 
millisecond of day. 


Author-Ron Cook 


M^STAT - Magsat Tape Status 

Called by MSDI to report and act upon errors encountered while reading 
Magsat tapes. 

Author-Ron Cook 


SBFIBM-Swap Bytes from IBM 

Called by MSDI to convert from IBM to DEC. single precision Integer 
format. 


Author-Ron Cook 


CVTIBM-Convert IBM 

Called by MSDI to convert from IBM to Dec. singie precision floating 
point foimiat. No source file for this one. 

Author-Lou Wadel 


2. Data Plotting Programs 
ORBPLT-Orbit Plot 

From the data 1n SATORB , DAT, ORBPLT plots the orbital section 
within 50 to 90 degrees noruh or south latitude, depending on user 
selection. In geomagnetic coordinates for each orbit In the file. 
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This program Is used to find the start and stop times to be entered 
Into MSPLT so that each vector data plot will cover the same area. 

Author-Dale Greer 


MSPLT-Nagsat Plot 

From the data In Magsat, MSPLT plots delta-B, the difference between 
the field measured and the field model (currently MGST680), the scalar dif- 
ference Is derived from the vectors. 

MSPLT also plots the orbital section In geomagnetic coordinates for 
the time during which the measurements were taken. The orbit plot Indica- 
tes northern or southern hemisphere through the drawing a solid or a dotted 
curve for each respectively . 

Since It takes about a half second to process each vector record and 
to find the values of the four points to be plotted therefrom, the user Is 
told how many points may be plotted from the chosen Interval, and asked how 
many of these points are actually to be plotted. Points not plotted are 
simply skipped over and no averaging takes place since this is done to save 
time. 

Author-Dale Greer 

MtLDG-Fleld Generator 

This program generates the field components from a model for sub- 
sequent subtraction from the measured values. 

Author-unknown 

Modified by Dale Greer 


TOD-TIme of Day 

Called by MSPLT to get hour-minute- second of day from millisecond 
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of day. 


Author-Ron Cook 

2.1 Subroutines Called by the Plotting Programs Only 
POSIT - Position 

This subroutine Is the main part In a package of three Interdepen- 
dent subroutines obtained from some U.S. Government source. 

POSIT was Initially written to get the five orbital data points 
necessary for Interpolation, send these to STIROB to accomplish 
such Interpolation, then send these results to SATPOS to get the 
Interpolated coordinates Into latitude, longitude, and radius. 

Now, POSIT not only does that but also conditions coordinates .jr 
geomagnetic coordinate Interpolation, flags northern or so«w.ierii 
hemisphere, and crosses day boundaries. 

Author-Unknown 

Extensively Modified by Dale Greer 

STIROB - ? 

This program does the actya'' Interpolation. 

Author-Unknown 


SATPOS-Satel 1 Ite Position 

SATPOS takes X, Y, Z coordinates from oqsIT and converts to lati- 
tude, longitude, and radius. It can accomodate a rotating or a 
stationary coordinat? system. 

Author-Unknown 
Modified by Dale Greer 
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TIME 


Converts hcur-mlnute-second to minisecond of day. 

Author-Dale Greer 

2.2 Subroutines Called by both the Tape Reading and Data Plotting Programs 

YMDDOC-Year-month-day; day of century 

Converts year-month-day to day of century. 

Author-Ron Cook 

DOCYMD-Day of Century; year-month-day 

Converts day of century to year-month-day. 

Author-Ron Cook 


B. Modeling and Model Plotting Programs 

The software described in this section was developed at the University 
of Texas at Dallas to model the magnetic fie'^ds that arise from distributed 
electrical currents flowing thorugh the near earth space environment. Care 
was taken to permit the calculation to be carried out with the fewest 
number of restrictions placed on the distribution of the input current 
system. All programs in this section were authored by Dale Greer of the 
University of Texas at Dallas. 

CURDIS-Current Distribution 

This program defines the model current distribution. 

The model comprises a large number of current filaments assemblac 
to simulate the Birkeland sheet currents. The filaments are lil^e 
current carrying wires in that they have thickness, but unlike wires 
they have a smoothly varying cross-sectional current density, i.e., 
the density varies as the hyperbolic secant of the square of the 


52 



distance from the center. Such a distribution shall be referred to as 
Kurtic (from the Greek "Kurtos"-bul ging or swelling. 

The basic element of the model Is the filament. The filaments are 
combined Into loops and the loops are combined Into cells. The base 
of the cell Is In the Ionosphere and comprises one north-south fila- 
ment and one east-west filament. These filaments are sourced and 
sinked by the field-aligned filaments which are tangent to the magne- 
tic field lines at thf Ionosphere. These filaments are straight and 
are three earth radii In length. 

A more complex modtl, In which the field aligned filaments curved 
with the field lines all the way to the equatorial plane was tested. 
The Increased complexity made a barely discernible difference In the 
final result and so was scrapped in favor of the former, simpler, and 
hence faster, model. 

The current In each filament is defined at the base of the cell. 

The north-south and east-west components In the Ionosphere give the 
current for the entire corresponding loop. The thickness of each fila- 
ment is defined separately. These parameters are changed by modifica- 
tion of the source program. 

The spatial parameters are defined interactively by the user. The 
user must give: 

(1) The number of cell rings, i.e., the numbef" of cells one would 
encounter on a trip from the pole to the equator. 

(2) The number of cells per 360 degrees longitude. 

(3) The latitudinal range in which the cells are distributed. 

(4) The maximum radius of the filaments, in meters. 
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(5) The extent to which latitudinal compression takes place. 
Finally, CURDIS puts the thickness of each filament, the current in 
each filament, the endpoints of the filaments, and the angles through 
which the filaments must be rotated in a disk file called DIS.OAT. 

The endpoints and angles are not given for each filament, rather 
advantage is taken of the symmetry of the model to decrease the size 
of the data file and speed up processing. 

AMPLT-Ampere Plot 

AMPLT shows the current flowing through the surface of a sphere 
just above the ionospheric currents for the Birkeland current model 
defined by CURDIS. 

Each circle represents a field-aligned filament and encloses about 
90 per cent of the kurtically distributed current. 

Each line in one of these circles represents 10 K amp (current 
equals number of lines plus or minus 5 K amp). 

CURPLT-Current Plot 

CURPLT shows the current vectors in the ionosphere for the 
Birkeland current model defined by CURDIS. 

BRKALC-Bi rkel and Calculate 

BRKALC calls MAGMOD to find the current density and field com- 
ponents of the Birkeland current model defined by CURDIS at points on 
a ci rcular orbit. 

User must select the orbital altitude, the orbital inclination with 
the magnetic pole, the orbital angle with the dawn to dusk line, and 
the number of measurement points in this orbit. 


The user must then elect to calculate the field of all the 
currents, the field of any one component bank of currents, all of the 
north-south currents, for example, are in one component bank or the 
field of all the currents, less any one component bank. 

Finally, user must select either a polar pass, whereby the orbit 
will go from 50 degrees latitude western hemisphere to the same lati- 
tude eastern hemisphere, or an equatorial pass from 50 degress lati- 
tude north to 360 degress south in the eastern or western hemisphere. 

MAGMOD-Magnetic Modeling 

Called by BRKALC and 3DBRKC to evaluate the current density and 
field components of the Birkeland current model defined by CURDIS. 

MAGMOD takes the end points of a filament as found by CURDIS, rota- 
tes it to its proper place, and calculates the field components and 
current density at point X,Y,Z given by the calling program. 

BRKPLT-Birkeland Plot 

BRKPLT plots the vaues found the BRKALC in XYZ, NEV, or SDV coor- 
dinates and shows the orbit. 

3DBRKC- 3-D Birkeland Calculate 

3DBRKC has the same function as BRKALC except that it creates a 
data set for a three dimensional plot with measurement points being on 
a sphere of radius "ALT" greater than that of the earth's and there is 
no provision for an equatorial data set. 

3DPLT-3-D Plot 

This program was not written by, but was modified by Dale Greer. 
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In Its original form 3DPLT simply plotted a 3-D picture of an 
array. The form used here does the same but with perspective, zooming, 
and some interactive data conditioning included. 

J3RKC- J Birkeland Calculate 

JBRKC is like 3DBRKC except that it calls JBRKFN to find only the 
current density. 

JBRKFN - J Birkeland Function 

Is like MAGMOD except that it only calculates the current density. 
JBRKP - J Birkeland Plot 

JBRKP plots the findings of JBRKC as an eleven color field. It 
does this by sending direct commands to a Tektronix 4027. Eleven 
colors are derived from the eight available through the use of pat- 
terns. 

BUFUN-Buffer Function 

From JBRKP the POP 11/45 sends commands to the 4027 so fast that 
the 4027 would crash and have to be turned off it is weren't for 
BUFUN. BUFUN just does a few calculations and has no effect on the 
program but to slow it down. 
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V. RECOMMENDATIONS FOR FURTHER WORK 


This contract has demonstrated that a much more thorough analysis of the 
Magsat vector magnetic field observations when combined together with a 
versatile modeling technique of the various contributions to the Magsat 
measurements holds the promise of yielding valuable new Insights on the 
subtle Influences of space currents on main field and crustal anomaly stu- 
dies. Further analysis of the Magsat data will also contribute Tundamental 
new knowledge to our understanding of the electrodynamic coupling between 
the Ionosphere and the magnetosphere. The low level of scientific effort 
allotted to the present study and Its restriction to the demonstration of 
feasibility has not permitted the type of In depth analysis required to 
address the above problems. This contract has permitted the Initial deve- 
lopment of an Inexpensive, versatile new tool whose ultimate application Is 
still ahead. 

A number of specific recommendations can be made for further studies: 

1) Recognizing the universal time effect described 1n this report, 
analyze the space current perturbations during the specific 
quiet days that have been selected out of the Magsat data for 
anomaly and core field models. Apply the modeling procedure to 
remove the average quiet time space current contribution. 

2) Perform an In-depth analysis of how the separate parts of the 
Ionospheric and magnetospheric current systems show up in the 
vector Magsat measurements. 

3) Use the modeling procedure to analyze the external and 
ionospheric current contributions at locations and altitudes 
relevant to proposed new magnetic field missions such as the 
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Geopotential Research Mission (GRM), low altitude tether 
satellite, and other free flying shuttle launched satellites. 

4) There has been no consideration given to the altitude distribu- 
tion of the Ionospheric currents. All modeling and data reduc- 
tion techniques developed to date have assumed that the 
horizontal Ionospheric currents flow In a highly localized 
shell around 110 km altitude. Yet It Is well known that there 
Is an altitude distribution of the Ionospheric conductivity and 
furthermore that the various terms of the conductivity tensor 
have different altitude profiles. These effects will beome 
much more Important In low altitude magnetic field measuring 
satellite missions. 

5) The problem of Induction has not been treated In our model to 
any order. To the extent that currents Induced In the earth by 
currents flowing overhead Influence the magnetic field at 
Magsat, then we have not handled them. There are modifications 
that can be made to the model that will make correction for the 
effects of earth-induced currents. 

6) Inversion: With respect to modeling the Ionospheric and magne- 
tospherlc currents, one would like to be able to solve the 
Inverse problem. That Is, from the measurement of magnetic 
field, compute the responsible currents. This may not be 
possible from a single point satellite measurements, but see 
Item 7 below. 

7) Ground level observations: Observations from a single 

satellite operating alone do not provide sufficient Information 
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to solve for a unique system of currents. With some clever 
techniques to combine selected sets of single point Magsat 
measurements and by additionally incorporating ground based 
observations taken simultaneously over a large portion of the 
earth's surface we might just be able to overcome the non- 
uniqueness problems associated with single satellite measure- 
ments and find a true representation of the distributed 
external currents. 
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VI. SUMMARY AND CONCLUSIONS 


This contract has demonstrated the feasibility of modeling the magnetic 
fields that arise from distributed currents In the near-earth geospace 
environment. The modeling procedure has been applied to the high latitude 
Magsat observations to show that substantial perturbations arise In the 
Magsat vector field, after subtraction of a spherical harmonic model of the 
earth's main field, that are due to currents flowing In the earth's 
Ionosphere magnetosphere system. The contract has also Involved data 
reduction and analysis of the Magsat data with respect to the potential 
effects of Ionosphere-magnetosphere currents on the application of Magsat 
data to studies of magnetic crustal anomalies. 

SPECIFIC RESULTS 

1) Developed Interactive data analysis software to permit graphical out- 
put of three-component magnetic field perturbations relative to a model 
geomagnetic field In different coordinate systems with Interactive control 
of time base resolution. 

2) Displayed and plotted Magsat vector measurements as perturbations 
relative to the Magsat spherical harmonic model magnetic field at latitudes 
above 50® geomagnetic latitude for ’all orbits during the first two months 
of the mission. 

3) Developed a new forward modeling software procedure that determines 
the vector magnetic field due to distributed space currents. 

4) Demonstrated that 3) could be accomplished efficiently, accurately, 
and with computational economy on a small (PDP 11/45) computer system. 

5) Used the modelling procedure to determine the separate effects at 
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Magsat orbit due to the currents flowing: 

a) In the Ionosphere along the auroral oval In the E-W direction 

b) In the lonsophere across the auroral oval In the N-S direction 

c) along the magnetic field direction between the Ionosphere and the 
magnetosphere. 

6) Pointed out that periodicities of the Magsat orbit with respect to 
the auroral and Ionospheric current systems can lead to contamination of 
anomaly and core field models by space current effects. 

7) Recommended that dawn and dusk orbits be treated separately to 
evaluate the effects of 6) 
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A TECHNIQUE FOR MODELING THE MAGNETIC PERTURBATIONS PRODUCED 
BY FIELD-ALIGNED CURRENT SYSTEMS 

O.M. Klumpar and D. M. Greer (Center for Space Sciences, University 
oTiexas at Dallas, Richardson, Texas, 75080, U.S.A.) 

This paper presents results of a computational procedure that utilizes 
various assumed distributions of Ionospheric and field-aligned currents to 
model magnetic perturbations observed at high latitudes from the polar 
orbiting MAGSAT satellite. The highly sensitive vector magnetometers on 
MAGSAT repeatedly observed magnetic field perturbations on essentially 
every transit of the high latitude Ionosphere. These perturbations, with 
field components lying predominantly In the magnetic East-West direction, 
are customarily viewed as the signatures of oppositely directed paired 
sheets of electrical current flowing parallel to the geomagnetic field. 
These paired current sheets are typically regarded as being highly 
restricted In latitudinal extent and elongated In magnetic longitude. The 
model developed under this research utilizes a computationally fast and 
mathematically simple technique that allows the magnetic field of a distri- 
buted current system to be calculated oy representing such a system by an 
arbitrary number of hypothetic^’ linear current elements. The facility of 
the technique derives from tht t'«e of an analytic expression for the magne- 
tic field of a linear current element having an extended and smoothly 
varying cross sectlona'' current density; thus eliminating unwanted discon- 
tinuities. Magnetic perturbations typical of those encountered at auroral 
latitudes by MAGSAT are produced by the model when realistic current con- 
figurations are chosen. Direct comparisons between the model field pertur- 
bations and those measured by the MAGSAT magnetometers permit more refined 
models of the Birkeland currents to be developed. 


1. University of Texas at Dallas 
Richardson, Texas 75080 

2. ER 

3. N. Fulcushima, S. Matsushita 

4. - (a) Oral presentation 

Submitted to Fourth Scientific Assembly of lAGA, Edinburgh, August 3-15, 1981 
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ABSTRACT 

A METHOD OF CALCULATING MAGNETIC FIELDS 
DUE TO SYSTEMS OF DISTRIBUTED CURRENTS 

D.M. Greer 

d.M. Klumpar, (both at: Center for Space Scien- 
ces, University of Texas at Dallas, Box 688 
Richardson, Texas, 75080) 

Electrical currents in the ionosphere and in 
the magnetosphere produce large amplitude magne- 
tic field perturbations that are detected by the 
highly sensitive magnetometers on the polar 
orbiting MAGSAT satellite. THis paper describes 
a computationally fast and mathematically simple 
method that has been developed and applied to 
modeling the magnetic field produced by the 
Birkeland current system. The technique allows 
the magnetic field of a distributed curr.nt 
system to be calculated by representing such a 
system with an arbitrary number of hypothetical 
linear current elements. The facility of this 
method derives from the use of an analytical 
expression for the magnetic field of a straight 
current filamenty having an extended and 
smoothly varying current density. The cross 
sectional current density profile of such a 
filament looks somewhat like a swuare wave 
pulse, of arbitrary width, with rounded corners. 
Thus the system is free from unwanted discon- 
tinuities and the field component in any direc- 
tion and at any point in the model space is 
easily calculated. Magnetic perturbations typi- 
cal of those encountered at auroral latitudes by 
MAGSAT are produced by the model when realistic 
current configurations are chosen. 
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MODELING THE HIGH-LATITUDE MAGNETIC FIELD 
PRODUCED BY DISTRIBUTED IONOSPHERIC AND 
MAGNETOSPHERIC CURRENTS 

D. M. KLUMPAR, (Center for Space Sciences, The 
University of Texas at Dallas, Box 688, 
Richardson, Texas, 75080) 

O.M. GREER 

The magnetic field disturbances resulting from 
distributed currents In the high latitude 
Ionosphere and from the Birkeland currents, 
which extend outward Into the magnetosphere, are 
computed using a highly efficient computational 
technique that has recently been developed. 

Using this technique the magnetic effects of the 
large-scale distributed Birkeland, Hall, or 
Pedersen currents can be computed Independently 
to ascerrtain their separate contributions to 
the magnetic perturbations typically measured by 
magnetic observatory arrays on the ground or by 
satellite borne magnetometers. Use of this 
technique Illustrates the complexity of the many 
contributions from various currents that combine 
to produce the net magnetic disturbance that Is 
measured. Such model ng analysis provides the 
basis for Improved Interpretation of ground and 
satellte magnetic observations In terms of the 
responsible currents. Such Improvements will 
subsequently lead to more realistic represen- 
tations of the true horizontal and field-aligned 
current systems than Is available from the 
customary "equivalent current representation" 
and hence to a better understanding of magne- 
tospherlc dynamics. We present a set of model 
calculations of the magnetic vector components 
arising from an assumed Ionospheric and 
Birkeland curent system and compare the pre- 
dicted magnetic signature to that typically 
measured from low altitude polar orbiting 
satellites and from high latitude magnetometer 
chains. 
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Model Magnetic Field Perturoations at Magsat 
due to External Current Systems 

0. M. GREER ana 0. M. KLUMPAR (Center for 
Space Sciences, Tbe University of Texas at 
Dallas, Richardson, Texas 75080) 

Significant magnetic field perturbations 
due to currents in the ionosp'ere- 
magnetosphere system are observed on virtually 
every Magsat orbit over the hign latitude 
ionosphere. We utilize a model jf distributed 
currents consisting of the horizontal 
ionospheric currents and f ield->aligned 
(BirXeland) currents to compute the 
perturbation magnetic fields along Magsat 
orbits. Tbe computer coce models the 
distributed currents by aecomposition into a 
large number of linear, fi*<lte cross section 
current elements for which che magnetic field 
can readily be computed. The perturbation 
field at each point .‘n space due to the entire 
distributed current system is then che vector 
addition of the appropriate contributions from 
each current element in Che system. We compare 
che model derived magnetic perturbations with 
chose deduced from actual Magsat measurements 
to iteratively determine the distribution of 
ionospheric and Birxeland currents for 
particular Magsat orbits. 
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A TECHNIQUE FOR MODELING THE MAGNETIC PERTURBATIONS ORIGINAL PAGE IS 

PRODUCED BY FIELD-ALIGNED CURRENT SYSTEMS Qp pQQR QUALITY 

D.M. Klumpar and D.M. Gre«r * 

Center {or Space Sciences. The University of Texas at Dallas, Richardson. Texas T5080 


Abstract . A computational procedure is introduced for calculating 
the magnetic {telds produced by virtually any distributed system of 
electrical currents. This procedure is being applied to the modeling of 
magnetic fields produced near the earth and on its surface by horizon- 
tal currents flowinc in the ionosphere and by the so-called Birkeland 
currents flowing along the geomagnetic field at high magnetic lati- 
tudes. This report describes briefly the principles that underlie the 
technique and illustrates the results obtained when the model is applied 
to the interpretation of perturbation fields being measured by the 
polar-orbiting magnetic fields satellite (MAGSAT). Even for a very 
simple assumed current distribution we calculate magnetic field resid- 
uals whose large-scale features are similar to those deduced from 
MAGSAT measurements. A predominately sunward magnetic pertur- 
bation is obtained over the region poleward of the Region I current 
system as a natural consequence of balanced Region 1 and Region 2 
currents. The model predicts the existence of low-latitude magnetic 
effects of auroral currents that represent potential sources of error for 
spherical harmonic representations of the geomagnetic field. 

Introduction 

The magnetic field measured from near-earth orbit, although domi- 
nated by the earth's main magnetic field, contains significant com- 
ponents that arise from electrical currents flowing in the ionosphere- 
magnetosphere system . In particular, at high latitudes near the auroral 
oval, currents flowing parallel to the geomagnetic field may cause per- 
turbations in the locally measured magnetic field in excess of 1500 nT 
directed primarily transverse to the main geomagnetic field. Such 
field-aliipie6 current signatures were first measured from satellite 
I9<33-38C and reported by Zmuda et al. [196C, 1967]. Since that time 
magnetometers on a number of low-altitude satellites (TRIAD, ISIS. 
AE-C. S3-2) have been used to infer the nature of the magnetic pertur- 
bations arising from field-aligned currents [e.g,. Hjima and Potemra. 
Id76a,b; Klumpar et al.. 1976; McDiarmid et al., 1978a.b: Klumpar. 
1979; Bythrow etal,. 1980. 1981; Doyle etal.. 1981; and others). In 1979 
a dedicated magnetic fields satellite was launched to make the first 
global vector survey of the geomagnetic field. 

The magnetic fields satellite, MAGSAT. was placed in a near-earth 
sun-synchronous orbit with the objectives of making precise magnetic 
field measurements to accurately describe the earth's main magnetic 
field and to map. on a global basis, ii^ fields caused by sources in the 
earth'.s crust (Langel, 1979], It was recognized early in the program 
that the sensitive magnetometers would also measure the magnetic 
.field produced by currents flowing in the ionosphere-magnetosphere 
system external to earth and that at some locations and limes these 
external effects would even mask the crustal anomaly fields. 

Analyses of these externally caused magnetic perturbations in terms 
of the responsible currents have generally assumed a highly idealized, 
local system of paired, infinitely long, planar, parallel current sheets 
oriented perpendicular to the satellite trajectory. Klsabeth [1979] took 
a major step towards eliminating these restrictive geometrical as- 
sumptions by devising a computational technique to determine the 
magnetic perturbations that would arise from more general distribu- 
tions of currents. The present work represents a new effort to model the 
magnetic perturbations resulting from distributed electrical currents 
flowing in space around the earth. To that end a method has been devel- 
oped for calculating the magnetic field at any point in space due to an 
assumed spatial distribution of electric currents. This paper briefly 
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describes the computation technique and discusses several important 
aspects of the magnetic perturbations that result from a simple large- 
scale Birkeland and ionospheric current system resembling that pre- 
viously deduced from the large body of near-earth magnetic field 
measurements. We conclude by comparing signatures derived from 
MAGSAT measurements with those predicted by the computational 
technique. 

Several aspects of the field perturbations derived from the modeled 
large-scale current system raise questions about commonly accepted 
interpretations of satellite-borne magnetic observations. They are; 

1. For a balanced Birkeland current system in which all of the field- 
aligned current closes in the N-S direction between the Region 1 
and Resrion 2 field-aligned current sheets, there still existsasun- 
ward-directed magnetic perturbation in the region poleward of 
the Region 1 currents. Thus, contrary to some suggestions, a 
polar “top-hat” field distribution does not necessarily imply a 
net field-aligned current in the Region I current system. 

2. Significant magnetic field perturbations due to high latitude 
currents extend to latitudes well below those normally associated 
with the auroral zone. 

3. The existence of a positive perturbation in the sunward com- 
ponent of the magnetic field over the polar cap does not require 
a cross polar cap current, but rather arises as a natural conse- 
quence of a balanced classical Region 1 and Region 2 Birkeland 
current system. 

The Field Modeling Technique 

The magnetic field compuution technique is based upon the additive 
properties of vector fields. In general, the field vector atapoint in space 
is the vector sum of the vector components arising from all of the ele- 
mental field sources in the universe. In the present case, the magnetic 
field at a point is computed by summing the contributions of all of the 
assumed currents that exist everywhere in space. The assumed current 
distribution is modeled by decomposing the actual currentdistribution 
into an arbitrary number of finite length current elements. The tech- 
nique itself relies upon the use of an analytical expression for the mag- 
netic field of a straight current carrying filament having an extended 
and smoothly varying cross-sectional current density. The cross- 
sectional current density profile looks .somewhat tike a square wave 
pulse with rounded corners. The use of such a platykurtic distribution 
has been found to eliminate discontinuities that exist in asquare wave 
representation and allows for easy calculation of the vector magnetic 
field at any point in the world space. 

The total current distribution to be calculated is represented by an 
arbitrarj' number of these finite length current elements. Typically 
several hundred such current elements are used to represent the hori- 
zontal and field-aligned current distribution over the high latitude 
ionosphere. By a suitable summation of the field at each point due to 
the contributions from all current elements, the magnetic field may 
be calculated anywhere, such on the earth's surface or along a satellite 
orbit. 

As an illustration of the technique, we show in Figure 1 a simple, 
hypothetical, ionospheric current distribution that is characterized by 
dominant north-south currents. A large-scale eastward electrojet 
current flows from noon across the dusk hemisphere toward midnight 
while a westward electrojet current is directed through the dawn hemi- 
sphere from noon to midnight. .All currents are confined Co a shell 
running between 60’ and "6° latitude. This horizontal current system 
requires, for continuity, that there be accompanying field-aligned 
currents, which are shown in Figure 2. The circles represent the loca- 
tions of possible field-aligned current elements and the cross hatching 
represents the direction and relative magnitude of the field-aligned 
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distribution of ionospheric currents 

Fi^r» 1. High latitude distribution of horizontal ionospheric currents plotted 
on a latitude versus local time coordinate ffrid. All horizontal currentsarecon- 
sirained to How inside a channel between Hu* and 76* latitude 
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currents. Vertical hatchinu depicts the presence of an inward current, 
which IS ixren to exist at low latitudes in the post-noon to midniirnt 
sector and at hiirh le’.itudes in the morninK hemisphere. Horizontal 
hatching indicates the outward current at high latitudes on the evening 
hemisphere and at low latitudes in the morning hemisphere. These 
high and low latitude field-aligned currents represent the Region 1 
and Region 2 currents deduced by lijima and Potemra 1 197ba.b| from 
a study of the TRIAD data. 

In addition to the sheet-like currents discussed above the chosen 
horizontal current distribution requires tnat there be an additional 
downward current near noon and an upward current near midnight. 
These currents partially feed into the eastward and westward electro- 
jets that flow across the dawn and dusk hemispheres in Figure 1. 

The primary question we seek to answer is. what are the magnetic 
fields produced by such a current system’ Figure :}dispiay.s the "-esults 
of the computation in the upper three panels as latitude profiles of three 
components of the magnetic field that would be measured at a satellite 
at 450 km altitude moving along the dusk to dawn meridian The cart- 
wheel plot at the upper right depicts the path along which the field 
IS computed. The bottom panel displays the field-aligned current den- 
sity profile as a function of latitude along the satellite orbit, which, as 
anticipated from the previous figure, passes through only the classical 
Region 2 and Region 1 Birkeland currents. As expected, the major per- 
turbation appears in the East-West component with the steepest gradi- 
ents occurring at the location of the local field-aligned currents. 

Smaller but still significant, magnetic field contributions are found 
in the region equatorward of the low latitude termination of current 
flow The magnetic field strength at .50° latitude, a full lO' of latitude 
equatorward of the auroral currents, are of the order of 10 to 20 nT 
and decay only slowly with 'lecreasing latitude. The presence of such 
mid-latitude magnetic effects in satellite measurements may. if not 
properly attributed to the external current system, contribute toerrors 
in a proper spherical harmonic representation of the mam magnetic 
field. Concentrating now on latitudes poleward of the high latitude 
currents, it is apparent that there is again a significant magnetic per- 
turbation due to the modeleil currents. This .so-called "polar top-nat ‘ 
field perturbation is directed primarily sunward. This moilel snows 
that It arises as a natural consequence of a balanced current system in 
w nich the upsvard and downw ard currents along a meridian are equal 


Such level shifts observed in satellite data have in the past been inter- 
preted as evidence for a net field-aligned current in the Region 1 system 
(Sugiuraand Potemra 19T5|. or as a result of cross polar cap currents 
[Fujii etal.. 19811. Finally, we note that the modeled currents also pro- 
duce a notable perturbation in the vertical component of B. Such an 
effect has been detected in the .MAGSAT data. 
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In Figure 4 are shown, for the same current distribution, the mag- 
netic field profiles along a somewhat different orbit where the satellite 
passes on the davside of the pole. The main features of the magnetic 
profiles observed in the dusk-dawn meridian are preserved w'lth the 
primao' difference being a reduction in amplitude of the N-S com- 
ponent and a widening of tne E-W profile, as the satellite makes a mure 
oblique pass through the current system. 

.Although only profiles at satellite altitude have been shown, the 
modeling procedure described here also allows the field components 
to be calculated on the surface of the earth. Such a model will permit 
further understanding of the external sources of the magnetic fields 
measured on the ground and in space and. in particular, of the complex 
magnetospheric-ionosphere electrical circuit. 

Comparison with .MAGSAT Data 

From the vector magnetic measurements made by MAGSAT it is 
possible to derive a difference field by subtracting a model represen- 
tation of the earth's main magnetic field from the measured field. This 
difference field is presumably the resultant perturbation that arises 
from the combined effects of externally produced fields due to currents 
in space and induced in the earth, crustal anomalies and inaccuracies 
in the spherical harmonic model representation of the core field. If in 
the first approximation we choose to Ignore the latter two contribu- 
tions to the difference field because they are small, and assume a steady 
state external current system, the difference field will represent only 
the effects of external currents. Figure 5 shows such a difference field 
for a dusk to dawn MAGSAT pass over the northern latitudes on 
Novemoer 13. 1379. In producing this difference field a thirteenth 
degree and order spherical harmonic representation of the main mag- 
netic field referred to as the .MGST 16/8O1 model (Langel et al.. 19sO) 
was used. For this orbit the satellite passes just to the dayside of the 
dusk to dawn meridian. The largest deviations from the model, up to 
950 nT in the East-West component, occur as the satellite passes over 
the dusk auroral oval between 10:56 and 10:59 LT. A deviation of 
approximately one-third of the E-W component is also present in the 
N-S component. Furthermore the N-S component has a ISO nT resid- 
ual deviation extending down below 41° geographic latitude The 
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Figure 5. .Measured difference fields along a .M.AiiS.AT orbit on .S'ovemoer 13. 
1979 .All three components are plotted relative to the MUST lO boi field model 
of Langel et al. 1 19Hil| 


absence of such a low-latitude residual in the E-W' component is some- 
what at variance with the predictions of the model discussed in the 
previous .section. Two possible explanations may account for this vari- 
ance. The first is that the auroral current model di.'cussed in the pre- 
vious section may not accurately portray the real currents during this 
pass, and that the real current system is producing no E-W field com- 
ponent at low latitude. .A second possibility is that the 13th order 
spherical harmonic expanison has been contaminated by the low lati- 
tude fields due to polar currents and as a result has these effects built- 
in as a part of the main geomagnetic field. The actual resolution of the 
discrepancy may rest in some combination of these two possibilities 
and will be one of the objectives of further modeling efforts. 

Further comparisons of this .MAGS.AT difference plot with the per- 
turbations calculated from the simple model and shown in Figure 4 
reveal gross similarities in the large-scale features and substantial 
differences in details. The latter arise from small-scale variations in 
the actual current sytem that were present during the .MAGSAT pass 
shown in Figure 5 for -.vhich no attempt to model r.as been made in the 
current distribution discussed here. This comparison serves to illus- 
trate the complexities that exist in the real Birkeiano current s.vstem. 

Conclusion 

A general calculational procedure has been developed to compute 
the magnetic field perturbations arising from distributed ionospheric 
and lonosphere-magnetosp.here coupling currents .A simplified cur- 
rent distribution nas been chosen to illustrate the technique and the 
resulting magnetic perturbations have been compared to actual mag- 
netic perturbations measured from M.AGSAT The balanced Birxeiand 
current system produces non-negligible low-latitude magnetic field 
perturbations .A sunward magnetic perturbation is also produced at 
latitudes poleward of the high latitude current sheet. 
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CONTINU 


C 

c •*** 

c 

22000 


22R30 

C 

^ * * * * 

22200 

i **** 

C 

C **** 

z **** 

c 

22431 


22432 


C 

c **«« 

c 

22433 


C 

^ ** *< 

24900 


DISK write vector DATE AND TIME EILTIR 


Ff JVMJD.LT.lBMjpiGO TO 24999 
r JvMjD,6t.|EMJD)G0 to 90000 
r JVMJO.NE.lBMjoiCO to 22030 
F JVMSOO.Lt.lBMSOplOO TO 24999 
Fi JVMjo.CT.If MJDJGO TO 22200 
F(JVM3ni^.Gt.IEMS00)60 ro 90000 


SOMETIMES THE FIRST OBSERVATION TIME IS NEOATIVEl 
IF IJVMSOD'.LT’.R) JVMSOObJVM3O04S64 00000 
WRITE MAGNETIC DATA ONTO DISK 

5 E««§ 5 ’ 

IF ri JVMSOO-JVM9P2)‘,LT,0) GO TO 2202 
If n JVMSOO-JVMSOal.Ot. (AV0ELT*R0Nf9S3.)) 

GO to 24900 

WRITE REAL DATA TO DISK 

NMEAS«NMEA5*1 

NMaNMEASflORB 

WRITEC4»NM) NMEAS, JVMSOD, (RVSTAT(3, J), J-1.31 
GO to 24920 


WRITE DUMMY DATA TO DISK 
NMEASaNMEASf < 

NMiNMEASilORS 


24901 

1 

24902 


1 VM90?wJVHSD2FRV0ELTwRQN 
F I JVMSQ2 .6 E. 80400000) J VMSD2W JV^^SO?- 
F (Jvm9o6-jOmsd2).lt,0) go To 2420I 
F f JVMSOD-JVMSOE) ,C e, ( ftvOELt*RON)) 

0 0 24910 

GO to 24902 

IF ( ( JVM8DD-jVM8D2fB6400000) ,LC>. (RVDELTwRON)) 
GO 0 24910 
RVStAT(3, 1 )«99999, 


88400000 




Ln 

Lo 


FOR 
mSd 

R247 

0248 

0249 

0250 

Sin 

0253 

0254 


0255 

0258 
0257 
0255 

0259 

0260 
0281 


028 

0264 

0260 

0266 

0287 

0268 

CO 

0269 

0278 

0271 


M IV. PLUS V02- 
TN 




10104109 


1O.OCT.02 


PABE 6 


RVSTAT 

!R?TE(i^ftM 
CO to 94000 
24910 NMEASa^MFAS.l 

NMaNM.l 

24920 JVMSOPaJVMSOO 

24999 CONTINUE 


■999||* 

*NMMS* JVMS02, (RVSTATO, J), Jal,31 


Ol 


60 TO 12000 

^0000 TYPE 00101 

ACCtr'T 0001 1, NRWN 
IF Inrhn.gt,!) go T 
CALL_REwlNf lu.TSJ 
NRITE(4I I )N0,NMEAS, 


90010 


09990 


\\:p 




wn, 


I ) 


00001 FORMAT 
00011 FORMAt 


TO 00010 
KK,NH, lORB 
on Tb mil 

FORMAT statements 




00100 format 

00101 FORMAT 

00120 forma! 

• 7*. • I 

forma 


00130 
00200 FORMA 
06903 FO 
C 


40A1 ) 

* I0MAGSAT OISKt) 
•HOrm FpR_0R8II_RECgRD 


)RMA 


iH 


, 

ENTER TFOR HEXAOE 

UNIPENTIFIEO record TYFEl 



return for regular • ) 


0272 


END 


POOR QUALU 



s 


0O-- 

p» 


FORTRAN IV-PLUS V 
RACTW.Ff*^ /TR| 


/HR 


19104197 


19-0CT-8? 


PAO€ 1 


0001 

0002 

0003 

0004 

0009 

0009 


0007 

0008 

0009 

0010 

2®h 

0012 

0013 

^0t4 


00 

00 

00 


0ig 


0010 

0020 

0021 

0022 

0023 

0024 
0029 
0029 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 
0039 

0036 

0037 

0038 

0039 

0040 




RACTH • 1980 APR 18 R.H.COOK 

PROGRAM REQUF9T9 aAO CONVERTS TIME NINOON. 


SUBROUTINE RACTh 

lugicalai AtNsrn 

INTEGERa4 I0OOC,! 
REAL*4 0IN9(1) 


BMJD, IBHOO, IB90D, lEOOCf lEMJO, TCH0O» IE900 


COMMON /M0N/MO9EO. AIN8 
IBVOC, 


COMMON /90TG/ 


BHOO, 

EVOC, 

EHOO, 


BMOYtlBDOM, 

BM0H,fS90M, 


CCCCC 
1 0000 

c 

1001 1 


EQUIVALENCE 


[FMOH 
(AIN 9 (I)fBTN 9 ( 1 )) 


BQOHiIBDOVflBOOCflBM J3. 


I? 


ON*IEOOY*1EDOC»1FMJO* 

00 




10100 


iKiifnBffiiillllffl''**- 


10100 MOSEOal 

TYPt 00100 
GO TO 10011 

OErOOE (13.00002. AT NS, ERRal 01 00)1 yor.TMOY. lOOM, 1 SHOD, TBMOH, I B90M 
CALL YMAoDCdYOCf IM0V,ID0M,l00Yf ioOCflOOH) 


1 0101 


YOC 

MOY 

DOM 

DOY 

DOC 

004 


BYOCa 
JBMOYa 

I BOOMa 
BOOYa 
BDUCc 
iBOOWa 
IBMJDa 

Ibsuob 
Type 00109 

10200 M0St0a2 

TYPt 00200 
GO TO 10011 
OErOOE(13 “ 


1 0201 


BDOCfl9020 

6H5D*3600AlBMQH*60tIB9OM 
109« IBDOy. IBMJD* IB90D 

200 
0011 

ErOt)E(l3,00002, AIMS, ERRal 0200)1 YOC f I MOY, 
ALL YM006 c( lYOC. IMOV, lOOM, 100Y,l00t, lOOH 
EVOCal YOC 


RS/.KRR!l0?00nyOC,IMgY, |DOM, lEHOD, IEMQH, IC90M 


EMOYatMOY 
^DOMafoOM 
[OOYajoOY 
iDOC«|LOC 
lEOOHafoOM 

lisioallMSoAlJSiSlEMOHAejflESOM 

TYPE 00209, IEDOV, iNJD, lESOD 





Ui 


OF FCOR 




FORTRAN IV-PLUS V0 
RACTW.FTN 


RR 41 


0042 

0043 

0044 
0049 

0046 

0047 

0046 


/TRiAlI 
90000 RFTURN 


/WR 


19104197 


19-OCT»02 


PAOE 2 



9TATEHFNT8<K><P<»'4><><»<><WK><>fi 


0000L 
00100 FORMAl 

“IRMAI 

)RMA1 
■OR^'A' 
.-><>* 
END 


001 
001 
00209 


AND Tl'tei 


US"??, 


» loOY •,t3,'» EMJD I, 9,T, E800 T,I8) 


HHMM8SI ) 
IH89I ) 


Ln 


oo 

cn 


OKiG?NAL F/Xu - 
OF POOR QUALITY 



ISIO9I07 


15«0CT»S2 


p«oe< t 


KfS;!iE 


/WR 


00 

CTl 


00RI 

0002 

0003 

3004 

0009 

000ft 

0007 
0000 
0009 
00)0 
00)1 
001 ? 
0013 
00)4 

00 5 

001 6 
0017 
00)0 

0019 

0020 
0091 
0022 
E023 
0024 
0029 
0026 
0027 
0026 

0029 

0030 

0031 

0032 

0033 

0034 

0039 

0036 

0037 
0036 

0039 

0040 

0041 

0042 

0043 

0044 
0049 

0046 

0047 


C 

90000 

91000 


91010 

91020 

91030 

91040 

91090 

01060 

91070 

91000 

91090 

91100 

91110 

91120 

99000 

C 

09000 

09010 

09020 

09030 

09040 

09090 

09060 

09070 

09060 

09090 

09)00 

09)10 

09)20 

r 


SUBROUTINE MT9TAT ( RUN, IS) 

• FORb'0', 

0)GO TO 99000 

09000) afor, run, is 


CO TO (9)010.91020.91030,91040.91090,91060,91070,91060,91090, 


CO 
NRT 
GO 
WRI 

nri 
60 
NRI 
60 
NRI 
00 
NRI 
GO 
NRI 
GO 
NRI 
GO 
NRI 
GO 
NRI 
RtTURN 


ifms.eq. 

CO TO (9)010.91020.91 

‘SJifc’in 

NF«0 
NRa-1 

CAlL SKIPR(RUN.NF,NR# 18) 
CO TO 90000 
NRltE(06.09020)AFOR 
0 99000 

E (06,09030)AFOR 
0 99000 

F(06.09040)AFOR 
REnInIRUN, IS) 

0 99000 

E (06.09090) IFOR 
0 99000 

E (06.09060)AFOR 
0 99jt00 

E(06.09070)AFO0 
0 99000 

E (06, 09060) AFOR 
0 99000 

E(06,09090)AFOR 
0 99000 

E(06,09100)AFOR 
0 99000 

E(06.091 10)AFOR 
0 99000 

E(06r 09120)AFOR 


FORMA 

FORMA 


» 'DEVIC 
,32X, 'P 

FORMA 

i|_4 


FORMA 


.32X, 'P 

forma 


»32X, 'E 

FORMA 

i‘ A 

»32X, 1L 

FORMA 

1 ' A 

»3?X, 'E 

FORMA 

r A 


FORMA 

1 ^ 

*3|x. •( 

FORMA 


»32X, 't 

FORMA 

1 ' 

*3|X, IE 

FORMA 

i5 

»32X, If 

FORMA) 

A 

32X 'F 


: 1|I2a* status INDICATOR ',12) 
IRMY ERROR') 

•OF' ) 

)A0 EXCEEDS REQUESTED LENGTH') 


lECOlil 
fOD' ) 


POINT') 


)T«) 
tlTE 

rAOR| 00 


[RRQRl 

-RRORf 


ATTEMPTED 
ODD BTT| 
LECT 
Tnva 


ELECT ERROR') 




UNDO 

RECO 




PARAMETER' ) 


ENTEO* 
92000 


BYTES' ) 


0040 


END 


iO 




m 

000? 

0003 

0004 

0005 
000S 
0007 
0000 
0009 


RAN IV< 
RM rtN 


PLUS V 


(0100192 


1 0«OCr*B2 


I!?S;WE/hr 

SURRnuTiNE 9RFTBM(BA0,BA1 ,NB) 
SWITCH BYTES PRHM IBM DATA INPUT 

?s»ut .»» js?prs6s.?i 

L0CICAL*1 BA0M),BA1M),B 

DO 000D9 iBl.NB.I,? 

■BA0(I) 


00099 


ColtlNU^ 

BETURN 

eRd 


PABC 1 



IBI06I5B 


ts«ocT»a:> 


PABC< 1 


A PORTRAM IV-PLUS V02-B1E 
O, SHPlHM.ptK /TR|AlL/WR 


00RI 


000 ? 

0003 

0004 

0005 
0000 
000 ? 
0000 

0009 


SURROUTIRE 8HPTRM(HA0,HA1 


SWITCH halfwords PROM IBM 1 
NHi number of ha 

INTEGER*? HARM) *HAt ri ) f H 

00 00109 I«l*NH«t*P 
H«HA0n ) 


INPUT 

LFWOR08 INPUT 


00109 

C 


HAi (l)«HA0rifl) 

HAjJU'.I-H 

CONTINUE 


RETURN 
END 




U1 

v£> 


Si J Tw.'NiiSiHO 




fortran 

orbrlt.f 


IV. PLUS V09-O1I! 15107104 IS-OCT.B? 

TN /TR|AlL/HR 


PA8E I 


0001 

000 ? 

0003 

0004 

0005 
0005 

0007 

0008 
0009 

iSIt 

0012 

0013 

0014 

0015 


00 

S 00 
00 


0020 

0021 

0022 

0023 

0024 

0025 
0025 
0027 
0025 

0029 

0030 

0031 

0032 

0033 

0034 

0035 
00 3 8 


0037 


C**************** *************•***•***. 

c 

c orbplt plots data in satorb.dat 


o> 

o 


c 

99B 

999 

C 


COMMON /ORBIT/ XX, YV. 

htstini ???{«!/?} 


ZZr SLAT, SLONC. salt MO A Y,RHEM 
EX(4), DARRAYfl4), PNUMf?) 


*4 ILAT, MIT 

CER*5 IM0SM2),1PDM(12) 

0|R*4 MMJD, fBYMD* IBMJO, IBDOC* IBHMS, 1B4SQ0 

e s:i ’'"***'’ 


. - JC ( HU t 

R*4 NUPO, iMftY 
IMDS. IPBM 


DATA IMOS/060, 031. 059, 090, 120, 151, 1 SI Ml 2, 243, 273, 304, 3 

1 • 1 0 ' / 

DATA PI / 3,14159265 / 

DATA ONE / 1.0 / 


34/ 

31/ 


PILE MUMBERiNN •) 


WRITE (S*99S) 

FORMAT ( I ENTER 

REAP (6.9?9> FNIJM 

FORMAt (20 

pARRAYI l)*»l» 

OARRAYi 2 -ENUMri ) 
pARHAYl 3 oFNUm|2J 
OARRAYC 4)k0 

OPFN(UNIT«2,NAHEbOARRAY.TYPEpIOLOI f 
1 FORM* I unform* TTEO t .ACCESS- 'B!RECT»,RtC0R0SIZE»5) 

8IJB liiji ss«Mr 


8 |ie 

READ 


reymd, i 

ROREFT, ^ 
rhms, RTO 


s?su. 

OS, c& 




Iembod 


:bm 8 
'em' 

DUM3 


RBHMB ■ RHMS 
CALL TIME fRMSOO, 
TS ■ RMSOO 
N ■ 4 Y 2*N0 
READ (2'N) RHMS. 
REHMS ■ RHMS 
NOM? ■ ?*N0 . 7 


RHMS, RTOS) 
RTOS, CO 


Cl 1 1 1 M 1 1 1 1 1 1 1 1 1 II 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 II 1 1 ) 
^ **•* REQUEST AND ACCEPT USER SELECTED PARAMETERS 
1111! CONTINUE 


OUlGIf'J.ViL Pr.GE IS 

OF POOR QUAUl'Y 




UD 

O 


15lH7lfl4 


15-OCT-B? 


e(MR 

0039 

0040 


0041 

0042 

0043 

0044 

0045 

0046 

0047 

0046 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 
0056 

0059 

0060 
0061 

006? 

0063 

0064 

0065 

0066 
0067 

0088 

0069 

0070 

0071 
007? 

0073 

0074 

0075 

0076 

0077 
0070 


ID4Y ■ 0 


NR 

01000 PO 


C 

01001 


ITP (6«010001 RBYHD. RBHM6, RPYMO. RE 

RHAT ( • THIS DATA Hl| CONTAINS MAOIAT 6r0I 
1 I YVfiMDO HHMRSB.T 4J0 •/ 

? • FROM i ,P7.0,2X,F6 .T»Ix, 15/ 

3 • TO •,F7.0,?X,F0,1,4X,I5) 




0100 ?) 


READ 

0100? formal 
C 

IF (lEX.EO.n 

c 


fR 1 
IPX 


TO EXITI RCTURR TO CONTINUE 


00 TO 20909 


02000 


02001 

r 




forma 




|h|m 


FOR NO 
FOR 80 


HERN 

HERN 


HEMISRH 

hImisph 


1/ 

) 


IF (IHfM-EO,0) them ■ 1 
rhem b iAem , 

IF (IHFM ,0T. 1) RHEM 


• -1 


02010 F0Rm1t^(*Pcn?^R plot 


READ (A 


CA 


Flh!)DBlD5cil^e?0 
IF (lYOc.EO,-- - 


YM 


.02030) lYO 
DDOC flYOC 


C 

020 ) 1 
C 


C 

02020 

C 

02021 


plot start 

OC/IMOV.IDO 

,1M0Y,160M, 


' ) 


START OATEiYYMMDD 
- DOM 

IOOY*IDOC»ID04) 


NRITE (6, 
FORMAT 
READ (6 


(RSHMS 

9 ■ 0 . 


0) FTMJ0«1BMJD 

T START TIMElHHMMSS 


6*02011) 

( T ENtpR PLQ 
,02029) RSHMS 


' ) 


,EO’.0’. ) RSHMS a RRHMS 


IF 

R T 0«* " D a 
IDAY B fImJD • IRMJD 
T0, RSHMS, 


RTOS) 

1400000. 

ra ♦ 60I00. 


CALL TIME (TB, RSHMS. RTOS 
T0 • TB ♦ rF)MJD«IBMjD)*66 
IF (RSHMS, EQ.RBHMJ) TB ■ T 

NRITE (6*020201 

FORMAT ( » enter PLOT END OATElTYMMOO I ) 
READ (6,02030) IVOC, IMOY.IDOM 
CALL YMDDOC (IYOC*1 MOY,!OOM,IOOY,IDOC»TDOn) 
,AMjoBioor 


If 

NRITE 


.0 

(lYOC 


Read ( 

02029 rORMAf 

02030 FORMAT 
C 


OC416020 

,EO,0) LAMJ0«IEMJD 

'"•nsw 


5eaD*(6!02059J 


ER PLOT 

rnhPs 


END TIMElHHMMSS • ) 


PA6E< 2 

HMS, lEMJO 
T DATA •/ 


• ) 


r '>! ;? 8 ; 5 ; E . 


| 5 IH 7 |(I 4 


I 5 * 0 CT*B? 


If (RNHMS.E 
CALL TIME*( 


.^ 9 .®.) RNHM 3 ■ REHMS 
(Tf, RNHMSi RTOS) 
iLAMJD.IBHJD)*R 640 R 0 R 0 . 
.EQ.REHMB) TC ■ TE - I 2 le 00 , 


“ TE ♦ aAMJO-IBMJD)*Bl 
IP (RNHM9, EQ.REHMB) TC ■ 

DLT ■ 70000. 

NOPn B IMTf {TC-TR)/DLT^1) 

I B 0 
13 ■ B 
{4 B 13 
IPL B 0 
Ipm B R , 

PACF B 0 . 


INDICATE USER BELECTCO PARAMETERS 


09 R 00 X B 0 . 

y B 0 , 


AOl B PAOE t t 

f :ll{ ?ltt ES 
tit 5f 

ALL CALMP T 0 .. 0 ., 0 * 3 ) 


PAUS 

CALCMP(K,Y, 1000,2) 


1 0010 


100?0 


10000 N B 1.2 

IP fN.EQ.?) cn TO 1RS10 

lDOCBnMJfi , ls 020 

PMJObMmJD 

RHMB B rShMB 
GO To .10070 
!DOCbLAMJd»I S020 

RMJOBI.AMJO 
RHMB B RNHHB 


MJO,0,,S) 


HX B 2 ^ ♦ (N. 1 )* 3 ,B 

CALL symbol (HX.,?,,14,0HVyMMOD,0,,B) 
call number. (HX.A,,. l4,RVMD«0.,»l) 

MX B MX ♦ I, 

CALL symbol IHX..?..14,BH MJO.0.,S) 

HX B MX ♦ -14 

CALL NJMBeA rHX,e«,,14,RMJ0,Bt,«l) 

HX B MX ♦ ,B6 

CALL SYMBOL. ( HX , , , 1 4 , BHHHMMBS , 0 , , « ) 
IP fRHMS.LT. 100000 .) HX B hX ♦ ,14 
If (RHMB. tt. 10000 ,) HX B MX t -|4 
JF (RKMS.Lt.l 000 -) HX B MX ♦ ,I 4 


If fRMMS.LT . 10 . ) HX, B HX ♦ 
!aLL NJmAcR (Hfi. 0 .,. 14 ,RHM 8 


• — It 

. 100 .) HX B HX 


T B I • 

All* 

M 8 . 0 .,< 



SSSJt-rFU"'-"’ 

^pesn CONTINUE 


15107104 


IP (IHEM .OT. J) 60 TO B501I 

CALL OYHbOl rO.fO.P, .28, tdHNORTHERN HEMISPHEOC » * I 9) 
60 TO 0501] 

1010 CALL SYMBAl (8..0.7,.2Sf 19H80UTHE0N HCMl SPHERE » 0. » 1 9) 
>01 CALL SYMBOL 10. ,0.4, •I4,BH0RB1T Iy.O.,B). 

CAlC symbol I 2.,0,4».14,4HfROM.O,.4) 
call symbol rp.S,A.4*.H.2HT0,l.,2) 



is^iiDH 

call SYMBA. /999.,999,,.14,TEXT,B.,2) 

if 355S?lt 

V\ 00 TO 20040 

IppPPPPPPPPPfAPPPPPPpPPPPPPAPPPPPPPPPPPPPPPPPPPAPPPPPPPPPPPi 

PLOT DATA 


20000 I ■ 


20001 


I TR • TS 

* ♦ 1 

M ■ fM ♦ OLT 
CHK ■ TM, ♦ JS • IOAY*6e4BfO0B, 
P (TCHK ,LT. BS400OB0.) 60 TO 
DAY • rOAY A 1 

i !l|hS;«!lfSt;{2S:lifS:e5;|j! 

P*"W4-n) iio|o!!ASo2O?A003O * 


28001 


PM,NE,0] 
$ r 0 0 3 A 


) 00 TO 20070 
) 60 TO 20900 
,sIdtm,tI,noh2) 


20020 

20030 

^0040 


20041 


r.n TO 20070 


|^L!i^-A!°A0o)0!2^0B0?200O0 

tJ^WjPL^^trj I2) 00 TO 200 

iri;!! I Si ;r iti'faft 

60 to obAoo 


SO to obAoo 

(QPG p f l3PL»l«4P((13PL-l)/4nP3,5 ♦ 1,78 
rORf. P 9,5-TNt(flAf t*l>^**J*5«T* 

;ALL CAL'mP (XdRO#YORO,0,S5 


Of Focf: v 




to 


7MR 


f BI07IR4 


t9«ocr»B? 


PASCi 9 


0174 


73 

J? 

73 

79 

30 

U 

03 

34 


a 89 

a 33 

a 37 


IS 

33 

li 

33 


fl(34 


39 
33 
37 

II 

0230 

0231 
0202 

0203 

0204 


0209 

0209 

0207 

0208 

0209 

0210 
021 I 
02 2 
02 3 
0214 
0219 
0213 
0217 
0213 


C 

20090 




i371 


. :2>x»/937i: 

[li 

■■ f-Qi.',, ,H, . 

CALCMP ()f,Y,l,l) 


0309. 


C 

20060 


C •*** 

C 

20070 


21000 

C 

c *•** 

c 


I PH 

[F (IP 


TM ♦ T9 • IOAV*36409000, 


0 20000 


IPM ■ 0 
CALL CALCMP 
GO TO 20000 


fx»v»0.n 


DRAH LATITUDE CIRCLES 

00 


21000 R ■ I0«40|10 

5*2"^.09JnU*pLi80‘ )*91M(JAP1/ 


B0.) 

/9i.) 

. 1 ) 


CONTINUE 

indicate orbit plot 


START AND STOP TIH£ 


DAY ■ IBMJp 
CALL NUMBER 
XTIM ■ .1.9 

n!r homL! 

CALL 3YMeOL 


A. lOAY ^ 

(.9B,l.l9,.14,DAV,0.,.n 


(TERT,TMB) 

(XTIM,VT1M, .14,TEXT»0.»6) 


YTIM ■ g 

Call hohJIe 
call 9YMR0L 




c •«** 

f 


DRAM MlT LINES 


IM.YTIH, ,14,TEXT,f .«6) 
AND INDICATE MLT 


22010 



V 


) 

.) 


O' 


> s 

n 33 

TlS 
D % 
C5 S 
X r- 

■O Z 

C * ‘ 




PAGE B 


) 


• ) 

O O 

"n :n 

"n '} 

t ; . • 

C' 3 

r- 

i' : 1 

,» • 
^ 

Ci 


oi 


IBIU7I54 


1B«0CT*6? 


FORTRAH 

ORBkr.f 




subroutine HOH] 


CTEXT,TM) 


BIBBO 


ENCODE 

FORMAT 

RETURN 

END 


cn 


( 


V 


PABE 9 



I 


cn 




15I0BI81 


1B«0CT*B? 


PAOe< 1 


0001 

0002 

0003 

0004 
0009 
0006 
0007 
0000 
0000 
001 


!ir 


001 
0012 


00 

00 

00 


00i 0 


00 

00 

00 


0020 

0021 


•••• MSPLT PLOTS DATA IN HA68AT.0AT 



ON P09 
ON *P 


;3),CO(3),Vf4),B(4),7MA0(4) 

!B0).V^L(4,60),7LOr4),B0OV(4) 

[i)2ATExrB)«DARRAY(14),FNUM(2) 


INUC IR*4 
NT '0 R*4 
NT .0 ;R*4 
InT :c|r* 4 
REA ,*B TBi 



t.NPR.NLR.N 


NSOD 

N90D 


JTB 


c^MfeoN IboJ.ikpfe^” 

COMMON /ORBIT/ XV, YY.ZZ. SLAT. 0LONB, salt. IDAV 
DATA IMDB/000,031 , 0 BD, 060 , 120,1 B 1,1b 1,2 12, 243, 279, 904, 334/ 


DATA OARRAY / jM , I A j , i 6 | , * B • , < A I , < f < ^ < . • , * 0 < , < A I , < T < , 

DATA MGBT / •Mi|l6l.lil,l|?,f6‘,'0<,‘i«,<,<,<D<,<A<,»T»,0 / 


DATA PI / 3,14199208 / 


0022 

0023 

0024 

002s 

0020 

0027 


0020 

00; 
00: 

0031 


129 
0030 


003^ 

0033 

0034 


0039 

0030 

0097 

0090 

0039 


*••• DATA ENTRY SECTION 


990 

999 


NRITE (6,990) 

FORMAT ( » ENTFR file NUMBERiNN • ) 
READ |b,999) FNUm 


^jrmai 


(ll)lFNUMjjj 


IRRAY 
DARRA Y ( ! 3)«fnum 


0PFN(UNITa2,NAMEB0 ARRAY, TYPEbIOLQI , 
FoBm«TunFORHATTED> ,ACCEB 80<D|RECTl,RECORO9IZEaS) 
READ (211) N0,NMFA9, JDUM.NM, ioRB 


read 

READ 

READ 

READ 


r^u # ■ a f a 

RBYMd,1@00 
revmd,! ■“ 
rorefT,*, 

JDUM, JM9 


• • • * 


TSaJM^O 
FIND time of first 



D,CO 


H 
H 
UM| 


► llN2,DlJM3 


AND LAST MA0NCTIC DATA RECORD 


NaSflORB 

READ ( 2 Tn) JDUM, JMSD, (V(I), Ial,3) 
CA|.L TOO (JM90 ,FHmS) 

READ (2'N) JDUM,jH9D»(V(I),lal,3) 


O' 


OF FOUR i^OAlITY 


C' 


'O 


fortran 

NSRLT.r 

R040 

0041 
004? 


0043 

0044 


0045 

0046 

0047 
0040 

0049 


0050 

0051 
0?5? 

0053 

0054 

0055 

0056 

0057 
0050 

0059 

0050 

0061 

005? 

0053 

0054 

0055 
0055 
0057 
0050 
0059 

0070 

0071 

007? 

0073 

0074 

0075 

0076 

0077 
0070 


I V*PLUS 




S|!lE/NR 

CALL TOO (JMSO.LHNS) 


16100101 


lS-OCT-0? 


PAOE- ? 


Mill CONTINUE 

OPE NfUNI Tal .NAME«M68T, TYPE* '0101 f 

iFOftM.TffiRMif TEOT 1 


WRITE (6* 

01000 FORMAT ( 

i ; FROM 
3 * TO 

0110? FORMST^f'^^SfJ 


000) RBYMD«FHM3. IBNJD.REYMO.LHMSf I EMJ 

TH15 DATA File contains haSsXt 5&ta t 

YYMMdO HhMMSS.T MJD */ 


{.F7:0,?X,Ffl, 
».F7.0,?*,FB 


.i«4X ,15) 


R 1 TO EXIT! RETURN TO CONTINUE •) 


READ 

01103 FORMAf \l) 

IF (lEX.EO'.n GO TO ?0909 

ft*** determine MSOD for beginning and ending plot POINTS 


C 

01101 

C 


ot 

D 


01 100 miwvvm. P. 

READ (6,01109) FIMJ 

plot 
hRS 


NR 
FO 

READ 


Jll/l'f'JSU 

' (5,01100) 


R , 
RSH' 


START MODIFIED JULIAN OAYlDDODD • 


START TIMEtHHMMSS < ) 


01110 FoSmSt ’ f *?*E»*fJ 


C 

01111 


I 

READ (6,01109) 


R PLOT END modified JULIAN OAYlODDDO 
LAHJD 


, rsiJSi'M'iJU R plot end TIMEiHHMMSS • ) 
READ (5.01106) RNHMS 
0110B PORMAf (F7.0) 

0li09 FORMAt (I5I 


IDAY^F* 
Iday?* 


BMJO 

HMS 




T 

R* 


mjd-ibmjd 
day 

ME (TB(RSHMS,RT0B) 
DAY*B5400000, 


IF (RSHMS.EQ.FHMS) TB*TB*3000. 

.AMJDalEMJD 


IF (LAMJD.EQ.0.) lAMJDalEMJI 

If (rnhms.eo'0,) rnhmsbihms 

CALL TIME (T?,AnHmS,RT 0S) 
TEaTE*(LAMJOrlBMJD)* 864 00000. 

IF (RNHMS.FQ.LMMS) TEaTE-l?0S0i, 


WRITE (6,01200) 
01200 FORMAT ( • ENTER 


I FOR SDVI 2 FOR NEV • ) 




VO 

00 


H IV-PLUS 

ptn 




0079 

0000 

0001 


0082 

0003 

0004 
0009 
0080 
0007 
0080 

0089 

0090 
009 | 

0092 

0093 

0094 


0099 

0090 


0097 

0090 

0099 

0100 
0 01 


LL/W 0 

ICO 

IF (icn.t&.ei ICOal 


19108101 


I 9 > 0 CT *02 


PACE* 3 


READ ( 0 c 9 )? 01 ) 
01201 FpRMAt M) 


* * • * 
* * * * 


DFTERMINE TIME RETREEN DATA POI< 4 T 3 f 'lUMJER 
IN INTERNAL AND RECORD NUMBERS FOR FIRST A 
DATA POINTS TO RE PLOTTED 


OF POINTS 
AND last 


Na 94 lOR 0 

READ ( 21 N) JDUM. JTl « (Vm, I 
NbNM 
READ 
T 


■ 1 * 3 ) 


lEAO ( 2 IN) JOUM.JTL* (y(l),I»l, 3 ) 
^L»JTL 4 ( IEMJD»t 6 mj 5 ) * 00400000 . 



)/DLT) 


**•* 

***• 


UO)*S 0400000 .)/DLT) 

determine number of points to 0 ^ PLOTTED. FRACTION TO OEi 
SKIPPED* AND AVERAGE TIME SETREEN PLOT POINTS 


NRITE ( 0 * 01300 ) NUPO 
01900 FORMAT ( • there ARE l*^* 
1 t DATA POINTS WITHIN ThE 


01301 

C 


J I 

- J., 5 , 3 » 


JllMAV Q 


OF POI 




irfExnii^. 


0 T«(T 


NaNuPO /4 

R.NfR* 1 )/FL 0 AT(N) 


0 
0 

0104 
0109 
0106 
0 


0 

0 

0 

0 

0 

0 

0 

01 


0100 

09 

10 

U 

9 


c ••** intitialize 
c 

02 DO 00000 lal .4 

" hi!’;;- 

R 6000 CONTINUE 
BXb 0 . 

07 BVa 0 . 

BZ- 0 . 

BF- 0 . 

i POl «0 
Pn 2 a 0 

OM2a2*NO*0 
IPnSaB 

y.v’ 

6 MMa0 

7 NEKT«I 


ON 

vO 


ORicrv'’!- la 

OF. PCOrt QUALITY 


•PLUS V0P-51E 

/TB| *LL/WR 


tsiesmi 


15-OCT-B? 


PAOC' 4 


NMXbI 4 
T|«|97 


TI«|979. 98770 . 

ckL c«lChp re. . 0 ., 2 * 0 ) 

CALL C4LCHP f0,.«,.0,2) 
call C4LCHP (.43, *93.0(4) 
CALL CALM** (.2, ,5. 0,3) 




**** data Plotting section 

DO 10000 IPOSal.3 
JM8D2a0 
Ioaybioay? 

IPLOT«0 


DO l^®J5|{85^^J^5jI!kND.riOAT,EO.NFR)) 8KP* ( NLP-NPPf I ) / 1 00', 

DETERMINE WHETHER TO PLOT OR NOT ACCORDINO TO THE SKIP' 
FRACTION aNd ThE NUMBER OF POINtS ALREADY PLOTTED 

fh^feeiB 


c 

1 0010 


( r j 

0020 


C 

C a** 

C 

1 0020 


IF {IP08.NE.2) CO TO, 10000 
IF (IDAT.NE.NLR) CO TO 10000 
CO TO 10050 


determine satellite POSITION 


ifkrn^isit ) JbUM,JM9D.(VrjJ),JJal,S) 
IF (JMS0-LT,JMSD2) lOAYBjOAVfl 
TMajMSO»TS^IDAYft8S400B0e. 


JHaO?BJM 80 
CALL' POSIT 


f |P08,TM,IBMJ0,R0REFT,STDTM,TS,N0M2) 


IF (ifoS^L^IS) CO TO 10030 
•*** PLOT orbit in geomagnetic COORDINATES 


XaYY/(637l,a8IN(40.aPl/lS0.))tl2.7B 




IBMJ0,R0REFT,8TDTM.TS,N042) 


IF (SLAT.LT.0.) IP02al 


o 


ORIGINAL PACC 
OF POOR QUALil'i 


19-0CT-8? 


p«oe< 9 


I 


o 

o 


FQPTRAS IW.PLU9 
‘* 9 *LT ,^TN 


PI 95 


01 9S 

0 | 9 ? 
0 
0 


2 ® 
99 

0)60 
0161 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


6 ? 

U 

69 

68 

67 

60 

69 

70 

71 


0 
0 
0 
0 
0 

0177 


J! 

74 

M 

76 


78 

79 
60 

5* 

60 

09 


0184 

0109 


0190 


91 

H 

94 


C • • * • 
C • • • • 

c 

1 0030 

c 


V0?-9|F 19100101 

/TR| *lL/0R 

Gc rn 10000 

OFTERMtnI GFOMaCNFTIC components, subtract from maqsat 

COMPONENTS, AND FTNO MAXIMA 

CALL' FOCr J,mu,NEIT, 9L«r , SLON6,0, T1 ,NMX,l« 0)',BY,B7,BP) 

[F (IPIOT'EO’I) TB*TM 
If IPOSjo.n TF-TM 
1 ; "Bx 
12 ■0'^ 

4 bBF 


1 1 000 


I 1 000 
C 

? •**' 



CONTINU^^ V 

IF (fC0,r0.2) GO TO 18091 
CONVERT TO 90V IP 8ELFCTE0 


p 000 

1 0091 


UT« 960, *JM90/B64 00000, 
SHmSINT (UT49|.ONC-90,y*P 
CH«Cd8ThJT4 9LONC*90. WP 

BS0V(9)« FLOrJ) 

B90yi4)« fCo(4) 

00 1^000 lal,4 
^LOrlliiBSOVll) 

IF MFOS.GT.p 00 TO 
MAX( nalMAli (ABBT^O 
CONTINUE 


li!! 


'-l.l 

H 

H 


Iff?; 


Fmax(D) 


IF (IPOS.EO.n 
If Moat.otInFR) 


60 


d' 


10000 
0 10040 


66 

FMAXM laAHAXl ( 

67 

FMAXTJpFMAxn 

60 

FMAXMIaAHAXI C 

09 

FMAXr4)aFMAXf S 


MAX(3),FMAXf4)) 


0195 


• • • A 

• • • A 
A A A A 


90|0 

3020 


NORMA' IZE P'^OT 
THAN 90NTi TO 
PLOT LIMITS 

00 


TO NEAREST ! 
nearest 10NT 


fISfc. 


TH GREATER 

tE thIbei 


J 3000 IBI.4 

F TFMaX( 1 ),LE, 90 .) 00 TO IS 0 I 0 
HAX( T)« 50 ,*TnT fFftAX(T)/ 90 , 4 n 
60 TO {9020 
FMAxrn« 10,*1 


<L«I .? 


INTfFMAXni/lR.Al) 





• FORTRA*< IV-PLU9 V0?.8lt 
( '^SPLT.FT^J /TR|ALL'/KR 


tBl0Rt01 


I 5-OCT-B? 


PAOE> S 


96 

9F 

96 

99 


0200 

0201 

0202 

0203 

0204 

0205 

0206 
0207 
0200 

0209 

0210 


HU 


02 

02 

»2 

02 


02 6 

0219 

0220 
0221 
0222 

0223 

0224 
0228 
0226 
0227 
0220 

0229 

0230 

0231 

0232 

0233 
02*4 

0235 

0236 

0236 

0239 


YL«(B.n*2,Af.B 

'i«iA4*r tco*»i ) 

CALL Symbol (kl«yl*.2b. 


ATrx(M),0., n 

XL*. 24 

NljMBFR (XL,YL,.14,FMAxm,0.«0) 


ih. 

XL*,>4 


I 3000 


COM 


YL*5,4*(B-n -,09 

CALL Number (xl,yl, . i4,fm,0,,0) 

FM.^FMAXtl ) 

H^iJ^'NuiSEi' ( xCJyL, ,14,FM,0,,0) 


? **** DETRMINE VALUrs OF POINT'! TO BE PLOTTED AND STORE UP TO 
‘ “ PLOT ALL AT ONCE 


C 50 TO 

C 

10040 IP 

XP 


1 


I 4000 
C 

C ***• PLOT DATA 


,,*r.tl0,A(TM-TS)/f TEwTB) 

no i40BB {•1,4 

P^( J.IPL3»2,4*(B-I)41 ,2*FL0(I)/FMAX(I) 
IPL.LT.B0).AND,fIOAT,NE’,NLR)) 60 TO IB000 


0050 


DO 


1 B01 0 


1B000 


1 6000 


X«XP EM 
Y»YPTEMrn 

call CALCMP(X,Y,0,n 

w?tlt ,> 

call CALCMP(X,Y,I, n 
5o 15000 IPNpMPL 


0)) 00 TO 18010 


X«XPL 

Y»VPl 



iSVJWSS'*''''"'” 


YPtiMtf )iYPL(I»IPL) 

J ONTlMUE 


10000 CON 
C 

CALL CALCMP (X,Y.0,-B) 
^ **** BACKQROUMO SECTION 


N> 


o 9 

~n 5 

“0 9 

O 72. 
O > 
50 r- 


lO X. 

c 

1 - o 

r >•'. 





19-0CT-8? 


(J^, I V-PLU3 


o 

ro 


V0?*9ie 
/TH I ALL/HR 


19108101 


PAOe 7 


0248 
024 1 
024? 

0243 

0244 

0249 
0249 
0247 
0249 

0249 


0290 

0291 

0294 

0259 

0299 


0297 

0299 

0299 

0260 

0261 

029? 

0263 

0264 
0269 
0269 
0297 
0268 

0269 

0270 

0271 
027? 

0273 

0274 
0279 

0276 

0277 

0278 

0279 

0280 

0281 

0262 

0263 

0284 

0269 

0286 

0287 


f 3 : 3 !iuSiii;S ::3 

SYMBOL 4f4H L«T«0.»4' 

Symbol ,43 

SYMBOL jlIl0H&K5?ES!!Aj|?S*lIt 



* * * * 


DRAW LATITUDE CIRCLES 

on 


20000 

on 20 

Xa 

Ya. 


^0000 

C * • * • 

c 


-ALL 

CONTINUE 


luu I*l<a 


• CALeii 


79 

4 


(X, V, 1, 1 ) 


DO 


. AND MLT 

30000 
SHaSIN 
CHaCOS 
Xa.25*SH>l_ 
y«-,29*CH4B 
call cal 
X a9H«I 2 



j 


> Y|0t 1) 


30010 


30020 


30030 


30040 

30100 

30000 

C 


CH4B,4 
LCMPfx, 

r'— wiiT I & «78 

Ya.CH48,4 

call CACCMPfx,Y,l*n 
RNUMa(I*n*i. 

GO TO. (30010. 30020« 30030, 30040) I 
XL«12.7? 

YL«7,2 

GO TO, 30100 

XL=13.92 

YLaS,53 

GO tO 30100 

XLal2,63 

YLa9,47 

CO TO, 30100 

XL»ll,49 

YL-8.53 

CALL number (XL,YL,.14,RNUM,0,,«n 
CONTINUE 


ES"*- 
Ei 


IE 




t 


_ Ca 

ALL 


SYMBOL 
SYMBOL 
SYMBOL 
symbol 
SYMBOL 
F tMJO* 
DOCYMD 







0S»6.^I9» , 

--0B»g«B0| I 
1.94.B,B,,Ii 

(lOor, lYOC, IMOY, IDOM, lOOY, loom 


o o 

Ti ,:0 

T 

C' 

c 






a W 


vag-BiE 

/TR| ALL/>*R 


15108101 


I.O.lf. 4 , 5 H 

'* 2 * 2 ** 9 * 2 i 


Mjo*e, 

T.B.*6 

N&.S..: 


RYMp«lVOC*10000.^IMOY*100,iIOO'< 



m I {:; :5:J:; 

call symbol ( 3. 32, 5,0., 0,3HeN6,8.,3) 
Call aYMSOL 1 3 , 05 . 5 , 5 ,, 4,6HyYMMofl,0,, 
10nC«LAMJ0-l5ll2i 

i DOC,IYOC,IHOY,IDOM,IOOY,I 

hi fwist 3:lh rl- , r 


IDOm) 


CALL 

call 

c*ll 


bbo TsHi-S S:h»frlS 4 r^'" 

C 999 ,i 02 ?*lAl<*Tf*T, 0..21 
ENCODE (2.02010,HyT) lf6c 

f09g.,906.,,l4,TCKT,B.,2) 

format !lSi 


02000 

02010 

c 

c • * • • 

r 


FORMAT 


TCKT,B,,2) 


TTC»MARK8 


INDICATCI ORBITAL' STATUS 


read T2INFR1 JDUM. JTB , f V ( I) , 1 ■ I , 3 ) 
TBBjT8tIDAYA86400000. 

READ (21NIR) JpUM. JTE. (Vn)-I«l,3) 

lE»JTEA(LAMJc*J 5 M 66 )*S 54 i 0060 , 

PDT«(TE-TB1/10. 


40000 I«l.ll 
DO 41000 JJ«1.! 
X« I 


4 1 000 


Ya2.4*(5.JJ)tl .25 
CALL' CALCMP(X.t,0, 1 ) 

Y«P,4*f5-JJ)*l. 15 
^CAU' CALCMP(X,t, 1, 1 ) 

JfBfn*l)*PDT-TS 

IK a TMfTS-J day *55400001 


CHK a TMfTS-J day *55400000, 

F (TCHK, ST, 05400000,) iP^T-IPAYAl 

all posit MPOS, TM, IftMJD.ROREf T,SIDTM, 

MbTM4TS»IDAy*S6400000, 

ALL H0MI5E (TEXT.TM) 

ALL IyMBOL (XL,1,». 4, text, 0,, 5) 

ALL number (XL, ,5.. 4,SLAT,0.,i( 
all number (XL, .6,, 4,SLONft,S.,h 
ALL number (XL, .4,. 4,8ALT,0,,f) 




ORIG'NAL P.-GsZ U 

OF POOR QUALITT 



PASe 9 



"c n 
o ^ 
o 



• ) 


13-ocT-e? 


PAQE 1? 


C 

c 

c 


SSii 

000 J 

SU0ROUTI 

RtAL*8 t 

LOCl'-AL* 

0004 

iNTECrR* 

0005 

RH»A|NY ( 

0006 

RHBAfWT ( 

0007 

RSb41NT{ 

0006 

lHM3BjfN 

0009 

EMCOOE ( 

00)0 

01000 FORMAT ( 

00t 1 

RETURN 

001 ? 

END 


HOMISF fTEKT»TM) 


) 


'<)/lS0J»,) 


O' 


o 

cn 


ORIGINAL PAGE VS 
OF POOR QUALTY 


!V 

^IfLnC.ftN 
00R| 


-PLUS V 




-OIE 
I ALL/MR 


19119112 


1 5-OCT-02 


PAOe 1 


0002 

0003 

0004 

0003 

0006 

000; 

0006 


c» 


subroutine fog (J,MM,N'fXT,Ot*T«OUONO,0,TN,NMX,L,*,V,Z»F) 

•• •• • •••••••••• •* rplative to ELLffFSOro' 


J.EO.0 

j’.fQ.B 


!S(^5?.if'J58S6.S.?eSr' 

Q^'^uT ne^o^gg“?ONeNt ? 


IN CEOOFT 


RDINATE9 


NORTH, EAST, vertical 


J.Nf ,0 

J.NE ,0 

MM.EO.0 
mm.Ne ,0 

NE«T,EO',0 DO NOT 


00002400 

>01002900 

01002600 

00002700 

00002000 

00002900 

00003000 


LAT.tL0*^0 IN SPHERICAL) COORDINATES. OrOEOCENTRIC' RADI JSfKM) i|eB3|00 

OUTPUT Flo components north,£ast, vertical' In sRhericali coo 

USE DEFAULT VAIUES *C»6376. 1 9iELATr290,29 
INPUT VALUES FOR AE»7lAT 6n FIrM CalL TO FOO 


70.19.FLA1 
5n FI 

FOR external FIELD PARAMETERS 


NEXT 

NEXT 


NEXT.NE 

DLAT 

OLONO 

0 


nmax 

NMAXT 

nmaxtt 

NMXTTT 

x’.fO.R 

K ,NE .0 

T7ERO 

ABAR 


DO NOT read input VAL'JES F 
WHEN L IS greater than 0 

,EO,0 DO NOT Evaluate ExtERNAL pielo from model 

>E.0 REAg INPUT VALUES FOR E* TERN»L* F I ElO PaRA 

external field model 


meters 4HEN 


LONGITUDE 1 * DECREES 

gE8?i;iS,^!Uy85 iS5! 3SK i:f 

MAxIMUM^OFCREJ *^?pJ5®55o®5 TESmR OP'^FIILO^MODEL 

" ■ " sIcONO ■ ■ ■ 

■ ■ • THIRD ■ ■ ■ ■ • 

PIIt8 2SSFt ?SFFFIFIF3l? 

EPOCH TIME FOR FIELD MODELI COEFFICIENTS 

MEAN RADIUS UJEP IN FIELD MODEL' POTENT I AL' EXP ANSI ON 
(DEFAULT ■ 6371,2) 


MODEL 


modext.fo.0 no external field solved MITH MOI 
MODEXt.NElo external MFLo SOLVED 0lTM MODEL 
L.eo.0 evaluate field 


R00003200 

00003300 

00003400 

00003900 

00003600 

00003700 

00003800 

00003900 

00004000 

00004100 

00004200 

00004300 

00004310 

00004320 

00004330 

00004339 

0O00433R 

00004340 

00004349 

00004390 

00004400 

00004900 

00004600 

00004700 

00004800 

00004900 

00009000 

00009100 

08009200 

00009300 

00009400 

00009500 

00009600 

81009709 

00009000 

00009908 

08006806 


000061 

000062 


EOIIIVAIENPE (SHMiTfl . 1)»T 
common /COEFF8/T6n6,l6) 

COMI'^nN /FLOCOM/sT,CT, SPH, 
ABAR,El£F"lt3l^EXtF 

8i:!Sil§3 


T6(i,in 

CPH,R,NMAV,BT,BP,9R,B, 

HMir(ie,|81,AID(30) 


16066300 


DATA IFRST/0/ 
DATA AF, 


FLAT/637B, 16,298.29/ 


08006408 

00006500 

00006600 

00006700 

00006800 

00006900 

00007000 

00007100 


rS 


FORTRAN IV- 
FTfLOG.r TN 


PLUS V0?w5lE 

/TR I ALLWR 


151 191 12 


18-0rT»B2 


PABC 2 


90R9 


00 

00 

00 

00 

00 

00 

00 

00 

00 


9 
1 
2 

4 

5 

6 
; 

. B 

0019 

0020 

002 1 

0024 

0025 

0026 
0022 
002B 

0029 

0030 
00: 
00: 
00: 

0034 
^ 0035 
o 0036 
0032 
0036 

0039 

0040 
004 1 
004? 

0043 

0044 

0045 

0046 
0042 
004b 

0049 

0050 

0061 

0952 

0051 
0064 

0055 

0056 

0052 
0056 

0059 

0060 
0061 

0062 
0063 


1 03 


2 

1 06 


1 02 
1 02 


10 

12 


DATA 


AST/0./ 


i 9t\ 


-.21.2/ 
1 10» i A 


ATA {linA 
F(1FR3T3 
on! INUE 

FRit*! 

LAT0I. -l./FLAT 

rg- 

- ?»0, 

e3«0, 

A?* Af **2 

A4« Af A *4 

82 »(*E*FLAT )*«2 


0.110 



NMXTTT,^00E1rT,K,TZER0, ABAR, 


UXNbB 
TEMP.g 
READ (1,6 
FORMAT 

IP V I uu I vtr 

MAXNalMllBrN,MAXN)3 
C(N,M)mCNM 
GT (N.MIeOTNM 
0tt(N,M)B6TTNM 

r 

Of (M-.' ,k)m . . 
Ott(M..l,Al)«HTTNM 


M,GNM,HNM, 0 TNM,HTNM, 8 TTN'<,HTTNM 


(N,M)mCNM 
HN.MIeOTNM 

1 1 (N|MJb4tTi, . 

fMpsAMAXl ( iFMP,ABB(OTNH)) 


- .N) I . 

,A|)eHTNM 


TT.fa’. 
6 


Ott (M.. 1 

00 To A 
If (N>.xT 

Ayrr.TtTNM 
.1) GO. TO 


0i 
6 



- IL 

JNUC 


IeHTTTNM 


1S> 

106 


^eacm»i 023 Ei*e2 eg 
00 12 Ne2ImaSa’*' 

,N 


DO }2 Ne2, 

00 12 M«i, 

HIeH.I 

8nt inA^ 

ONT iNUE 


1) GOTnl0 


00002200 

000022U 

01002300 

00002400 

00002B00 

00002200 

00002800 

00002900 

00006000 

00900100 

00006200 

00006300 

03008400 

00006900 

00000600 

00006200 

0O0OBO00 

00006900 

00009000 

00009100 

00009206 

00009300 

00009400 


600 

000 

000 

000 

000 

000 

060 

006 

000 

006 

000 

030 


0306 

0406 

6000 


000 
000 
030 
000 
006 
000 
000 
030 
(100 
000 
000 
000 
000 
000 
000 
066 
000 
000 
600 
000 
000|2300 


2401 

liii 

2200 

2000 

3200 

3300 

3400 

3900 

4100 

4203 

4600 


o c 

m 50 

tjC 
o z 

O T> 
50 r- 

o y 

5;* O 
r -1 


-< (/) 


00 





1BII9I12 


t5-0CT«8? 


^AOE' 3 


o 

03 


■*364 

4b59 

0OC6 

0087 

0068 

0069 

0070 
007] 

0072 

0073 

0074 
0076 

0076 

0077 
007B 

0079 

0080 
0(:8l 
0082 

0083 

0084 
008B 
0086 
0087 
008r 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 


1 3 
1 4 


I B 


1 ? 


210 

220 

230 

270 

240 

250 

260 

18 

19 




1 


MUf.lV. 

0 IB N«2,M 




M*lfi 


IB N*2, 

IHMIT(N, 1 J-SHMlTfN. 
SHHlTh ,N)-0. 

J J"2 


1 * 1) *7101 T(2«*<«3) /FLOAT fN-n 


0 15 M«2,N 

HmI? (N, *^J«9HMTT(N,M*n*9URT(FL0fcTf (M.Mfl )»JJ)/FL0AT(N4 6.2) ) 
HMIT (M-i ,N)BSHMlTfN,M^ 

U’\ 


6 N«2,MA)(N 
6 M«i,N 


M)49HMlTfN,M) 


?l!Nftx?T?IiTj0!ANn!ll,LE.8lGTTTfN,*n»OTTT(/J,M)65HMIT(*l,4) 

If fM,Eo,n ftotoiB 

fMNMX^T?rG}(0n^o!N?t!F.8loTlf ?i.l,N)«CTTTfM-l,»0»8H'<lT(H-l,M) 
CONTINUE 
TbTM.TZERO 
00 16 NalfMAXN 
“ 18 M«1,N 


00 18 Mb 
CXM0. 

HX«0, 

M.fo.l) go TO 
■ M6XTTT1 G 
iMJat 


FfM 

FfN 

GX 




HXbGTTTN-I.NIaT 
FTN^GT.NHAXf n 60 1° 22® 
GXarTG* ♦ Gt fN,M)J*T 


(N 
QXb 
HXBitHX 


BT 

8 


^Xb 

GXbTGXbGTn^M) 
HXBTHXfO|M-i ,Nj 
TC(N,M)BTdx 
TGTM-1 ,N)«THX 
Go TO 10 
ONTINUE 

F(N.GT,NMXTTT) go to 240 
GXBfiTtf (N,M)*T 
F (N,GT .NNlAxtTl GO TO 290 
CXBf fGx*GTT(N,M) )*T 
[F (N,6T ,N6AXT) GO TO 260 
'CXBflGi4GT(N,M))*T 
GXb tCX4Q(N,M) 
TGTN.MIbTGX 

OLATRbDLAT/B7,299779BO0 
SYNLAbSINTOLaIP) 
RinNGB0L0N0/B7. 2957 79600 


000|S|00 


000 

000 

000 

000 

000 

000 

000 

000 

S 00 
00 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
020 
000 
000 
000 
000 
M00 
B?0 


00 
300 
400 
5S00 
3600 
B/00 
5800 
S900 
6000 
6100 
6200 
6300 
6400 
6B00 
6600 

tlSS 

6900 

7000 

a®® 

7200 

7300 

7400 

7900 

7600 

7700 

7800 

7900 

8000 

8100 

8200 

8300 

8400 

8B00 

8600 

8700 

8800 

8900 

9000 

9100 

9200 

9300 

9400 

9900 

9600 

9700 

9800 

9900 


00020000 

00020100 

00820200 

00020300' 

00020400 

00020B00 

00020600 


>o 


;^;cod 'JQ 



rp5TR»H IV-PLU3 VBg-Bir 
►ifLnG.rTN /tRi*lL/rr 


0 ?0 
8 21 
0 ?? 


23 

24 

23 

28 


0127 


28 

29 

3a 

w 

34 

35 
38 

37 

38 

39 
4B 

41 

42 


0143 
0{44 
01 45 
01 48 

0147 

0148 


20 

C 

C 

c 

r 


21 


22 


C 

23 


15119112 


1 5-0CT-82 


:PHrr g8(RL0N8) 

IPH-ilN/prONGj 
:F (J.Efl.B) o6to20 

0 IS 3EOrENTRIC RADIUS RHEN Ju\ 


RbO 
CT«81NL4 

50 t6 21 

5|NL*2 


a. 


SLA**2 


0 IS QEOOETIC altitude RHEN Jb0 
ALTm9 




-EN2R?2ll2B}»i*^ 
OEN«SQRT(DEN2) 


INLA2 


F*C«(((0*OENWA?)/nO*OEN2}B2n**2 
CTbSInlA/SORT (FAC*C dSLA2t4lNL42) 
R>SOftT]Q*(Ot2,*DENW(*4-A4B428INLlAI 


ST«80lSt ( 

NMAVaHlN 

NEXTFbNF 

F»B 

Z«-BR 
RETURN 
TRAN8F 


nb|nmx,mavni 


L«2)/0EN2) 


22,23.22 


ANSFORHS 


field 

Sf *SQR 


SINp>8iNLA*8 

cosdbsorTm 
Xa.BTaCOSO*^ 

ZaBT*SlNO*BR*Cn 

return 

END 


JO^pgpOET Ip .DIRECTIONS 


. .-08L42)*CT 
•SIND**2) 

*StND 
D 


PAOE 4 


(((1020700 
0b020000 
00020900 
00020903 
00020905 
00020900 
00021000 
00021100 
00021280 
03021300 
00021305 
00021 
000f ' 

0001 
03021 401 
00021500 
00021600 
00021700 
00021 000 
00021900 
00022000 


>2100 

'2200 


0001 
000L ^ 
00022300 
00022400 
00022500 
080|?6O0 
00022700 
00022800 
00022900 
00023000 
00023100 
00023200 
00023300 
00023400 
00023500 
00023600 



J 1 V-PLU 3 V 02 - 91 E 

r FIPLnC.FTN /TB|AlL/**« 


19119154 


1 9 - 0 CT- 8 ? 


P40e- 1 


8001 

808^ 

8083 

0084 

0005 

0005 

0007 

0005 

0009 

0010 

00 I 

001 ? 

0013 

0014 

0015 
0018 

00 7 

001 B 
0019 
00^0 
00|l 
008? 
0023 
00?4 
00?5 

_ 00?6 
00?7 
o 0088 
0089 

0030 

0031 
00 3 ? 


0033 

0034 

0035 

0036 

0037 
003 B 
00 39 
0040 
004 1 
004? 

0043 

0044 

0045 
0045 


9 

r 


9URR0UTINF M40F 
tOMMnN /COEFFS/CMBtlS) 

COMHON/Fi 0C0M/9T,CT,SPH.CPH,».N'44K,9T,9P,BR.B, 4B4R.E1 i E2.E3, *<{*1 
DIMENSION 8M8,iS),6P(l&,l8),C0'49Tn8»lS),98fl8),C^(ll)»8Nf|9),F 

IF (P(l,n,EQ.l>) 00 TO 3 


FN(N)i 
DO 2 w«i,s 

C0NST(N,»4).FLO4Tf(»<-2)**2-(‘<-n**2)/FwOATf(2*N-33»(2«N-9n 

SP(2)«9RH 

CPf 8 )«CPH 

DO 4 Hij.NMAX 

5PtMl«9Pr2)*CPfM.l jACPf2)*9PfH.l 1 
CP]MS«rp^2)*CP?M.n.SP^2i*8P(M.I) 

IoAbabIr/Ii 

ARbAQR**? 

9 T* 0 , 

0Pb0, 

RR>0, 

DO 8 N*?.NMAF 

ARbAOR* A A 

DO 8 M» 1 ,N 

IF (N-mJ 8 . 9,6 

AfN,NlBST*PTR-l ,N-n 

DPfN,ft)B 9 T* 0 P(N-l ,N -1 )*CT*P(S -1 ,R-n 

CO t 6 7 

PfN.MjBCTBPl^J-l.Mj.CONSWR.m^PCS.J.M) 

NOTE I CON9Tf2.11*0 

IPfN,M)BCT*OP(N«!,'4)«BT*P(N*!,M).CON9T(N,M)«OPfN-2,'4) 


PaPbP (N, H) *AR 
IF (N,E6,h 00 TO 9 
fEMABft(N.M«CP(M) 4 G{M-l»N)B 8 PfM) 
BP»BP-?Of N,M)« 9 Pf » 4 )«C(M«i ,N)*CP( 
CO TO 10 

TEMP«C(N, ' 4 )*tPfM) 

0TbHT4IENP*OP(N,M) 

br»br«tenp*fn<n 5 *P 

8 PB 0 P/ 9 T 

IF (NEXT -CT, 0 ) CALL |*TFLO 
BBdURT(fiT*ATABP*R8«BR*BR) 

RETURN 

end 


PC'AITaFMfHlBPAR 


• AR 
•AR 


00023780 
00023800 
00023900 
M( 0 B 024000 
01024100 
00024200 
00024300 
00024400 
00024900 
00024600 
00024700 
00024800 
08024000 
00025000 
00029100 
00039200 
00029300 
00029400 
00029900 
00029600 
00025700 
00025800 
00025900 
00026000 
00826100 
00026200 
00026300 
00026400 
00026900 
00028600 
00026700 
00026600 
00026900 
000289 j 0 
00028919 
08026920 
00027000 
00027100 
00027200 
00027300 
00027400 
00027 B 00 
00027800 
00027800 
00027900 
00028000 
00028400 
00028900 
00028600 
00028700 


-< t. 


oo 


Sy«HOuTi«<r EBTrLf) 

T}.h*si 

4 :U:w 


aasi 

a?8? 

2283 
8884 
2285 
8385 
8288 
8385 
8888 
881 8 


't' 

-Tl*rT 

♦T 1 *5T 


I5I28I87 1S«0CT»B? PASS IP 

CPM,R,>(^ltf,BT,8P,SP»B.«B«R»ll,E?.r3 


P*BP45i*5PH.E3*CPH 
BTrBT»T 2 
Rf Turn 
E»iO 


B882SS0B 

88P|09«I 

8IP29BI3 

080291BB 

B8P29283 

880293P8 

P8P29408 

01029900 

01029608 

00029700 


^5 


C'K.GINAL page 
OF POOR QUALI 





?85!Pt:i 


0003 

0004 

0005 
007S 

0007 

0008 

0009 

0010 
0011 

0017 


IV-PLU9 V 


15120113 


lS-0Cr»B2 


P46E 1 


SUPROUTINE rODriMSODfRHMS) 

CONVERTS MlLISeCONDS OF THE D«Y TO HOURS, MINUTES, SECONDS, AND 

thousanoths Of seconds, 

INTERMEDIATE yARllBLFS ARE USED TO CONSERVE MULT IPUBCAT I ON 
THEREBY INCRFASINO SPEED, 

INTEGER* 4 IMSOD, IMSOH, IMSOM, IMSOS 

HOO-IMSOO/3600000 
MSOHa HOD*3000000 
MOHa( M90D«IMS0H) /60000 
MSOMs MOH*60000 
SOMa( M90D»IM90H«IMS0M) /1000 
MSOSb SOM*1000 

. TOS« I MSOO-ImSOH- IMSOM- IMSOS 
RHMaBlHOO*l 0000.4lMUH*100,*ISOMtITOS/ie00. 

RETURN 

END 


O 
- 1 

-r 

C 

o 


<. 


00 



t9l2Cll20 


lB-OCT-8? 


PAOE 1 


fORTRAM IV-PLU9 V02-91E 
POSIT, PTN /TRiAlC 


/WR 


0001 

000 ? 

0003 

0004 
0009 
0009 
000? 
0009 

0009 

0010 
001 1 

0012 
0013 
00 
0019 
00 ‘ 
00 
00 

0019 

0020 
0021 
002? 

0023 


0028 

0027 

0028 

0029 

0030 

0031 

0c:? 

0033 

0034 
0039 

0038 
0037 
0030 

0039 

0040 
004 1 

0042 

0043 


PDSfT . 0FT9 tOOROINATES, CALLS 8T1ROB AND SATPOS TO INTERPOLATE 

SUBROUTINE P03I T ( I POS . TM, I BM JO, ROREFT . 9IOT M, T9, N0H2) 

COMMON /ORBCOM, XT EM ( 6) , HE M ( 9 ) . Z T EM ( ST 

COMMON /ORBIT/ XX,YY.Z?,9Lit,9L6*40»8ALt,lOAy,RHEM 

dimension COf31,P6S(3T 

REaLoS T#TM,tS 

lNTEr.|R*4 ISmJo, JMSO, joum 

DATA PI / 3,J419R299 / 

DATA ONE / 1.0 / 


C 

10 


K M INT(TM/60000 
IP (K .LT. 31 K 
IF fK ,QT, n6i 


SM2) 


♦ 1 ) 

I 3 

K ■ N0M2 


JJ 


■ 0 

JJ ■ 
KK ^ 
READ 


JJ ♦ 1 

'MSiS' 


C 

c 

c 

1 1 


I 2 


RAD ■ 90RfTcflcn«*2 ♦ CO(2)«»2 ♦ C0(31*(*2) 
KK B KK ♦ i 

R^AD P^KK) joum, JM9^jFLAT,MR,0LAT 

GO TO 100 
IPOS B 4 

CONVERT TO MAGNETIC INVARIANT COORDINATES 


RMEM) ) ) 



JM9D 


T9 



S; 


o 9 

m ?o 


•33 0 

o 

o 2 

» r* 


o y 

r rn 


“t — r 

■< (A 



r3»T<»*'< iv» 
5 1 T . r tm 

»L09 

/ ' ® 1 *LL/"* 

1 9 1 ?B 1 ?8 

1 9-OCT.S? 

P49C 8 

8344 188 

8049 181 

88 4 9 

IPOS • 0 
ft f f 0»*« 

E sn 





c n 

•» , j 


X 

■y» 


l uOli QUALIFY 



r'- 


V-PLUS Vg2» 
N 




15120149 


1 5-0CT-8? 


P40C i 


aaai 

0082 

aaai 

0004 

aaas 

aaas 

29ir 

aaas 

2809 

aai 0 


aa 

00 

ae 

aa 

aa 


0015 

aai 7 

00 5 
0019 
00?0 
aa?i 
aap 
00M 

0024 


0025 

0025 

002 ? 

0025 

0029 

aai0 


•*** STIPOB - iNTEPPOLATfS POSITION 5IVEN PIVC COOPOINATES 


20 

C 


Subroutine stipob 

COMMON /O00COM/ XT 


Eommon /ORslf/' xJi vylill^Lif I J l 651 jsI^ 

6l MENS ION OFx(5,5J,0^t(S,5),6P?(5,%) 


I0*Y 


r LC ■ 0 
■ DT/A0a00. 

U2 • U/2. 

U3 > U«u5 

U4 ■ (U2/1. )*(U3-.B) 
US ■ U2*U4 

Erv! I 

CONTINUE 


00 28 L ■ 2, 5 
N2 ■ 5 - L 
DO 20 N ■ < 

0F*fN,L 

CONTINUE 


) ,0T 

■ 

• YtEMi 

• ETEMf 



9999») 00 TO 50 


, N? 

■ DPX(N4 

■ OFY(Nf 

■ DPZTN* 


l!0i : iijliltlil 


V7 • DFx f 3* n ♦ U2*(0F« (2 
» U4*f0FK(I »4) ♦ 




♦ u 


univ 




.1) 


yy ■ DFyfS.n ♦ U 2 *(OPY( 2 i 2 )* 0 Fyf S,2)) ♦ Ul«0Fy(2,S) 
♦ u4»roFyfi»4)40f y(2*4)) * uS*Dryn,6) 


ZZ ■ DFZfSin ♦ U2*(0FZ(2,2)t0FZ(3,2) 1 ♦ U340FZ(2,3) 
1 ♦ U4*fnPZf 1.4546 f?( 2,4)) « u5*0FZfI,B) 


RETURN 

jfifS.; ' 

f NO 


a* 

a 


owiGiNAU r..Gi- y 
OF POOR QUALIFY 





PLUS V02-61E 

/TR| *LL/»<R 


I si2sn 2 


15-OCT»02 


PAOE 1 


0001 

0002 

0003 

0004 

0005 
0005 
0002 
0005 

0009 

0010 
001 1 
0012 
0013 
00t 4 
0015 
00 5 
00)7 
0015 

^0019 

0020 

SSU 

ISI 4 ’ 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0034 

0035 

0036 


SATPOS . FINOS LATITUDE, lONOTTUDF, AND RADIUS PRON X,y,Z 


10 

C 


SUBROUTINE SATPqS f 1 POS , POS . DEL T DT , UP R AD V , S 1 01“) 
COMMON /ORBIT/ TV, SLAT, SC0N6 , salt, I6ay 

DIMENSION P0S(3) 

DATA R360, RIS 0 / 360., 180, / 


ROT ■ DFRAQY06;3003e94 4 DEL TO Y*. 0 1 7202 1 4 ♦ 8IDTM 
F (I^OS ,CE, 5) ROT ■ 0, 

IT ■ SORTJpoSI) T**2 ♦ P0S(2)a»2) 

IF (TT ,nE. 0,T 00 TO 1 
RLONG 
SPN ■ 

CPN ■ 

Co TO 
CSA ■ 

SNA ■ 

SNL ■ 

C3L ■ 

SPN ■ 

CPN ■ 


s ®* 

I- 

lO 
PO 




COSJROI 
SNAaCSL 
sna*snl a 


CSAaSNL 
CSA4CSL 

RLONG ■ 57,295779B*ATAN2f SPN,CPN3 

F (IPOS ,LT' 3) CO TO 2 
F IAlONC .Lt. 01 RLONG ■ 360. ♦ RLONQ 


I 

IF f ABS(RLONG) .LE, 100.) 00 TO B 
RLONG ■ RlONG • sIgnIRSBB, RLONG) 
CO TO 2 


0 TO 5 


^LONG ■ RLONG 
R ■ S( 

R2 ■ Pnsj3)/R 

RCaT • AilN(RZ1*57.295779B 


■ SORTfPOSn )**2 ♦ P0S(2)**e ♦ P0Sr3)**2) 


GO 

SLA 

SAl 


»LA 
GO TO 7 


ABSfRlATI ’lF. :3,1 60 TO 10 
■ SlONiRl A0,flLAT ) - RLH 


RLAT 

R - 8371 .0 


RETURN 

END 


00 


• l) JO 



1 Bt23|?l 


I5«0CT«B? 


f''' wr.n. 


IV-PLU3 V02-BIE 

/tB| ALL/HR 


PABE i 


*•** time - CONVERTS HOUR-MINUTE-BECONO OEl DAY TO MILLISECOND OFi DAY 


0001 

0002 

subroutine 

realms RMS 

0001 

RH ■ AInT 

0004 

RH ■ AINT 

00ES 

RS ■ AlNT 

000B 

RMSOD ■ 36 

000F 

return 

000B 

END 


time fRMSOD,RHMS,RTOS) 


RHXl0000.)/ieB 
RH*i000e. • A 
B00eA,*RM 


♦ 10I 


>(RS4RT08) 


00 

00 


o o 


TJ ^ 

o . 

O 

» r* 
O '0 


5^1 

r~ ' T 


-< O) 



r^, '‘M-.Vd 




0002 

0003 

0004 
0009 
0006 
000r 
0006 

0009 

0010 

00 t 
00 ? 

001 3 
0014 
0019 
0016 

001 7 


/NR 


1BI23I28 


lB-OCT-8? 


1900 JAN 01 ■ OOC 1900000 
tNTEGER*4 lOOC 
COMMON I0 PM(i5i 

lOOY«IOOMtIOPMf JMOY) 
YOC«lYOC 
AUEAPbVOC/A.B 


. S181I 

S0O 


01900 

01 SI0 
01920 


IdoCmI YOC* 36941 LP A P4 100 Y -00000 
GO To 01920 

lDOC*lYOe*3694lLPAP4lpOY-00001 
iDONBiooc-Tnooc^i )/M»>40 
ftETURN 


lETURN 
ENO 


PAse> 1 


CD 


0 

1 


o 


c 

c . 


( • 


/r 


I I 


* 1 . • 

-< c. 


00 

vO 


k tki ^.w 


FORTRAN 

OOCymO. 


300? 

0001 

0004 
0309 
0006 
0007 
0006 
0009 
00 0 
00 t 
00 ? 
00 3 
00 4 
00 9 
00 9 
00 7 
0016 
0019 
00?0 
00?1 
00 ?? 
00?3 
00?4 
00?9 
00?6 
00?7 
00?6 


FlN 


PLUS V 


10123136 


|B«0Cr«8? 


V02-O1F 
/TP 1 4LL/XR 

9UBF0UTISK OOCYMOfTDOC.IYOC, I00'*7*P9V»II)0'<I 

CONVERT 0»Y OF f»NtuRY fOOC) T 0 v|*R.MONf H«0» Y (YMD) 

1NTEGER«4 lOOC 
common IORM(i?) 

IOOK-IDor.( (IDOC^I )/7)*742 
5oc*rIooC4000e] )*i0e 
I YOCpOOC/30S?9.0 
YOC"! YOC 
ALE*PpYOC/4,0 
"LEAPbYOC/4,0 


PAOE- I 


32000 “ 


nYeIOOCi0ft001.lYOC*3a8-lLt*F 
F f ALEAP. 6 T, ILEaPjCO to 32000 
F (ALE*P,E 6 , ( 0 , 0 ) UO to 32000 
^nrB^Do940At 


32900 MPy»MPY-01 

AP.fiT.I 


f 


0 TO 32600 

si'fs j ? Jr ‘ 

M(MPY) 4001 ) )60 TO 33200 


ap:eq: ( TQ.32S00 


Y 

..OYICT. ( I 
0 3$90A 

F (igOY^OT^lOPM(MPY) )00 TO 33000 

33000 IDOM.iOOY-IOPMfMPY) 

CO TO 33900 

IF (If'OY.LT, 061)00 TO 33000 
lD0M»10flV-( T oAM(MPy) 400 1 ) 

I MOYbMPy 
return 


32600 


33?00 

33900 


O c 

•n , • 

-o2 
O 2 
O S 
73 r“ 

O ^ 

c 


FNO 


00?9 


r 



tSl28|2S 




rO^TRAN 

CUR0T3.F' 


Vr?LUS V| 
N /’ 


-sii; 

1 ALC/RR 


I5-0CT-8? 


PAOE 1 


* • * • 
• • « • 


fuTS* r)A?K'*FILE *NO 


IN> 

O 


8081 

0082 

0083 

0004 

0009 

0008 

0007 

0008 

0009 

0010 
001 1 


0012 
031 3 
001 4 

0019 

0018 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 
0029 

0028 

0027 

0020 


0029 

0030 

0031 

0032 

0033 


ON 

§5 J 2 fl 

L4NB04 

T4 PI / 3,14189280 / 
T4 RF / 8571000,. / 

TA aLTI / l5000A. / 




TP(4,8,72) 


01000 

C 

01010 

C 

01020 

C 

01040 

01090 

C 

0 1 080 


NRITF (9|01000) 

FORMAT (' FfiTER NUMBER OF CELL RINGS '3 
READ (9i«) NlIMT 

format 'en in NUMBER or CELLS ^fR 380 DEO LONOITUDE' <) 
READ (9,*) NUMl 

RRITE (9*010201 

FORMATS ENTER INNER AND OUTER COLATTTUOES OP RIN0 l) 
READ (9,*) CLl. CLP 


R MFX RADIUS OF CURRENT FILAMENTS AND •) 


FoJmATM *LAT? fuDlNAL THICKENING EXPONENT •) 
READ (9,*I RF, DF 

RRITE (9,01060) 

FORMAT ( T ENTER MODEL NUMBER » 1 
READ (B,*) NCODE 


• • * • 
* * • * 
* • * • 
• • * * 


LOOP 10000 DEFINFS THICKNESS OF CORRECT FILAMENTSl 
FOR A-A CURRENTS SUPPLYING E-R CURR^tS ON I ■ 11 
FOR N.S CURRENTS ON I ■ 2f 
FOH F-R CURRENTS ON I ■ 3, 


DO |^0000_I_* .2, 4 


10000 J ■ 1, NUMT 

no 


1 


10000 CONTINUE 


10000 K B 1, NUMl 


(RFb*2>) 


Cr I'UOR QUALITY 



FORTRAN TV-PLUS V0?-9lE 
f '' rjnoiS.Ft!^ /TP|*LL/XR 


1 5I2SI79 


I 5-0CT-87 


PAOE 2 


0834 

0030 

0036 

003F 

0030 


0039 

0640 

0041 

0042 

0043 

0044 
0049 

0046 

0047 
004b 
0049 
009t» 


*•** LOOP 20000 DEFINfO THICKNESS OF CURRENT. F !L 4ME NT S 
•••* FOR F-A CURRENTS SUPPLYING N-S CURRENTS, 

N ■ NUMT » I 
no 2g000_J_a I, N 
0 


I 


20000 CONTINUE 
C 


90800 K a I , NUHi 

TP(|,J,K1 a j JJjJJj|J5^**CLl)/NUMT ♦ CLD 


') 


C •••« LOOPS 30000 and 40000 DEFINE CURRENT PER LOOP FOR N-S AND E-4 RPSP. 

C 

DO 30000 J a 1«NUMT 

00 30000 K R. 1 * NUHL 


if" tKlft(!ftUML/2, ) S a 
AHPi(f.J.K) a S 


10000 CONTINUE 

00 40000 J a 1,NUM7 


no 40000 K a 1, NUML 
i a 0 . 

'.CT. 12.#NUHL/24.) 3 a 0. 
lPS(2,J,il) a S 
40000 continue 


0091 

0032 

0039 

0034 

0053 

0056 

0057 
005B 
0099 
0060 
0061 
0062 

0063 

0064 

0069 
00 66 
0067 

0066 

0069 

0070 


* • a • 

• • * * 


llNn^?M?fR"L0CA?l8N*lS SrACE^ LENGTH OF CURRENT FILAMENTS 
2 

. N 

pJ-n/d^NUMTl ♦ CL'l j«FI/lS0. 


I 


00 B0000 I 

00 B0O00 J a I, N 


TMU 

PZF 

R§( 

RT? 

B0000 CONTINUE 


a a » r» 

: 

a ATA 


♦ .D/TMUd, J)1 


RId,K) ■ PM 

RI(2,K) 


NUMT 

t ALjn JJC09(CLd.K)-CL(2»K)))/CCLf?»K) 




^tl>i5i|NfCLl(f*lf)- CLd»K)) 

(K)Ua 


RE 

RZE(K) 


IN(PI/NUML) 


'O 
I J 


'll ^ 
o V?: 
o > 

33 f 


O “O 
C j’ 


o 

<•3 




FQ!TR 45 < IV-PLI'S V 02 * 5 jir 


f J 

r j 


RR 7 I 


aa^? 
ea/) 
aar 4 
aefs 
aare 
aa/r 


asTB 
837 3 
8398 

aasi 

aas? 

8893 
8094 
8899 

8890 
8897 
8899 
8899 

8891 
aa^t 
889? 
f '^291 

8894 
8899 
889e 
8897 
8899 


LL/"R 


1 9l28l?S 


1 5-OC r-92 


8A8C' 3 


68888 


corttnije 


».*• Lnof 70888 OtPTNES THE »INC CUR9E»*T 
8INGA«4.*RE 

8INC8B4. •8E*TAS(8I/»<U‘«L) 

00 ? 0 f »0 

TPRM)»l.M9/f 1808*»2) 


78800 

C 

C ••• A 

c 


AMpdjf 1 >■! . 

CONTINUE 


OPEN (UN?Ti 


HR 

HR 

hR 

HR 

HR 

HR 

mR 

H? 

HP 

hR 

HR 

,.R 

HR 

HR 

HR 

HR 

HR 

HR 

HR 


CLfiS 


^NAME««g^9,OAT» ,TyP|»INEH* ) 
CLl^ 

tiU IhI»huhM 

i: 

Ja 

• J" 

a a n a B 




8099 
81 88 


STOP 

ERR 



£ 




fqhtrah* iv-PLua 
AHPL'T.riW 




/MR 


1 M2nifl9 


1 5*0JT*S? 


p*(jet 1 


0001 

0002 

0003 

0004 

0005 

0006 
0007 


0008 
0009 
00j0 


ro 

CO 


00 

00 

00 


0012 
I 3 


0014 

0015 
0018 

0017 

SS I 

0020 

0021 

0022 

0023 

0024 

0025 


0026 

0027 

0026 

0029 

0030 
003} 

0032 

0033 

0034 


AnPLT 3H0W8 CURRENT 
IT A0OVF THE I0N08 


PLQWINO THROUGH 
PHERIC CURRfN 


THE 3MRPACE OF * SPHERE' 


• ^4CH C?8cLE represents * f|tlO CURRENT flLAHENT AND 

• ABOUT 90 PERCENT OF THE PL A T ICURT 1 C * CL Y OISTRTBUtEO 2URRENT 

• IS THPrEiN ENCLOSED, 

• Each line in one of these circles represents one ahp rNUM.^ER 

• of liSes ♦ or - ,5 li currIntT, 


VO 

o- 


dimension AMPr?,9 72). TPf4, 

^pGICAL«1 PARflAY(in»‘NUM(2j 


EAL*4 INCL, Ml T 
DATA D4RRAV / »0 


0 » , 


1 0 

DATA PI / 3.141O920B 


• n, 

/ 


• S' 


8,72) 


I . I 


•O', »Al, IT', 'I 


I t 


DATA CF / 1,2394F-8 / . 

DATA INCL', THTA / 'iNCL*, 'THTA* 


WRITE (6|01000) 

01000 FORMAT (» ENTER FILE 
REAP (6,01010) FNUM 


l|‘"’ 

DARRAY($)bF 


numberinnI ) 


01010 forma 

DARRAyIiB) 


NUM(l) 

•fnumJ?) 


WRITE (6,01020) 

P1020 format (< Enter incl and theta ofi orbiti 0,0 for no orbit •) 

READ (6,*) INCL, THTA 

8INCL ■ SIN(PI*INCL/180.) 

CINCL ■ CQ8iPl*|NCL/l80,i 
ST ■ siN(P|*THTA/|l0.) 

CT ■ CO8(PfoTHTA/lB0, J 


YORG ■ 5*5 

CALL CALMPfXORO, YORO,0,3) 

**** LOOP 10000 PLOTS THE ORBIT 

.AND, (THTA 

ii 6 n 


IF ((INCL* ,EO, 0,2 .A 
00 10000 I B I, 160 

0 MP . (I-l.lpP./H 


,EQ'. 0,)) 00 TO 10010 


SMP 
CMP 
XI 


N(30.*RMFBPI/j 
.J3(30.bRMPbPI / 
STbSMP ; CT*CMP*i ^ 
ST*CMP*i 


■ ■CT*SHP • 

■ CMP*CINCL 
AT B lfl0.*(ACO8(8ORT(XLBB24YLBB2)))/PI 


OF POOR QUALITY 


I 




1V-PLU3 V0P-B1P 
• MPLT.fTN /TR|*LL/»<P 


t5l?B|59 


1S«0CT»82 


PAOE 7 


0035 

003S 

003r 

0038 

0039 

0040 


0041 

004? 

0041 

0044 

0045 

0046 

0047 
0046 


0049 

0050 
005} 
005? 


0053 

0054 

0055 

0056 

0057 
0056 

0059 

0060 

006} 

006? 


0063 

0064 

0065 

0066 

0067 

0068 
0069 


{all CACcMAfllrM^i) 
CQf^IINUE 


L/SORTf 

XL**24YL 

L/SORt ( 

XL«*2tY[ 

CALCMP{ 

X, T,0, 1 J 




10000 ^ ^ 
10010 CONttNU^ 
C 


OPEN 

READ 

READ 

READ 


;UNlT»liNAME«DARRAy,TYPE*»OLO» } 
NCODE 


tUNT 

!!;:! 6 W»;'fc 5 -L 

N ■ NUMt t I 

5lig !i;: i!!iEii)^3;j),":ii:wj/3iirw;:'»,,„ 

CLOSE fUNtT«l) 


N ■ NUMT ♦ I 
DO 20000 ST ■ 
DO 20000 I 

no 20000 


1. 2 
■ l» 

J ■ 


N 
I * 


NUML 

C 1 1 1 II 1 1 1 1 1 1 1 1 1 1 II i II I II i 1 1 1 • I M M I m i I 'll 1 1 'll 1 1 i 11 1 ii ii 1 1 1 It I i 1 1 1 1 It ) 
*••• THIS SECTION calculates THE CURRENT PER FILAMENT 

and, . f I _ .EQ._N) I 60 TO ?0000 


XL ■ -4.8Sr*c63(2,* JJ*Nt/2 ,-| ,)aPT/^UH 
YC ■ -4,5*F'*SIN(5,*(j4NT/2,-l .)*PI/NUM 
Mmp . jiMP(I, 

IE((NT ,E0, i .and, fJ -LT, NUMl 
FAMP ■ |mP?2,1,J4J1 . 

,AN6. fJ ,FC, NUML 

#!j • AMPfSjf.’syMLT 


i5p(»^ 

jjsr: isiipf!,!? 

I^((»<T .EO. I) A 


I ) 

n 


jur: »5?i. jp.-.u.i fi j, 


FAMP ■ -AMPn ,MUMT# JI 


c 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i n 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ) 

C LOOP 21000 calculates THE RADIUS OF EACH FILAMENT 

C •••• AND DRA68 THE CIRCLE REPRESENTlSe IT, 


NT,:,ji) 


RF ■ SORT (1 ,009/TP( 

DO 21000 A 1, 01 

X ■ XL ♦ RFCCF*8IN(K«PI/2B,1 
'l.T Rf •CMtpSiXAPf/?®,} 

Mr(X,V ,0, 


21000 

C 

C I I I I I I I I I I I 


CON 


V ■ YL • RF«CF«CQ9(X0P! 
IFTK ,EO. 1 ) CALLi CALCM 
^aJ^l CALMP(X,Y,T,1) 


1 ) 


I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I C I I I I I I ) 




sO 

uy 


I 


h { 


- 9 IP 


19128199 


C • • • • 

c • • * • 


LOOPS 22000 AMD 2S000 OPAM TMf LJMfS IN CACM THAT 

9HOM Th? CURRPMT TM each FILAHeUt, 22000 FOP I M AMO 23000 FTP PUT, 


IFf ABS(FAMP) 
MLIME ■ IMTfA 


IF(FAMp it, 0.) op TO 203 
5n 22090 L & If NLIME 

iTN • pNcF*&opTn ,.(t .. 


9) BO TO 20000 
•10.) ♦ .9) 

TO 20300 


L/(fMLIMr*n/ 2 ,))** 2 ) 


K ■ ML • 10 

y . vr ♦ RF*CF*( 1 ,.L/( f MLIME a1)/2. ) ) 
CALL CALCHP(«,Y,0, I J 


CALL CALCHP(« 

V ■ VL A X 0 


22000 

20300 


y ■ YL ^ PF*CF*(l,.L/((MLINCAn/2.)) 

call CALCMP(*,y, 1 , iT 

comTimuf 

1F(FAmp .6T, 0.) 60 TO 20000 
OO 23000 “ • If HLIME 

yp • Pf*CF»S 0 PT?l,.(l,-M/(fNLlMEAn/ 2 ,))** 2 ) 
X ■ XL ♦ RF*CF«n ,• 0 /? (MLINEaI )/ 2 , ) ) 
y ■ yr - yp 
call calCmp(x. y, e, 1 ) 

X • XL ♦ PFPCMd ,.M/(fMLl*<EAn/ 2 ,)) 

y ■ YL ♦ yp 

call CALCPPiJ^y, I » 1) 

COMT IMiiE 


y ■ YL ♦ yp 

call CALCPPiJ^y, I » 1) 

23000 COMTlMiiE 

C . .... 

^ M i I I I I t I I I I I I I I I t I I I I I I I I I M M t I I t t I I I I I t I I I I I I I I I I It I I I I I I I t I I I I I ) 

20000 COMTIMUE 
C 

CA**P»»P#P« 000 P» 0 A»P 0 *P 0 PP«P« 0 *P*»»F#PP»AFF**FPFPF 0 »A»PF 0 »P»FF 0 PPF»**) 


c LOOPS 30000 AMO 40000 pPAP THE L 0 JITJOE CIRCLES 

C *••• PlT limes AMO LABFLO PE 9PE C T I VE Cv . 

C 

00 30000 I B 1. 3 

on 30000 J • I, 91 

lEfJ ,E0, n CALL CALCHPfx,r,0,I) 

Call* CALCMp(x,y,i, i ) 

30000 COMTIMUE 

c 

no 40000 I ■ 1, 12. 


40000 I ■ 1. 12. 

it : i Sij: 

CT B qOsitHi 
XL B 4,B*CT 
YL ■ 4.9*iT 
< ■ XL - 2,*CT 

jALL^CALCM^rJIy , 0 , I 1 
X ■ xD 
y ■ YL 

call CALC*«Pf X, Y, 1 * 1 ) 


Xf 

J 


ORIGINAL 

OF POUR (^UAuH 



!5-0CT-l? 


40^09 

c 

c 

c * • • • 

c • • » • 

c • • • • 


“LT • 2'. *(1.12 

?N ■ <L ♦ , y*tT ♦ ,1^9 

TM ■ ♦ ,P*5T - ,07 

^“1-* i5l|9) YN ■ YN - ,14 

C4LL )T<l, YN, ,?l ,MLT,90, , .1 


LOOP S0000 ')R4**S TuO Pf PBESENT4TT Vl 7ILA“ 
ONf 8 h 54IN6 CURPE^^T IN, THf 3TH?#, CUPREN 
i BiOIuS OF 400000 M^TlRft AND 4 CURRENT 0 


«. a 0 1 p 


10P20 

•'0030 


81000 

80040 


RQP ■ 400000, 

DO 50000 I ■ 1, 
GO TO (90010, 
n ■ 8,3 

G^ To'iftaia 

n • 9,} 


80020) I 


82000 

80080 


5 

Y ■ Yu •, ROP*CF*cfi8 I J*Pl /28, ) 

iFtj ,ro, 11 c*ll 
Call* f I » 1 ) 

COnTtQjE 

GO TO (80040,50091) 1 
00 82000 J ■ 1. 10 

t0 m RDP.cksftPT f 1 ,-M ,• J/9,8) **2) 
F ■ YL • *P 

Y ■ Y* ♦ PQP*CF* ( 1 ,.J/9,9) 

Call* CALCMPf », Y,0, I ) 

Y ■ XL « YH 

Y ■ y[ ♦ P0P*CF*( 1 .•J/8,8) 

CALL. CALCmP( Y, Y, 1 , I) 

CONTpJE 

GO TO 90000 

DO 83000 J ■ If ii 

Yh ■ P0P«rk5fiPT(l,.(l,-J/9,9)»*2) 

Y ■ YU ♦ OOP*CF<*M ,»J/k8) 


J C ^ J " It 

Ym • PDP.cksftPT f 1 ,-M ,• J/9,8) **2) 
Y ■ YU • YP 


8 3000 
80000 
C 

C #••• 

r 


Y ■ Yf - Y0 

Caul **i.r‘^P(»f t, 0, i ) , 

Y ■ YL ♦ POP*CF* M ,-J/8,5) 

Y « YU ♦ YP 

CAUL CALC“Pf Y, Y, 1 , 1 ) 
CONTINJF 
CONTINUE 


CALL 9YM90L(6 
Call 9ym93L(9 
CALL 9ym95C(9 
CALL 9ym90L(9 

»coof ■ Ncoor 


OL(6,7.-9,07,,20,22HOIfT0IIUTlO' 

3L(959.,9jB,..20,UHAri8N!O tuRl 

5 l 9 ,k*l, 9 » 

OL I9.8.0.k .?*» 3H0UT,0a,, 


PA9E 4 



9»sEeri0N9 

ICrt Hl9 

^1. 


1 



) 


I 


• FIELD »99,»92) 
9,91., |&) 


t 9I2BI99 


19 «aCT -02 



FORTRAN IV-PLU8 V0?*B1E 
AMPlt.FTN /TR|*LL/WR 


CALL NUMBER f S. 07 «.HtRp ODE, 90 .( 0 ) 

CALL number 1 9 ^ 9 ., 9 $ 9 .,, I 4 , CL 1,96., A) 

JH|j_^N[jMBERi999:.9O9;,tl4,cCJ,90.,0i 

CALl'*NUMBER( 990 '. , 999 .,'. 14 ,RUMT ,9 0 , , 0 ) 
Ruml«numl 

call number ( 099 ., 999 .,. H,RU'<L, 9 0 ., 0 ) 

Call paub 

CALL CAlCMP(*,V, 1000 , 2 ) 




PAOE 9 


vO 

00 


) 



^ (A 


r‘- 


rj 

CO 


FORTRAN IV-PLU3 V02-B1F 15l29|B8 IB-OCT-B? P«0E< t 

curplt.ftn /tRiXuL/tr 

********************** •••***••**••*******••*«*•**••*•***•• •*****«) 

i *••• CDRPLT SHOWS ThF CURRENT VECTORS IN THE IONOSPHERE 
C 

? 

0081 OlMfNSION AMPf?.B,72), TP(4,B,72) 

0002 LOGiCALwl OARRivT I t ) ,ENUH{?) 

0001 REAL*4 Mfl . ^ . 

0004 DATA DARRAY / «n<, 'S«» t0'» 'At, <T'» 

1 0 / 

0005 data pi / 3,141892BB f 

? 

0000 NRITE (6,01000) 

0007 01000 FORMAT (' ENTER FILE NUMRER|NN I) 

0008 read (6,01010) FNuM 

0009 01010 format (2A) 

00)0 DARRAY(e)aFNUM(t) 

0011 OARRAYn0)*FNUM(M 

C 

00)2 WRITE (6,01020) 

00)1 01020 format ( T |ntN current MAONITUOE- ' ) 

0014 read (6, a) cur 

00)5 OPEN (UNITa1,NAMEbOARRAY,TYPE»(OLO») 

00 6 read {),*} NCOOE 

0016 N ^Nuif* Shlf/buML 

002? CLOSE (UNlTan 

§§13 ' m 

0025 XORG ■ B.6 

0026 YORG ■ B.S 

0027 CALL CALCHP(X0RG, YQRG,0, 3) 

^ 

^ •••* LOOP 10000 DRAWS THE VECTORS AND THEIR HEADS 

0026 DO 10000 I ■ 1, NUHT 

0029 DO 10000 J ■ 1, NUML 

C I II II 1 1 M 1 1 1 1 1 1 M II I i 1 1 'll 1 1 M 1 1 1 1 1 III II i I M i I i I i I ii i I M 1 1 1 li 1 1 1 M M i li I'l s ) 

c THIS PART DRAWS THE VECTORS, 

0 

0030 F ■ (CLI ♦ (2.AI-1 ,)A(CL2-CL1 )/(2,*NUHT))/30, 

0031 XL ■ -4,B*FaC08(2,*( J«,B)*RT/NUHL1 

0032 YL ■ .4,5*F*S1N(?,aT J-,8)ap!/NUML) 

0033 X ■ XL 


O O 

-l ; 
t 


L 


v£> 

<£> 



I 

t 


am 


0034 

003s 

0038 

0037 


00: 

00' 


0040 

0041 


004? 

0043 

0044 

0045 
0048 

0047 

0048 

0049 

0050 

0051 


r J005? 

■o 


0053 

0054 

0055 
0058 

0057 

0058 
0050 


0080 

0061 


0063 

0084 

0065 

0088 

0087 

0088 
0089 

0070 

0071 
007? 
0073 


/HR 

Y| 


1 5I29ISB 


15-0CT-B2 


P40CI 2 


IF 

XF 

YF 

X 

V 


Y,g,n 

• *2 

r , 

[♦CUP* ?AHPr 




♦ *»<P(2,1,J)**2) .to'. 0.) CO 


TO 10000 
XLP*2^VLP*21 
X[P*2* YL'**2i 


C4LL C*LCMP(¥, Y, 1 , 1 ) 

1 1 1 1 II 1 1 1 i M 1 1 1 1 1 1 1 i i 1 1 1 1 1 1 1 1 1 1 1 1 1 li 1 1 1 1 1 1 1 1 1 *1 1 i I i 1 1 1 1 1 II 1 1 1 1 i 1 1 1 1 It I M 1 1 ) 
THIS part draws the apron HEADS’. 


* • • * 


th ■ •Acnsr* 

TH ■ ATAH??( 
X ■ XF c , 

Y ■ YF ♦ . 
CALL* CALCH 
X ■ XF • , 

Y R YF ♦ ; 
CAtU CALCM 
X ■ XF 


(YL«YF)/80PT 

♦cosMh • 
x,Y*l<t) 

H ♦ 

H ♦ 

«1) 



2* 
(X,Y,1 


nXLwXF)* 


,*PI/100,1 

•pi/is0; j 


•2t(YL-YF)**2) ) 


Y ■ YF 

CALL' CALCMPfX, Y, 1 , 1 ) 

C . 

C I I I I I I I I 1 1 I I I I I I I I I I I I t I I I I I I I 1 1 I I 1 1 1 1 1 I I I II I I 1 1 M I 1 1 I I I I 1 1 1 1 1 1 I II I I 1 1 1 1 1 ) 

10000 CDNTIHUE 
C 

I 

^ •"** ’ULT^LlSES^IwD^LAllECs^R^SPFcKSrklY^*^*^'^^ CTRCLC8 AND 
^ DO S0000 I ■ 1 


DO 50000 


000 J ■ 1. 61 

FfJ .ro, 11 CALL calCm 
ALL' CalCmP(/,Y,T,1) 


pf X, Y.0, n 



00000 CONTINUE 
C 

DO 60000 

I? : 

CT ■ . 

XL D 4,5*CT 
YL ■ 4,5*St 
X ■ XL- • 3,*CT 

Y ■ YL • .3,*ST 

CALL CALCM^fX, Y,0» 1 ) 
X ■ XL 

Y ■ YL 
CALL C 


St 




ALCMPfX, V, 1 • 1 ) 
,7iCT ♦ ,106 


2.*n-i 1 

([' ♦ ,7*CT 


o 

o 


G 

• a 

i? 

2 a 


FOUTRAM 1V-PLU9 V09-91E 15l?9lSB IB-OCT-0? 

CURPlT .fTN /TR|AlL/P» 


000RR CON 


i N ■ YLi ♦ -7*9T - ,0F 
F (MlT -of, 9) YN ■ YN - ,14 
*1^^ *ljMftpRfxN,YN,,21,MLT,0»l.,-l 


CALL number ( 990.. 999. *.14, CL2, 90., 0) 

RUMT ■ NUMT 

CALL NUMBER (999, *999.* .14*RUMT*PB.#e) 
RUNLrNU'4L , 

call NUMBEA(999..999., .14*PUML'*O0.,B) 
call PAU9 

call CAlCMP(X,V. 1000*9) 


CAlCMP(X,V. 1000*9) 


2 o 

*c o 

O 2 

-O r- 

•p V 
f - x> 
>-• o 
n: 


o 


1 51 3 (f 120 


I 5 - 0 CT-B 2 


fgHTRAN IV-PLU 8 V 02 -BU 
BRKilC.FiN /TBiALL/KR 

BRKALC CALLS MiCMOD TO FIND CURRPNT OENSI 
FIELD C 0 ' 4 P 0 NENT 5 OF tHE BIRkIlaND CURRENT 
BY COROTS AT POINTS ON A CIRCDlAR ORBIT, 


PASEi 1 


0001 

000 ? 


0003 

0004 

0005 

0006 
0007 
000B 

0009 

0010 
00! I 
00 ' 
00 
00 


; I 

li 


0015 

0016 
001 7 

0016 

0019 

00?0 

0021 

00 ?? 


0073 

00?4 

0025 


0026 

0027 

0028 

0029 

0030 

0031 
003? 


* • • * 
* • * * 
* • * • 


TY AND MAGNETIC 
MODEL' defined 


|jOGlCAL*J OARRAYI 1 () .FNUM(2J 

datIi* DA ' I ' * • 8 ( f • i I » ( D T » • A i » <7 

10/ 


AMPfP-S, 72 ) , 


in, 1 8 ( f • n * i D i » • A i » i 7 • « • f • « i if I I 


8ijJ ?’llSI§j'// 

DATA PI / 3 , 141 B 9 ?Ab / 


RRITE ( 5 | 01000 ) 

01000 FORMAT (' enter file NUMBER|NN ») 
REAg^jB, 0}0I0) FNU« 

RRAY 


01010 format f?A) 

DARR 4 Y(i)«FNUM(i ) 
DARRAYnBI^ENUi^f?) 


1 1^010) 

^NUMi 


(5,*) alt 


NmE AS 

01040 format (f^HTER 


FOR TH 
FOR 


WRITE fB|010201 
01020 ^NTER 

C * 

WRITE (5, 

01030 format (' 

read (B,*) 
WRJTE_(5|01040) 

1 
1 
• 

4 I 

read (B,*) IFLO 

c 

WRITE (3|010B0) 

010B0 format (1 ENTER 

i : 

I . . ...RfiR) 

0 1 060 

^ .EO',0) 

!**PI/1 »?. j 


A|.TITUOJ|tJ^^‘-^'^*^*°^» and theta I) 


*jNTfi NUMBER OF MEASUREMENT POINTS •) 


FOR 

FOR 

FOR 


FOR 
FOR 
)R 


THF FIELD OF. ALLI C 

wlvS:jtj?3^?N5hj 

EAst.WEST ONLn/ 
RING CURRENT ONLY*) 


URRENTBI / 

i / 


4 FO 

IPASS 

IPASSbI 


f^oSGAf"??!' 

IF (IPASS, I 

ST ■ SlNfTHj 
CT M COS|tHli""rj/|nw«# 
SlNCL ■ SlNf INCL*Pt/|i 0 , ) 
CINCL ■ COShNCL*PT/l 80 . J 


POLAR 1/ 

eoDatoriali WEST'/ 
EUUAtoRl aC« east I ) 


o 
r j 


o o 

3 

“0 52 
o S 
o > 

aj r- 

O '3 

c 

£I 

-i ^ 
-< j* 





FORTRAH 1 
8R<il c.rt 


V-PLUS V0?*B1E 
'V /TR|ALL/WR 


16130190 


13-0CT«B? 


PAOE< ? 


(NJ 


0033 

0034 

0035 

0036 
0039 
0036 

0039 

0040 

0041 
0049 

0043 

0044 

0045 

0046 
0049 
0046 

0049 

0050 

00?l 

0053 

0054 

0055 

0056 
0059 
"056 
0059 
0066 
0061 

006? 

0063 

0064 

0065 

0066 
0069 
0066 
0069 

0090 

0091 





WR 

WR 

NR 

HR 

NR 

DO 


NtTNl.NAMEaiHAG.DAT* rTYPEB'NENt ) 

‘ ncooe.alt 

INCLiTHCTA 



I L • I nt 


h 


0000 


NR 

CONT 


-n^A5S.n*?90.*Pl/IB0, 

SMP • lINfMPj 
CMP ■ C0S?MPi 
XL ■ (RE ♦ AlT)*("8TaSMP ♦ eT*C 

It : SI : Jtt 

CA||j^MAeMOp{XL,VL»ZL»PLD,IFLO) 


• n)*Pi*4e'./i60, 


'•'ssjfwev’ 


(!,•) PLft 


2 Q 

.*n ^ 
•r o 

o J 

o . 

^ r* 

0 T3 
CT 

1 

EZ "3 
H 

■< 03 


009? 

0093 

0094 


CLOSE (UM1T»1) 

STOP 

END 


o 


C«J 

CO 




LL/-^b 


1SI31 197 


lB*0Cr»B? 


P* 8 E 1 




R0«? 


0003 

0004 
0009 
000A 

0007 

0009 

0009 

0010 

2®1 • 

0013 
0011 

0014 

0015 


001C 
00 7 
00 
00 
0000 
00T 
00; 
00: 

00?4 
0035 
00?6 
0037 
00?8 
0039 
E030 
0011 
00 12 
0011 
0014 

0019 


C 

c ^ * 


9UBR0UTI*JE MAGMOO f XL » YL » ZL , 7L D » I PLO 1 


CALCULATES rURRENT DENSITY AND NACNETTC 7IILO CONRONENTS 
OF BIRkEL*N 6 current model defined IY CUROll, 


o 


U /FILE/ CL1.CL2iNUMT.NUMl>RF.TF(4.fl,721,AMRr2.9,72), 

B ,t ),RB( 2 .fll,ftT(f», 6 ) .RM (BT.Rz|(9j,ftEJ(ll ,RI (3,5), 

DIMENSION rkH^(6,H), FLO(4) 

OlMENSlON 9F(i'l, BTF(3), 6tl(3), BTE(3), BTR(3) 
5?*J*^,L4M$D4, MP, Jt . 


DAT* PI / 3,1'^ 

DO 10000 I m 1,3 

flTMn ■ 0 , 
BTI fj - 0 ; 


10000 


It 00 TO 03000 

If nFCD.QT.2) GO to 03000 


; h 8 ?E ?5rf.5”SSo’l!'./l5liUlilS??.'’t5«?F.?}! ’“"-LriN. 

> POSITIVE CURRENT IS VERTICAL'. 

1 ■ numt 
00 2000 A 


NUM| ♦ I 
2000i M ■.1,2 
IF (M .EO; 3) I ■ NUMT 
on 20000 A ■ 1 , I 


DO 20000 J ■ 1 


FampTNOMt^I , J) ■ •*mPM .NUMT, J) 

Ifl 3 : S ?; [ !”. 5 i 


2001 0 


20020 


F (F*A^(ft,i| 

GO TO 20020 


.EQ. 0.) 60 TO 


• *MP(1 ,N« 1 , I) 


FAMRfN.l),- AMPf2.N,NUMLl • *MP(2.N, 
(J aBT. I) FAMRt^.Jl ■ *MP(2.N,Ji 
If (FaRptn.J) .EO, 0) 60 TO 20000 

i AMflDA p 2. • f J-,5)*Pt/NUML' 

^ T** .fO. ■ 2 ,fJ«PI/NUML 

L* ■ BinIlambda) 

CL* ■ CPSfl AMBOA) 

RMfi,!) p CL**CMU(M,N) 


li 


• *MRC2,N,J) 


ORIGINAL PAGE IS 
OF POOR QUALITY 




-PLUS V 


f5;5iE, 


ii 

PMf 3»3) 


t5l3t 197 
- 90 U(M,N) 

'III 


1 9«0tT-8? 


»eP0n CONTINUE 

c 


»HM,n ■ Cf4*JMUfM,Ni 
RHM.n ■ 9 LA« 9 HUM«NI 
PMfJ,3) ■ C«U(M,N) 

i mm l itisjiiili : 

RC^ ■ *F«*? ♦ YF«*? 

R0F ■ SORTIRCF ♦ ( 4F-RB ( M, Nl 1 
RTF ■ SORTfRCF ♦ C ZF-RT ( M» '<) ) •*2) 

RXP ■ f AMP( n. J)*(TANH(RCP»TP rM|N. J) ) )*(YF/RCF) 

» 1^2 ♦ gYK*RM(?,IF) 

STFeiF) ■ BTF(I^) ♦ BF(!P) 

?F**nii''(M,H, J)*RCF) ,GT, 49.) 00 TO OBBOB 
Jt • Famp/n JUtR(m,n, jJ 

/(PI«nCOBHf.TRCM,fc, J)«RCF) )««B) ) ♦ JT 


IF (1(1-0. C°.?) CO TO B7BBB 
B3BBB IF {iKd.F0.3) 00 TO B4BBB 


00 30000 H ■ l.NUMT 

DO 30000 J ■ l.NUNL 
LaMBOA ■ 

K« ■ BINfrAMAol) 
CLA ■ cOS(LambDA) 


)*PI/NUMt' 


m 

RM(3, 3 


If lifts 


CLA*CCL(2,N) 

9LA*CCU2,N) 

feiSI*® 

fO.=3) 00 TO 3B1BB 
PQ.4) CO TO 3810B 


’’nri 



fO^TBAH IV- 
( M 40 vinO>fN 

r 

E 


OJ 

<n 


PLUS V 


?> 9 J 7 

007 S 
0079 
0000 

0001 

000 } 

P 003 

0004 

0000 


1007 

0000 

0009 

0090 

0091 
009? 


0091 

0094 

0095 
0090 

0097 

0090 

0099 

0100 

0101 

010 ? 
0101 
0104 
0 05 
0106 


u 


000 

010 


15131 167 


* * • • 
• ••• 


IPJ!,?5rce86?N?"5.’‘S5J?i;32KJr 


15-OCT-B? 0ABe< 3 

lONOSPHE.Vie CURRirNTS'. 


lF(*MPn,N,j) ,£Q, 0.) 90 TO 30010 


o 

(T 


tf 

YF 

IF 


XL*RHr 

XC*RM' 

XL*RM 


iiil! : mm : ilti31lisiil5s!l:li 


in 


YF 4 *? « ir *»2 

$O?T?0^f ♦ j*F 


ryf 

BZF 

00 


SORT^RCF ♦ (XF - R 
-R 








H 

,Nn/Ri 


(if!> ■ 


•|ZF/RCF) 

•>lvFi/RCF) 

2 , 


• • « * 

101 00 


fgi?.? 8 !'cB§l?J" 5 .'«lT 3 n 5 . 


FLFCTR0JPT8, 


IF(4MP(2.n, J* ,C0* 0.) 00 TO 30000 


XF 

YF 

7F 


• RH 
XL*RM 
XL*RM 




ilili : mm l 


RCF • XF**? t (IF - RZFlNn**2 
RU ■ SQRTfRCF ♦ tYF ♦ NEJC'I) •*?) 
RE? ■ SQRt^RCF t (YF - REJ(Nn**}) 


BXr • *MPf?- 

55 


32000 
10000 rONTINIIE 
f 


*XF* f (YF 4 rF jI 

no 32000 iP ■ 
roMlNul 


bi-riFtiir*""'’*"’ 


Bxli^ 


) ■ PTE 


0107 IF (IFLO.LE,!) go to 04000 

0100 CP»»#p5»»#i»p5M#>pl#S#*i###JSS5< 


0.09 


C 

04000 


IF (IFLO.EQ.-B) 60 TO 07000 


OKKilNAL PAGE 13 
OF POOR QUALITY 



F09TRAH IV-PLUS Vt»?-5tE 


t9l3t ISF 


15-OCT-0? 


PAOC' 4 


a\ 


Sll? 
W 


9\\t 


Ptia 
(?! 3'. 


0j 3; 
0i 3a 
e 39 
014a 


0H1 
01 4? 


• • • * 


LOOP 
DO 4 


40000 calculates flELO DUE TO SING CURRENT 
I»1,NUML 


^AMBDA.P.i 


ILAaSIS 
CL AbCOS 


(LAiJol® 

aAMBDA^ 


9>«P1/NUML 


CLA 


RP 

RH 


SLA 

e< 


""'iilliia 


RM 
RU 
RM( 

RM 3, 

RMb, 


/?}* 


0 . 

0. 

?; 


:f f ampr(I)',fo,0.) go to aseia 

Hiii : i \[ 


XF«aL*RM 

YF«XL»RM 

2F«XL*RM 


• RM 

rRM 

»RM 



?7 

RCFi 

« f XF®RI 

0!l 

iB 

RRji 

ij 0RT(A 

01! 

)9 

RR2i 

•SORT (R 


Cf ♦! Yf ♦« 




•?) 


§iii 


00 41000 

? HI 

a! 33 

41000 

CONTlNjf 

01 38 

40000 

CONTINUE 


J * * * 4 

LOOP 70000 


BXF«APpR(n*(TmMfRCf4TPRn)l WrF 

M (yf mgsi/rAi F(niNofl.vn/RR?)/Rcr 
BZF»-AmPR( )•/ t ANH|RC*«TPR( t )) )4fxA«R!NOA) 
*t (YF^PImgS ) /RR1< IRIM68-Yn/RR?)/RCF 

■ BXF*R''(1» IP)tBZP*R'<T3. IP) 

')»B7R f IP/ ♦BF (IP) 


° ff8?|p! • 2iA?*OTF(IP)48TI (IP)4BT| MP)*8TR(IP) > 
F0000 CONTTNUf 

FLTmA) ■ JT 


RETURN 

CKO 


O 


Cr I jL/{» 


FOi^TR^H IV-PLUS V02-81F |B|32|40 t5-0CT-8? P*8f 1 

( ' BR4PI T.fTN /TR|4LL/»<R 


0001 

300? 

0003 

0004 

0005 
0008 

0ear 
0009 
0009 
001 0 


00 

30 

00 

00 

00 


001 8 


0017 

0018 
0019 


03?0 

0021 

0022 

0023 

0024 
0029 
0026 
0027 
0029 

0029 

0030 
00 31 
0032 


0033 

0034 
0E39 
0038 


0037 


* • • • 
• • • • 


BRKPLT plots CUPRFNT density *N9 FIELO CO'«PONENT9 
FOUND BY BRKALf, 


DIM 

iImfnsIon 






NSION 
NStON 

5IUU?iin«55|f, 


FTEM(3), FACTf4), FL0T4) 


I 0 / 

DATA PI 

DATA RE / 8371000, 


MP, 

<0» , 


NA 

I 




•T*, • 


• 1 


I I 


C 

01000 

01010 

r 

01020 

C 


J 


DATA ALT 
DATA ate 


3-1419928B / 

8571000. / 


/ 140000, / 

* / 171, Iyi, 'Z», »J«, 

•N«, *E'» *Vt, • Jl, 

•S«, «D*, *V»» IJ* / 

write I6|O1000) 

format (> ENTER file NUMBER|NN«) 
READ (8,01010) FNUM 


FORMA 
OARRAY ( 
OARRA Y ( 


Pl 

10) 


NUM( 1 ) 

• FNUM(?) 


WRITE (8,0 

format (' 

READ (8,*) 


MODE 


1 FOR XVZ) ? FOR NEVA 3 FOR SOV •) 


I 0000 

r 



c •*** CONVERT TO NEV ON 21 SDV ON 3 


o 

OS 


original pauL is 

OF POOR QUALITY 



A. ^9‘»TR*N 


003B 

0039 

0040 

0041 
004? 

0043 

0044 

0045 
004S 
004T 
004B 

0049 

0050 


0053 

0054 

0055 




1 51 3 ?l 40 


DO 11000 J. ■ I , NMr AS 

1 •nAAis.h 4 ?; 0 , 

SMP a SINfMPJ 


/(NMEAS-1 ) )*PI*40./180. 
0,«PI/1B0. 


CMP ■ CQSfMPi 

XL ■ -ST*SMP 4 CT*CMP*8INCL 
YL ■ CT45MP 4 ST*CMP431NCC' 
ZL ■ CMp*rjNCL 

I il! : Itiyi 

W : 11 !$-= : 


FL(3,J) ■ .FTFMn)4XL - FTfM(?)*YL » »TrH(3)*ZL 
IF (HOOF '.EO'. 3) 60 TO 11010 


11010 F 
1 1 000 CONT 


CO TO 1000 


**** FIND maxima 


nil DO 


12000 CONTIN 


I 20 I 0 l.NMEAS 
FMAXm ■ AMAXI (ABS(FL(I 
CONTINUE 

FMAx!i)*FHAx}[)'^^’‘^*^'^^^^^^^^ 

CALL CALCMP(X,Y.2,0) 
call CACCHPix,Y,0,2) 


BS(FL(I,J)),FMAX(in 


4 *** PLOT BACKGROUND 

KORG ■ .15 
YORO ■ ,1 

CALL CAlCMP(XORO«YORO»0«3) 
DO 20000, I ■ 1 .9 

Y ■ I !• 1 ) *1 *2 ♦ . B 
CALL CALCmA^X, Y* 0f 1) 

Y I ♦ I® 

CALL cal 6 m^}x»y« I t n 


ORiGir.’AL PA3!Z IS 
OF POOR QUALITY 



(*' 

T 




'/HR 


t 51 32 I 40 


1S*OCT»02 


PAflE 3 


0077 

007B 

0079 

00B( 


0082 

0083 

0084 

0085 

0086 

0087 

0088 
0089 


0090 

009t 

0092 

0093 

0094 

0095 

0096 

0097 


099 

0100 

0101 

0102 

IISI 

0105 

0108 


0107 

0i 08 
0 00 
0110 
01 1 
0)12 
01)3 
01)4 
01 5 
01 1 6 

iiU 

0119 

0)20 

0121 


in 


9) 

9) 


M I 

XL 


9) CO TO 20000 

N . M • 

IX (J ’OT. 

* ■ XL 

Y ■ fN/10. t 1 - 1)81,2 

CALLl C*LCMP(X, Y,0( ■ 

X • ,I ♦ XL 

Y 5 fft/10. ♦ I - 1)81,2 
CALL' CALCAP(X, Y, 1 , 1 ) 

20000 CONTINUE 
C 


J - 9 
• 10,9 

I# 1 ) 


*888 indicate GEOLAT. 

1 f 


INVLAT, MLT 


DO 30000 I 
DO 30000 


J 

f45< 


5 

» 0 ,. 


SJ ■ l,8j ♦ 5, 

?o II cSbI?! 

G*CINCL 


(IPABS-l )*27B>PI/1BB, 


Jt 


C' 

sSrt f XL 8*2 > YL8*2) 


^8 180* * TACOS (SORT (XL8828YLA82))) /PI 
■ "l8Tl.*(Afoi( (§QRT(RE/(RE4ALT) ) )8C0S(CLATR) ) )/PI 




)) 

)) 


10010 


H 8 6, (G4(IPA8S«1)8270.8PI/1B0,)/. 13962641 

IPTI ,E0. 6) 00 TO S 00 I 0 
H 8 (l.n*?,4 ♦ ,8 

call" i ALCHP(X,Y,0, 1 ) 

X 8 H 

Y ■ ,05 ♦ H 

CALL CALCMPTX, Y, 1 , 1 ) 

GO To 30000 


IP 
H 

X N 

YN * H - .24 

CALL' NUMRf R(XN,YN, ,14,OLAT,0,,2) 


( ABBTAMODTGJ, 10,) ) 
• 8 

I H • .29 


,0T, 0,) BO TO 30000 


»— • 
o 


C 


FUOR QUALITY 




15139140 


! B-0CT.8? 


XN • K • .99 
VN ■ H - ,45 

CALU' NUMBf Rf IfN, YN, .|4,IL4T,8.,?) 
XN • N • .?9 
YN • H - .66 

CALC' NUMBM(XN,YN,,14,MLT,I,,2) 
10000 CONTINUE 


>*•* dram latitude circles 


40000 

C 

c 

c • • • * 

c 

1 

04900 ( 


CNnRHB3lN(40.*PI/100. ) 

IF (IPASS.GT.n 60 t6 §4900 
DO 40000 I • 10,40.10 
CRADaSlN(T*P|/lM.l 
DO 40090 } ■ 1, lAl 

Xb CRA0*SIN( J*pI/B 0, l/CNOR'<t!2,) 
Yp-CP* 5 *C 03 ( J*Pl/ 8 B|J/*NO 044 l A 
iFfJ Jo. 1 call ClCCHP(X,Y,ft,l 
__ ^^Aj.tj^CALN^|x,Y,Tf 1 ) 


OR FOUaTORIAL VIE" OF THE EARTH 


IF (IPAsa.EQ.n 60 TO 05000 

CALL CAiC^P (X,Y.1«*5) 

DO 42000 lRl.9 

YBCPAO/CNORMAi.e 

call calchp (X.y.0,1) 

Ya*CRAO/CNORM«0,6 
CONTlfcuE^**'**^^ fX,t,l,l) 

?m 5 cJSl?!l?!!i^/i 80 ., 

CONTINUE 

CALL CAlCHP fX,V,0,«5) 


0RH41?,F5 


42000 

C 


43000 


•••• indicate time 


05000 IPb! 
IF 


IPASS-GT.n IPr2 
0000 I R 1.4 


c 


PAOE 4 


T 


c r> 

**: I 

i 

c 

O 

:r ■ 


■< c.’ 





90^01 

90010 


R00?0 


90021 

90030 


90040 


19132140 


19-OCT-02 


Y ■ .,29*C ♦ 0*9 
CALL iALCMPfx,^ ,0, 1 ) 

X ■ 9 ♦ 12, Y9 

Y ■ -C ♦ 9.9 

ssbb 

GO TO ^9^|l0.9A020»90030, 90040) 1 

IL ■ LKflP-I)«.9 

dlLL ‘JOM0fAfxL,YL..14,RNU^,0,,.l ) 

Cn To 90000 


9A020»90030, 90040) 1 


XL • p.8? 

"iMsJ^OE'.P) 00 TO 
CALL NyM0fR7xL,YL» .14, 


CALL NUM0E 
GO TO 9000 
call 9YMB0 
GO To 9000 


90021 

RNU^,0.,«1 ! 


(XL.YL, ,14, lHW,0,,n 


XL ■ t2,93fnP-l )*,«01 
YL ■ K974ilP-l)A,6 

call MQMBKRfXL,YL,.14,RNUM, 
GO TO 90000 
XL . 

1^ .f i..PO T9 


RNUM,0.,.l ) 


Call nijmb 
CO To 900 


XL,YL» ,14,RNU4,0,,.l ) 
CXL,YL,‘,14,1H9,0.,1) 


A*** PLOT 0R91T 


09000 00 90000 
1 




BMP ■ SlNfMP) 

CMP ■ COSfMP) 

XL ■ -|f*9MP ♦ CI*CHP*BINCL 
YL ■ CT*9MP A STaCMPaSINCL' 
ZL ■ cmpacincl 


ZL ■ CMpACINrl 
XA-YL/CNORMATJ.rO 
Y« xr/rNORMAe.A 


r f IPASS’.Gf .S) X«ZL/CN0RM^12 29 


60000 continue 

IP (MODE ,0T’. n go to 07000 

c 

C»PPP#PPP»PPPPPPPPPPPPPP»P»PPPPPPP) 

c 

c *•** SHOW X,Y VECTORS 


PAGE 9 



O O 

*T1 *: 

■c 

c ;j; 

O > 


) 


X r* 




.U F0RTP*H 
O 8RKPLT> 


■ PLUS V 




t5|32|40 


lB-OCT»e? 


PAOE< 6 


« ■ 11.58 

Y ■ 0.87 

CALL CaLCMP(Y, V, 0, t ) 

Y ■ 

CALL CALCHP(X,V.|,1) 




e indicate maxima 

c 

07000 IPLT«4 


70010 


70020 


70030 


70040 


70050 

C 

70080 


[PLT«4 

\i 

>0 70000 I ■ I.IPLT 


f pilihi 

FACtMJBl, . 

F (MUjh.G 

r jF^Sxlijlk 

GO TO 70010 

fSS li UFiJIi 

00 to 70010 

IF (INTIFMAxr 
FMAXM )■( AINt 

on tO 70060 

FMAX(I)«AINT( 


r ifmax 

F FMAX 
F I FMAX 
F F'lAX 
F iFMAX 

?t : hi 

N B I 4 


AX(I)*1000000. 

EC«1 


i ! ii i!il II II lisii 


III!: 


?ST(FMAXmt‘,5)/FACTm 


00 TO 70000 


000.) ITYbI 
00.) 1TYp 2 
0.) ifYBl 

il 




YL B 



XL B iTYb Ha, 4 
YL B IS«n*?,4 ♦ .6, 

CALL NUMBER ? XL , YL » , I 4 , FM A X f I) , 0 , , I DEC ) 

FM ■ 0, 


113 


PO<»TR*N 
C ' HRKPLT. 

0257 

025S 

0259 

02A9 

02S1 

02A2 

02A3 

0284 

0265 

0266 
0287 
0268 

0269 

0270 

0271 

0272 

0273 


0274 

0276 

0276 

0277 
0276 

0279 

0280 
0281 
028? 
0283 




1V.PLU9 V 
FTN / 




I ALL/iXR 


1SI32I40 


1 8«0CT*82 


PAOE 7 


02 6/1 

0285 

0286 
0267 
0205 

0289 

0290 

0291 

0292 

0295 

0296 

0297 

0298 

0299 

0300 


1302 

0303 

0304 


X I ■ ^ 0 


YL * (5-l)*2,4 - ,47 

CALL MUMifR ?XL,YL,,14,7M,0.,.n 

M*-FHAX(!), 

XLbITY*. 14f 26 

Yl ■ f5-I)*5,4 - 1,55 

CALI NJM0ER ?XL,YL,,14,FM,0.,IOEC) 


70070 



70070 

(XL, YL, ,14, IBHNANOTESLASfQe., 10) 


(X 


70080 . ^ 

79,,, COnIinuI' 

f 




4,9HH1CROAHP8,90.,9) 

7AAXn)/1000000', 


**** PLOT riELOS 

DO 


)LT«ALTD) 00 TO 80000 


80000 CONTINUE 
C 


!f ??i!e 0.4J lAho,(ALT, . 

DO 60000 J » I.nA^AS 

MP ■ ( J»n*2./(NMEAS-1 ) 

CALL CALCfiPfX, Y, 1,1) 


, J)/PMAK(I) 



v9*#0A#0i9) 

0,#0,?4,,l4,3HMLT,0.,3) 


21 ) 


ORY-7.2 
IF (IAA3S,GE,2) 0RY*6,5 
CALL SYMBflu? f2,55,0RY, ,H,5MORBI7,0.,B) 
CALT * ALf/1000- 

I mil: ’ 

E • NCOOP 
NUMBEr 


Sg'iils. 

I 

lU^<L«NUMLr " ' 

:AlC NUmBeR(999.,999,, ,14,RU'1l,S,,0) 


KM,0,,|7) 


CALL NUMBER(12,.0,. .14,RCOOE,0.,0) 
CALL NUMBER(999.,d99,,.14,CLl,A,,0) 
CA .L NUM6ER^99!,999,, . 14,CCi,0.,0) 
RUMTbNUmT 

CALL NIIMBER(999.,999,, , 14,RU'<T,0,,0) 
PUMLbNUMLi 
CALL NUmBI 

rfld-iflo 


Cn.^i^AL i"AG£ S3 
OF POOR QUALlfY 



1Bi3?|40 


t9-0CT»e? 




CALL NUMBER (9®**. • 999, » ,l4,RfLO,0,»R) 


CALL CALCHP(X,Y.1000»2) 


PA8EI B 




^ fortran IV-PLUS V 0 ?*SlF 

TD 3 RKC.FTN /TR|AlL/KR 


15134117 


1 S-0CT-S2 


PAOE. 1 


1 n 


0001 

0007 


0003 

0004 

0005 
008ft 
0007 

0005 

0009 

0010 

001 1 
001 7 

0013 
00i 4 
00 5 
001ft 

0017 

0010 

0019 

0000 

00?1 

00?0 

0023 

0024 
0075 
002ft 

0027 

0028 

0029 

0030 

0031 
00 1 ? 

0033 

0034 
00 35 
003ft 
00 37 


5 3O0RFC CALLS M/| 
r FTELD COMPONENT 
t CUPDIS AT point 
C •*** THE EARTH'S. 

GMOO TO FJ 

\ 81, I"5pS 

ND CURRENT QENS] 
irkElano Curren' 
ERE OF RADIUS 'i 

ITY AND MAONH] 


C 


o 


dimension FL0(4), FTEM(3) 
common /Flit/ CLl 
RZF(2.fl),ftB(2#ft).ft ‘ 


.... 

1.9 / __ 


yMT,NUML>RF'iTP(4-fl, 7 23 , ANP(2,5, 


F2) 


•D» 


' A I 


Sta re / 

DATA Pi / 
DATA ALTl 


< T ‘ , • I ' » » ' • 


8371000 , 

^•Ml 

/ 140000 . 


01000 

C 

01010 

r 

01020 


^P^N |UN^ jal .NA mEb ' 30BRKC. OAT ItFORNi' FORMATTED' « 

*REAp"f 1*01000) ALT 
FORMAt (l«.F9,n 

?6SBAfMfS;iii 

READ (1101020) FNUM 

format MA) 

DARRAY(4)b^NUM(| / 

OARRAY(l 0 )«FNUMr?) 

CLOSE (UNIT«1) 



2 ) 


CW’.G’NAL PAGt !S 
OF POOR QUALITY 



r^ 


CTi 






1BI34I 17 


IS-OCT-S? 


PAGE' 2 


003 B 

0039 


0040 


iSU 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0090 

0091 

0092 

0093 

0094 
0099 


0096 

0097 

0098 

0099 
0060 
0061 


0062 

SS83 


0069 

0066 

0067 

0068 




OPEN (UNIT* 2 ,N«ME«I 3 DHA 0 
5 TYPE-Ine^I ,RfCOR'i 9 I 7 E* 4 ) 


. 0 AT'«ACCC 86 a>DIRECT', 


R^ppr_» NCODE 

HR 
DO 


^?E ALT. 

0000 16 ■ i, 48 
M ■ .?Ip*e 0 ,/ 47 . • 
CM ■ |!n(Gm 1 
CM ■ cosrcH) 


RCODC, RPLD» or 

40 . . B 0 ./ 47 . )*PI/|B 0 . 


DO 10000 IE • 1* 46 
6C ■ (|E«»80,/47, 

: hii ! 

?h : iSf : stmiy-' 

IL ■ ("t ♦ ArT)»CfiCoCO 
CILLI MAOMoD rxC»YL»ZLi 


- 40 , • 8 B,/ 47 ,)*PI/ 180 . 


OCftSCM 


ftCQM 

FLO#IPLO) 


*•** CONVERT PROM K.Y.Z TO 9 #D,V 


imw : ftSiii 

OQ ■ 909T(xL«*P i YL*42 ♦ 

im\ I \mm : 


ELP*2) 


TEM(3)*ZLi)70O 


NllltE T 2 <H) Pl{ 
10000 CONTINUE 

CLOSE (UNIT*?) 

C 

CP»PPPPP*0PP0P0P*PPPPP090#0#PPPI 

C 


1 )*48 


0089 STOP 

0070 END 


T 

O 

C 

^ f 


lO 


i. : 



FORTRAM TV*PLU9 VI02-BIC 19|34|4S t9-OCr«0? 

3 D**LT.MN /TR|*LL/WR 


PAOf I 


0001 

000 ? 

0901 

0004 

0009 

0009 

000T 

0008 

0009 

0010 

001 1 
001 ? 

0019 
0016 
00 7 
0019 
0019 

00?0 

00?1 

00 ?? 

00?1 

00?4 

00?9 

00?6 

0007 

00?B 

00?9 

0010 


3DPLT PLOTS APRAYO IN 3-0 PfRSPfCTIVE 


• ••* 9UBR0IITINF CRFAO READS AND CONDITIONS DATA 


SUBROUTINE' CREAD ( M A T VEC • NRP . NCOU» NERR ) 
dimension MiTVEC(?304),RATVEC(4y 
LOGICAL*! FNAHf]4) 

R|AL*4 H4xV|C#MfNvEC 
COMMON /FILE/ 


NFRR 
IT( 


'OP 


/FIL 
0 

NRR*NCOL 


VFrMULfZMN,FNAM,*:«XVEC,MINVEC, IVEC, ALT.RCOOE 


DO 


1 


, ITOP 
0 


10000 I ■ 

MATVFCU) 

10000 CONTINUE 

nPfNfUNf !■?, name *F NAM, TYPE«* OLD' » ACCESS* <OIRECT») 


I 1 000 


IF (FNAMf 3) ,E0’. 'Mf : >11 TO 01000 


AlAf)?2^jT ) 

RATVECll )»RATV e 6 (i)fZMN 


il 


jjONJINU 




ATV|C( J) 
0 01010 


RATVECmiS, 


01000 READ (?'l ) 

MAXVrC«0. 

MINVFC-0. 

C 

DO 1?000 

ftEAO JZ'J) RaTvEC 


ALT,RC0DE,RUMT,RUML' 


AO 

maxvei 


hinvec«amin 

M*TVEC(J 


12000 CONTINUE 
C 

01010 RETURN 
END 


AMAX 

IN 
1 ) 




CRlC'.r'iAL P"GE 13 
OF. POOR QUALITY 


IV-^LUS V0'-6lE 
N /Tft I 


1Sl^4lSn lS-OCT-8? 


PAOe 3 


r* MAIN PROGRAM 

Dimension ocas.ar) . ym(20h) ,xm(7B0) , vmr(2R0) ,kmp(?0r) 203), 

FNAM(| 4 ).NAH( 6 ),NU« 4 (?), 0 AT(B) 

REAL*4 MAXVEC- MiNufr 

?WE 8 !! 5 » 55 E >’38655 

isiiuitHE iri;!; isini, 

I fssss’fijsis'ji 8f :vjr,EiaEii’i?5S!W) 

RRITE ( 9 , 80000 ) 

I format ( I ENTFr FTLF »-AME IAAAAA 4 NMl) 

REAP ( 5 , 80010 ) MAM, Rum 
I FORMaI ( 6 A, 1 X,?A) 

FNAM( t 4 ) *0 

write ( 9 , 80020 ] 

I format (' enter VECMUL, ZMIN, XVIIW, rvicw, IVFC n 
READ (9,*) VECMul ,fMN,kV,YV,IVEC 


800Q0 

80010 


00020 


80030 


C 

c 

c 

09099 


WRITE (5,80030] 

FORMAT (' ENTFR FQC AND 019 I) 

READ (a,*) FnC.DI9 

OPFN(UNiT«l ,NAHp«i 30PLT,0AT« ,TYPEw»OLO*,FOR'("«FORMATTEt)i ) 

• initialize plotting parameters 


DIN 

UN1T(3) 

CALCHP(0,,0.,2,01 
CALCMP 0 ., 0 ., 0 , 2 | 
CALCMP(0.,0 ,0, 3J 
CALfMPh.,.^,4,-4) 


00001 

80)00 


VY3«X9C*(NX.O*SAj a9A2/(902/FOC.J ,) 

1 ♦YSC*(NY-1)*CA{08A2/((8024804)?FOC»I,' 

VVbVX2-(VA2-VX4)/2. 

YV»yY3/i, 

WRitF (3,00100] MAXyEC. MINVEC 

format (' MAXVEC ■ •,1PEF,1,TnT MINVEC ■ ',1PE0.|,IST') 



K) 

O 


. ..JINAL PAGE IS 
OF POOR QUALITY 



FORTRAN IV-PLUS VBg.’JlE 
f ' 3DPLT.FTN /t6| AlC- 


0 01 
0 04 
0 0S 
0 0fl 
0 0F 
0101 
0 0; 
0 0 

2 ‘ 
§' I 

0 14 
— 3 15 

<-n 0 i 6 

o 0 j 7 

0 10 

2 I* 
0 ?0 
0 ?\ 
0 ?? 


01000 


NXMP«N 

*|e‘ 

? nvi.OE.NV) 

60 TO e?000 
IXI«IXU1 . 


15134150 


cn TO 01010 


Xj.LT.NX) 00 TO 07000 
PAu| 

CALCMP(X,V, 1000,71 


C 

02000 


P701 0 


01010 IXI«IXU 1 . 

IF (fxi.LT.NX) 00 TO 02000 
Call piOB 
CALC PAul 

JAL^ CALCMP(X,V, 1000,71 

c 

02000 IXhIXI 

IybIyi 

|xn»i 

IxPrio 

02010 call PTRAnX,IY,Xl.Yl) 

IP (IXD.LC.fl) 00 to 07110 

IXTptXtl 

CA||L^PTAa(IXT,TY,X2,Y2) 

Call ^JtJ a;!XT,IyT,X 3 |Y 3 ) 

V TbT 1 tIX?-X 1 j*lYS.Yl)/f X 3 *xi ) 
ip (Y 2 ,Ct,YTi 60 TO 02100 
CALL MaOMX, lY. INOG) 

IF ( 1 NOG.LP. 0 ) 60 TO 03000 
X 2 *X 3 

i 2 »Y 3 
XaiX7 1 
XDB 0 
YDb -1 

0 TO 03000 

YT-Y 17 CX 4 -X 1 J*(y 7 .Yl )/(X 2 -Xn 
IP (V 4 ,LP.'<tl 60 TO 03000 


C 

02110 


IP (V4,L 
CALL d!I 
IP (!NOG 
Xl«X4 

SJ-?J 


4,LPvn 60 TO 0 
Dflpnx.lY-lNOO) 
NOG.dT.M 60 TO 


03000 


03000 


YT»IY«1 

ALL PTAAnx. IYT.X2. Y7) 

P (CAl .Of,. .04201 66 TO 03000 
XTBlX.f 

ALL ?tAa(1XT, IYT,X3,Y3) 

TbYI ♦(X3-Xn*(Y2-Yl 1 /(X 2 -Xn 
P (Y3,LE-,Yt) gn TO 03000 
ALL OlAOnxf ,tv. tNOO) 

P (INOO.GT.P) o6 TO 03000 
2*X3 


X2«X3 

Y?«Y3 





'/HR 


lBl34tSfl 


I5-0CT-8? 


P*8C ft 


43 

44 


0145 
01 4S 


4? 

48 

49 
90 


0151 

0}52 
S3 
94 
99 
9ft 
57 
9ft 
99 
60 

§1 

93 

64 

65 


01ftft 
0167 
0 68 
0(69 
0170 

01 M 
0 72 
0 73 

0174 

0175 

0175 


03000 iBl 

03001 !■! 


:r nrj'Mzwii ?8 JS.SJ 

03010 C2»(V2.»1 )/(X?-xn 


03010 


03020 


CiBXi«C2»y 1 

tee; 

XDFNBCr'-eS 


X'4Pf I-D) 


IF (’ 00l6*B8f XNUM)’.GF, AB9(X[)FN) ) CO TO 03100 
XT»xftUM/XOFN , . 

(XT.Xll.L5..00n,CO TO 03100 
60 TO 03100 
1 GO TO 03100 

XT.OT.XflPd )! 


03100 


(xt-xj) .0 

^xt.ot.xAp 


0011 
60 TO 03100 


03110 
03110 


0 04000 

, i*2,LP'.XMP(I)l GO TO 
> d.CL^XMR) c6 TO 03 

;g TO 03020 
03110 XfaX? 

Yt«Y2 
C 

04000 XTT«.6*(X14XT) 

|F (XMPni.OT.XTTl 00 TO 04020 
IF n .L'T.nJImp 5 00 TO 04010 


04010 


C 

04020 


YTTByMP(iw(xTT-xMPni )4jrvMPf T.n- 

IF ' .P*(Yl4YT) .LT.vtf ) 00 TO 09000 
N. ■ •NXPLit i 

NXPtitP.?®01 GO TO 04030 


YMP(i) )/f xMpn-n-xMRd) 1 




o ^ 

.m X 

T3 ? ’i 
O :2 

O 5 

50 r- 

•o ^ 

C > 
T» O 

E oi 

■H — • 
-C <Jt 


c 

77 

04030 

XPL (1,NXPU1 

■ XI 

0 

78 


XPL (2,NXPL) 

■ xt 

0 

79 


yplm.nxplJ 


0 

60 


YPL <2,NXPl.i) 

■ vf 

0| 

181 

05000 

IF ((XT-X21 

,ce’. 

0l 

1 82 


X1«XT 


•.001) GO TO 05010 


C 

05010 


n 03000 

?:i?8 


?b1 YO 
YDb-IXD 
X0»-1T 

F ((Ix*!X0).GT.NX) 


00 TO 05020 




FORTRAN IV-?I.U 3 V 
30 (»L'T.FTN / 


/WR 


1 5l 34 ISO 


15-OCT-e? 


PAGE 1 


i\v, 


95 

9S 

9 ? 

98 

99 


cn ' •! 
ro ) 


0200 

0201 

0202 

0203 

0204 
0208 
0206 

0207 

0208 

0209 

0210 
0211 
0212 
02 3 

|ih 

0216 

S!U 

0219 

0220 
0221 
0222 

0223 

0224 
0228 

0226 

0227 

0220 

0229 

0230 

0231 

0232 

0233 

0234 
0238 

0236 

0237 

0238 

0239 

0240 

tV 4 


DO 


05020 

C 


18010 
1 5030 

1 5030 


15000 CON 
C 




15000 J:’l.N)rPL 
|7 MPOR.et.Ol 
3IPLb2 


60 TO 15010 


fO 15020 


ipl«nwpi.*jh 

To 

XPLNX-XPLfllPL 
YPLNXwVP 

tn 

X? 

SfiihE 


PLRx-XPLfllPL, JJPL) 
PLNX*VPL(I 1 PL, JJPL) 
r (((xplnx-x*lst)** 27 (vplnx- 

itiiis iiltilpli Slltr" 

A I C*LCMpfxPLST»Y^LST,l#n 


rPL9T)**23 ,LP. t0*0l ) 80 TO 18030 


05‘ 00 


05110 


05120 
051 30 


051 40 


05150 

C 

06000 


IPOR«* 

NXMb 0 


IPOR 


XL8T««1000. 

SiMPjJ 


"xpT^iS^.NXPi,) go to 05120 


XMprlXM^ll 

F (ixMp.i 00 TO 05120 

F 05150»06050£060B0 


Ir (iXPL.GT.RxPL) GO TO 06040 
IF ( f XMPMxApi.xPLM » 


« I \ % r* • • r . 

NXMaNXHf j 


iiSf 

XI STp 
XH(NX 


00 


ixPL)T.oe'.-', 00 i) 

TO 06000 


TO 05010 
140 


00 TO 06020 


XMP( IXMP) 


I(NXM}bXL8T 
YM^NXMjByqpnxMp) 


CO 70 05130 
08010 NxMbNxM>1 

I MNXM-200) 06030* 06030) 05 150 
06020 NXM«NXM*2 . 

IF (NXM M.200) 

XM(NXM-f )a9Pl 
VM^NXM-li 4 YPI 


06030 XL8T0XPL(2fIxPL3 
XMfNXMlBXLSt 
YMlNXMl.YPL ( 2 , 


00 TO 06150 


C 

06040 

06041 


BTi 


TO 05110 
XMP 
XMfl 


TXPU 


I -j 


PACJR m 
CF POOR QUALITY 





PORTR*»< IV-PLUS VBP-Oie tei34|50 

■ " “ ■ /TR| AlL/WR 


SOPLT .PIN 

0243 

0244 
0243 

024 5 


024 r 

0248 

0249 
0258 
025 i 

92^2 

0293 


IP (NXM, 8 T* 200 ) r.O m 091 B 0 

XM(NXM)*Xvi^(I 1 


YM(NXM)«yMP(|) 

;p (i.lt.nxmM 


CO TO 06041 


060 B 0 NXMPhNXH 
C 

DO i 600 
XMPf 
YMP? 

16000 CONTINO 


!■! |NXM 

;s!i! 



GO TO 
tut) 


1 9«0CT»8? 


PAGE 8 


ro 


oi 

CO 


JKM PAGE IS 


rOUTRAN IV-PLUS V0g-?ir 


30 PlT , r 

C 
Cl 

c 
c 
r 


l*LL/w« 


15i35l?B 


1 S-OCT-82 


PA9C II 


C 00 


0001 

0004 

0009 

0009 

000T 

0009 

0009 

0010 
00 ! 1 
00 ‘ ‘ 
00 
00 
00 
00 
00 
«0t 9 


SUBROUTI'^r PT»* determines SCREEN PLOTTINO POINTS 


§yB?ssir' 


01000 


PTR4(JX, lY.kP.VP) 

Ks{^Y&C.i9C.CAl,B*l(C«2,SA2,NX,NV,FOC*DI8»XW,YV 


JN 0 

CSmMoN 0 
X»XSC*n ... 

i}, 

NaX* 9 *|fY*CAl 
S 0 »H*Cl 2 tOT 8 ^POC 
XTmX*C* 1 *Y* 94 | 
yI>w«SA 24 Z 0 CF^ 

IF (FOC,EQ.0.) CO TO 01000 

00 to 01010 
XP«rxi-XV)<»DTS 
YP»(Yl-YVj»0f S 
01010 RETURN 
ENO 


cn 

f- 


r J 


CRfGirJ'M, PAGE fs 
OF POOR QPAI ITY 




un 


FORTRAN IV-PLUS 
50PLT.FTN 


V0?-SIE 

/TRiAlL/WR 


1SI39I31 


I5-0CT-8? 


PAOE 13 


0001 

0002 

0003 

0004 
0009 
0006 
0000 
0009 
0009 
00i m 

00 ) 1 
001 2 

0013 

0014 

0019 
0016 
001 r 

0018 

3019 

0020 
0031 

00?2 

0023 

0024 
0029 
0026 
0027 

0029 

0029 


*•*» SURnuTiNEi OTAn PI0URC9 OUT 90METHTNS RUT I DON'T KNOR RHAT 


to 

cr 


SUBRnUTlNE OIARflX. IT, INDO) 
dimension D(4B,46) 

COMMON 0,KSC,YSC,Hc,CAl,BA|fCA2,9A2,NX,Ny 

CEN«0. 

NCFN«ft. 

gS J8 IlfiS 

•D(IX.IV-1)-,9*D(IX-1,IY«2) 



01000 

C 


01010 

r 


01020 

C 


ENbNCENI 
IF (IY.CE.NVT no TO 01020 

CENbCEN^I , 9*D(TX, IY)-,8*D(lX-l,IY4l ) 

ncfnbncenIi 

IF (IX.CE.fNXtl )) GO TO 01030 


(lY.GE.NY) 


60 TO 01020 

(TX4i,iy)-,n*0(ix42,iv4i) 


IF 

CENbC.p^ti »y 
NCFNbNCENM 
IF ( lX.OE. (NX41T ) 00 TO 01030 


<1I#I.E.2) GO TO 01030 
CENaCENf 1 ,R*0(IX41 , IY-1 )..BikO(IX42,IY.2) 
NCFN-NCENfl 
01030 CFNbCEN/NC^N 
C 


DC 1«,9* 
, 9* 

Ibs 


01040 

c 


jUi, 

RETURN 




2)) 00 TO 01040 




return 

END 


.AL. p:>ge is 

OF HOOR QUALITY 





cn 

(Tk 


FORTRAN IV-PLUS 




LL'/WW 


1 Bl 351 3B 


f 5«0CT«8? 


PAOE IS 


0001 

000 ? 

0003 

0004 

0005 

0006 
0007 
0006 

0009 

0010 
001 1 
001 ? 

0013 

0014 

0015 

0016 

0017 

0018 

0019 
00?0 

00 ?| 

00?? 

00?3 

00?4 

00?5 

00?6 

00?7 

00?8 

00?9 

0030 

003} 

003? 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 
004? 

0043 

0044 

0045 


**** 8IIBR0UTINF DIN INTT1ALI7C8 ALL' SORTS OF PLOTTING PARANETERS 


SUBROUTINE OIN 

DIMENSION 0(48(40) .NAR(?304) 

COMMON 0. XSC. YSC(ZSC»CA1 (8Ai(CA2(8A2(NF(NY 
COMMON /fill/ VECNUL 


NERR 


3F10.B) 

iklsS?*), 


BAD COMPRE 


90000 format ( 16H1BAD. DATA 
90010 F0RMAt(4l8,“ - ‘ 

90020 F0RMAt(2l5 
90030 F0RHAt(|4H 
90040 F0RMAT(l7H 
90050 FORMAfriOH 
90060 F0RMAt(l7H 
90070 FORMAtnOH 
90080 FORMAT(40H 
C 

NRN ■ 48 
nth ■ 46 

CALL CPEAO(NAR,nCH, NRN, NERR) 
NQHbNRNaNCH 

GO TO 01000 


1 ?) 


h!s'?j!i,"eisn? 

compressed arra 


ON RARANEITER) 
ON) 

ray too LIAROE 


OR TOO SMALL) 


IP (NERR,EO.0) G( 

.UoP 


01000 DO 


10010 

C 

1 0020 
10030 
1 0040 
1 0060 


1 0060 


1 0070 
< 0080 

1 0090 

C 

10100 


Rf A0( 1 (90&10) NYIN.NY1,NYF,NCY.Y8C 
IF (NYI) 10010, 100?0> 10030 

Type 90070 


STOP 

NVI»1 
IF (NVF) 
NYF«NY jN 


\l 

i?sf ” 


100101 100401 100B0 
00 


CT.NYI 


) 00 TO I 
N) GO TO 


010 

iw 10010 

Y) 10060, i0070f 10060 
900B0 


-NYI 


NCY-(NVF. , ^ 

NY«(NYp*Ny IH )/^CT 


♦ 401/48 

1)/NCY 


IF (Ny-CT*40) 00 TO 1 
/F lI.ftT.f) GO TO 100 
NXINbNVIA 


10100 


IN 

ny!»nyi 

NCX-NCY 

NXhNY 

XSC-YSC 


0090 

00 


o r? 

•Ti . 

' 

c 

c* 



I 


r ’ 
—! 

■< 


ro 



FORTRA'J 1V-PLU3 V0g-9lE 
30 «»LT.FtN /TR|*lL- 

R04S 1000R CONTINUE 


1BI35I3S 


1 B-OCT-S? 


00F0 
00F i 

00F4 

007B 

00F« 

0077 

0070 

0P7IJ 


01100 XSCb 
Y sr* 
01110 IF ( 


01100 READ 


\i 00 TO 01100 

If (ysc.Bt,0,!1 go to 01 I 10 
xsc-i. 

Y9r»?fiY.l,)0NCY/f(NX«il,)*NCX) 

IF (NCH,|0,NxIn*NYTN) ftO TO 01120 
90040 

Z9C«X9C 


01130 00 11000 J-l.NY 

00 11000 IbI.Nx 
NTa0 

IPaNxl-Uf l-n*NCX4TNYl*lA( J-1 )»NCY)*NXIN 


DO 12000 JCb 
00 12000 


1 2000 
C 


CONTIN 


12000 iCBl.NCX 

JTbNT»NAR(IP4IC*NYIN*JC»NXIN) 


NAR( IfNXBj-NXTaNT 
11000 CONTINUE 

_ • 


01140 


IF (f(ALl-RB.)* 
90030 


rALI-9.)),LC.0.) 00 TO 01IB0 


011B0 JF *JT*LJ-8B‘^JajAL2-8,)),GT,0,) 60 TO 01140 

ciUclji{AL.|? 

9Al*8IN(Ar{j 
AL2-AL 2/5F-29B70 
C42«rOS(IL5) 

SA2B9TN? ArpJ 

IF (NROt*(NftOT-3l .LE.0I GO TO 02000 
'VPE 900B0 


02000 IF 
IF 
C 

02010 IP 


(NROT’.LE 

iNROt.Jl 


01 CO TO 02010 
02020*02030,02040 


Pi ■! 

P2-NX 
P3B-NX 
0 TO 03000 


02020 NTbNX 
NXaNY 
NYbNT 
9CT«*8C 

xScbY8C 

YSC"9Ct 





PAOE' 17 


O C 
■n jc 

T ?2 
O - 

r 


o --J 

o 

r; r 1 

H — f 

-< W 


K) 

'O 


1 5-0CT-8? 


P«QE 18 



( 


FORTRAN IV-PLUS V 

idplt.ftn / 



/MR 


0140 

SIS) 


ZSC*ZSC*VNA)( 

Returh 

FND 


19139135 


o 


IT\ 

■U I 


OF POOR QUALITY 



CT\ 

O 


POKTRAil IV-PLUS 
JBRkC ,flN 


RPR I 


0002 

0001 

0004 

0000 

0006 

000T 

0006 

0009 

0010 
0011 
001 ? 
0011 
0014 

00 5 
0016 
0017 
0019 
0019 
00?8 

0001 
00?J 

0023 

0024 

0025 


0028 

0027 

0026 

0029 

0010 

0011 

0012 

0011 

0014 

0015 

0016 
0017 
0016 
0019 


V02-51E 

/TR|*lL/WR 


1 Bl 37| 32 


16 -OCT -02 


P49e 1 


•••* JBRKC - Cal 

*••• Of 


KC - Calls jbrkfn TO nsp current dens 

THF BtRKELANB CURRENT ROOEL OEFINEO BY 
CURplS AT POINTS ON A SPHER [ OP RAOIUS »A 
•••• OrFATER Than that op THE> EARTH'S. 


ENSITy 


ALT' 


S OMMON 
ZP (2.6 

DATA P! / 
DATA ALTI / 


/PILE/ CLl .CL 2 .NUMT.NUHL, RP.TP( 4 . 6 , 721 ,ARPSf 2 ,B, 


72 ) 


3. Hi brSb^ 

■ 140RRR. 


ALT ■ 490000 *. 

OPEN 
READ 
READ 
READ 



■ 1 , name ■•pis', DAT ',TY PE ■'OLO '1 
nCooe. op 


O O 

^ 5? 
tj n 
o ^ 
o 
:c 

i- T 
C “i 

•< Ui 


OPEN fUNl T" 2 ,NAMEb'OR 0 iJPIL,OAT' , TYPEb'NER' , 
1 ACCESS* 'DIRECT ',REC 0 RDSI 7 £«I) 

( 2 'i) AIT 

1 10 • 1, B0. 

-■■SiSrss! 


"RITE 
00 1 ID 

> ■ -nP*00,/49. 
SCM ■ SiN(CmJ 
COM B cosrcM^ 


• 40 , - BB,/ 49 ,)*P 1 / 180 ', 


DO 


ic 

SOC • SIN(GC) 

CCC ■ CPS(GM 

YL I (re ♦ 

ZL ■ (RF- ♦ ALT) •CGC 
BRA a FUN(YL, YL,!L) 
IT ■ 1 ♦ IE ♦ (!0 • 


IE 


■ 1 * 60 
( 11 * 80 . / 49 , 


I 2 


• 40 , . BB./ 49 .)*Pl/< 60 , 


*SQH 
• COM 
1 )*B 0 





,r 


ri^ 


IV-PLUS V0 


/T 


i:Vd 


LL/HR 


15137132 


SSI? 

RP42 

0043 

0044 


coniin!!?”' 


END 


{ 5-0CT-8? 


PAOE 2 


ORIGINAL PAGE IS 
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,1 FORTRAN IV-RLUS 

JSRKTN.rtN /TR|ALL/R« 


16I37IB0 


I 9-OCT-0? 


PAOe I 


Ct\ 

ro 


0RRt 


903? 

0003 


0004 

0005 

0006 
0007 

0000 

'^009 
0 
0 

00 ) 
00 
00 
00 
00, 
00 ) 
00 
00 


I 

00^0 

00?1 

00?? 

00?3 

00?4 

O0?9 

00?8 

30?7 

00?R 

00?9 

0030 

003) 
003? 

0033 

0034 
0039 

0036 

0037 

0036 

0039 

0040 

004) 
004? 


rUNCTION F«»Nf XL.YL.ZL) 


•*** JHRKFN CALCULATF9 CURRENT DENSITY IN E«A CURRENTS 


1? 
I 1 


DIMENSION RM(3,.3), EAMR(6 
COMMON /FILE/ ClI 
1RZF(?,6),46(?. 

?9C' i?,6j,CCLri 
REiL*4 lambda, 

DATA RE / 6371 
DATA ALTI / |4 
PI / 3,14 


DATA 


1000 .. / 

.J 000 A, / 

14199?60 / 


MP(6.F?) 


N ■ NLMT ♦ I 

][^."nuAt ♦ 1 

Do 10 M p 1,? 

If fM.EA.?) I 

Do 10 N ■ 1, 


NUMT 


IS 




J? ^M%ij?)^3!i TO 1 ? 

FAMpSIumKi# J) 

If 1MmP?N» J) .EQ.f. 3 80 TO il 

ILA ■ SlfiK^ABDA) 

CLA ■ COSfLAMBDAj 


AMRn,N-l, J1 


AMP(?,N, J) 


i/numl 


RM 

RM 

RM 

RM 

RM 

RM 

RM 

RM 

RM 


J:? 
?»?' 
?»31 
3, ) ) 
3i?) 
3,31 


CLA*CMUfM,N) 

8LA*CMUiM,M) 

•9HU(M,N) 

-SLA 

CLA 

CLA*?Mlif M, N) 
SLA«9MU?M,N) 
CnU(M,Ni 


Itlll ■ 


<F ■ FL«RM 
YF ■ XL*RM 
ZF ■ XL«RM 

RCF ■ XF**2 ♦ YF**2 

RBF ■ 90RT(RCF ♦ (ZF-RBfMfN) )••?) 


u> 

u> 


ORIGINAL PAUu fc 

OF POOR QUALITY 


1 

( " 

I 


IV-PLUS V02.5lf I5l37|50 19-0CT-8? PAOt. 2 

JBRxrN.rtN /TB|ALL/»<R 


0043 

004 4 

0049 

004S 

004 / 

0040 

0049 


C 


r 


1 

CONTINUE 
FUN ■ JT 

return 

end 


RTF 


90RT(RCF ♦ (ZF-RT (V(,N) )*•?) 

40.) CO TO 10 
i)*RCFT)«*?n ♦ JT 




03 

Ui 


original page !S 

OF POOR QUALITY 



FQRTRAV 

jbrkp.ft 


JV.PLUS 


ywR 


1BI3RIP9 


!S-0CT»8? 


PAOCi 1 


I 


* * * ft 
C ft ft ft ft 


JBRKP • PLOTS DENSITY AS ELEVEN COLOR FIELD 

USING direct commands TO TEKTRONIX 4027 


0001 

0002 

000S 

0004 

0009 

0000 

0007 

000B 
0009 
001 0 

001 
00 ! 
0013 

00*4 

00l9 
0016 
0017 
001 S 

00I9 

P020 

0021 

0022 

0023 

0024 
0029 

002fl 


0027 

0028 

0029 

0030 

0031 
0022 
0033 


0034 
0039 
0036 
00 37 
002b 

0039 

0040 


DJMENSJON jPOLie)# !PC(3)» IPA(14), IVf4) 

^eIlm MAxIK, MaXOUT 

DATA IPA / l>0tB9« 170,B9» 1 70,89, 170, 06, 1 70, SB, 1 70* 89, I 70, 89 / 
DATA AE/637100ft,/ 

DATA PI/3,141B<»20B/ 

DATA R^ / iRFt / 


iS. 

c 

1 1 0 

r 


rtRlTE (9,1001 

FOPMAT ( I select COLOR RANGE FACTOR • ) 
READ (9,*) Rr 

WRITE (?»110> . 

FORMAT { I select delay BUFFER < ) 
read (9, a) IDEL 


CMAX ■ 8, 
MAXIN ■ A,. 
MAXOUT ■ &. 
CPI * 1,. 

CP2 ■ I, 
CP21 ■ I. 

CP31 • I, 
CPA ■ 1 , 

CP41 ■ f, 
CP0 ■ 1, 


0 (?'l) ait 

CMaX • AMAX! 


read 

no 


2901 
C _ 


MAX] 

MAXI 


200 

continue 


C 

ftftftft DEFINE COLORS 

r 

1 02 

WRITE (9,102) 
FORMAT MIH IPON 
WRITE (8,103, 

4 H K) 

103 

FORMAT M0H iCRA 
WRITE (6,1000) 
fORMAf (20H IMAP 
WRITE (9,1001) 

1,30) 

1 00W 

CO 0 


0,100,100) 


ORfGliViAL PAGE 13 
OF POOR QUALITY 


PO^TRAM IV-PLU3 
JB'IKP.FTN 


0054 

0055 

0056 
0037 

0058 

0059 

0060 
0061 
006? 

0063 

0064 

0065 

0066 
0087 
0066 
C069 
0070 
007 1 
007? 

0073 

0074 

0075 


1 007 
C 
C 
( 




16135109 


1 5.0CT-8? 


FORMAT (?0H IMAP Z\ 180» 60,100) 
WRITE fS.l00?) 

format (^0H IMAP C2 140, 80,100) 
WRITE (5,1003) 


WRITE (5,1003) 

FORMAT (?0H IMAP C3 00, 60,100) 
WRITE (5.1004) 

format (?0M IMAP C4 220, 60,100) 
WRITE (5,1006) 


format (20M IMAP C4 220, 60,100) 
WRITE (5,1006) 

FORMAT (?0H IMAP CB 300, 40,100) 
WRITE (5,1006) 

format (20H -MAP C6 330, 40,130) 
WP'TF (5,1007) 

format (20H IMAP C7 0, 0,i00) 


(20H -M 


IMAP C7 


OEFTNF PATTERNS 


WRITE (5,2000) TPA 

format (15H IPAT P0 C0 CO , 13(13, IH,) 
WRITE M,2001) !PA 

format M6H IPA PI Cl Cl , 13(13, IH,) 
WRITE (8,200?) PA 

FORMAT (15H IPAT P2 Cl C2 , 13(13, IH,) 
WRITE (5,2003) IP* 

FORMAT nSH IPAT P3 C2 C2 , 13(13, IH,) 
WRITE (5,20^4) TPA 

format (IBH iPA* P4 C2 C3 ,13(13, IH,) 
WRITE (1,2005) IPA 

format ISH IPA PO C3 C3 , 13(13, IH,) 
WRITE (! ,2006) PA 

format I I5H iPA? P6 C4 C4 , 13(13, IH,) 
WRITE (8,200h PA 

FORMAT (15H IPAf P7 C4 CB , 13(13, IH,) 
WRITE (5,2008) IPA 

FOftMAT (IBH IPA P8 C8 CB , 13(13, IH,) 
WRITC (8,2009) PA 

format (Ibh IPAT P9 CB C5 , 13(13, IH,) 


WRITE (5,2006) 
FORMAT (I5H iPA 
WRITE (8,200h 
FORMAT (IBH iPA 
WRITE (8,2008) 
FOflMAT (IBH IpA 
WRITC (8,2009) 
FORMAT (18H IPA 


FORMAT (I6H IPAT P9 CB C5 , 1 3( I 3, IH, ) , 13) 
WRITE (5,2010) IPA 

format (i6H IPAT P10 C6 C6 , 1 3 ( I 3 , 1 H , ) , I 3 ) 

***• graph data 

DO 300 10 ■ 1, B2 
NF • 90. 

IP (10 ,CT. 80) NE ■ B 
JY ■ ID*B » 4 


on 300 IE P 1, NE 
IE (10 .Of. B0 
JX • IFaB 


B0) 00 TC 32 


**** oepine pixel parameters 

ISStHi : J! : J 

POLi'B) " JX ♦ 4 

®0L(7) ■ JX ♦ 4 

PnL(2) ■ JV • 4 


PARE' 2 


O C 

■n 



O TJ 
C Jt> 


> O 


r- rn 


H — 
■< w 



136 



. FORTRA'i IV-PLUS 

r^t jbrkp.ftn 


cr 

C^l 


00SS 

00B9 

009{| 

0091 

0092 
0L9 3 

0094 

0095 
0095 

0097 

0098 

0099 

0100 

Wi 

0103 


0104 


05 

06 

07 

08 
09 

0 


0115 


11 

8 

9 

0 

01 

?? 

23 

24 

^6 

27 

28 

29 

30 


C 

32 


C 

51 

30 

C 

c 

c 

c 


34 

3B 

38 

37 

33 


30! 

302 

3301 


303 

300? 


304 

3003 


305 

.3004 



16138109 


15-0CT-8P 


PAOEi 3 


« « * * 


V ♦ 

)Y ♦ 
jy » 


570 t nD*90)»10 

he-i)«»2« ♦ 21 

870 ♦ (IO*S0)*ia 

It*28 ♦ 21 

570 ♦ (ID«49)*10 

i/5^t (lSi49J»10 
♦ 21 


ME-t)*J8 


JyU^'bij^uJtfkEL) 


no 

CON1 

IF. (10 '.CT', 50) CO TO 33 

CONDITION D4TA - ADJUST DYNAMIC RAN9E, CONVERT TO INTE<iFR 

II ■ 1 ♦ IE ♦ (ID - 1)*50 
read (2»lt) C 
IE (C/CmAxI 34,39,36 
jPZ ■ «|NT(ALr 


iiz]i i 


hn\ 


•RF*C/CMAK) ♦ ,B) 


IJ fig ;??; IS! n 


B0h*IE) 


« * * * 


* * * * 


3FT COLOR FOR PIXEL 
IE (IP?) 308,307,301 
POSITIVE 

5) 


l5'?3«l'Se 

,TE (8.300 

Imat (8H 

AMIN! 

TO 314 


303,304, 


IE 

BO 

WRITE 
EORHA 
CPI 

CO IV j I 

WRITE f9.302?) , 

format (8h KOI, P? 
CP2 ■ AMTNl (CP'?,C) 
CO TO 314 
NRItE (3,3003) 
FORMAT (8H KOL P3) 
CP3 ■ AMlNl (CPS,- 
CO TO 3i4 


00 TO 308 
309), (PZ 


C) 


EO 


P4) 


CRiGI.'iAL P/A'ic 


'S 


CTki 


FORTRAN 1 V-PLU 3 
JSRKP.rTN 

II 


34 

38 

38 


cil 3 r 


in 

(314 


8 

9 

4(3 


n 

45 

48 

4 ? 

48 

49 
5 P 

2 > 

52 

53 

54 

55 
58 


0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

04 57 

0158 

0159 
0180 
0181 


ini 

0184 

0185 


0188 

0187 

0188 
0189 


306 

3008 


307 
3000 

C 

C 

c 

308 

309 
3008 




07 


3009 


313 

3010 


C 

?S1 

c 


108 

300 

C 

C 

c 

c 


'‘M'.in 


18 I 3 RI 09 

ll/wr 

CP 4 ■ AMIN 1 (CP 4 »C) 
00 TO 3|4 
PRITF ( 5 . 3005 ) 
format ( 6 h ICOL P 5 ) 
CP 5 n AMINI (CPI,C) 
00 TO 3 l 4 


18 - 0 CT- 8 ? 


PAOE 4 


00 


•*** ZERO 


foJIat^ P0) 

no TO 314 


314 

***• NEOATIVE 

IF- (IPZ ’.L^i *91 09 TO 31 J 
RO to ( 3 ^ 9 , 510 . 311 . 312 ), -TPr 

foJmJt^JAm icJl P6) 

CPI P AMINI (CPI , ABS(O) 

00 TO 314 
NRltE ( 8 . 3007 ) 
format ( 6 h ICOL P 7 ) 

CP 2 P AMINI (CP 9 , ABS 

format (6H TCOL P8) 

CP 3 p AMiNt (CPS, ABS(C) ) 


(C)) 


»4RltE (8,3000) 
format ( 6 h 1 C 6 L' P 9 ) 
CP 4 P AMtNl (CP 4 , AbS 


i 


0 TO 314 

HRItE (8,; 


(O) 


8 , 3010 ) 

format (OH ICOL P 10 ) 

CP 5 ■ AMtHl (CP 8 ,AftS(C) ) 

***• SET PIXEL SIZE 

» 1 H,), 13 ) 

IE no ’.RtI 50 ) GO TO 300 

IF. ((C .LT. MAXIN) ,AH 0 , (C , 0 T, MAXOUT)) 00 TO 300 
**** PLOT pixel 

5 H ILIK) 

continue 

rL 05 F (UN 1 TP 2 ) 

**** PLOT LATITUDE CIRCLES AND MLT L 3 NES 
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6. 


CD' 


JB9kp , 


V-PLU3 V 


15l3ni09 


si;? 

SIH 


^4 

5? 

78 

79 
89 
8t 
82 
81 


SISS 

0103 

a 1 8? 
9t 8n 
0)89 
0)99 
0 91 
0)9^ 
0191 

0)94 

0191 

0198 


0197 

0)98 

0199 

0?00 

0201 

020 ? 

0201 

0204 

0201 

0208 


40 

C 


4000 

4001 

4002 

4003 
400 

C 

9000 

C 


10 


9001 

C 

900 

C 

C 

r 

8000 

8001 

8002 


8003 

c 


;?5;5JE.«r 

no 400 I ■ 

on 40 K ■ 1 . TOIL 
§jr ■ rufum(Ioel: 

tINUP 


19-0CT-82 


P40e> 9 


■ I. 4 


CONI 


J ■ 1l* 90 
WRITf 


WRITF (5.40001 _ 

FORV44T (7H IlIN P) 
WRITE (9.4001) 


FORMAT 


?I7H 


WRl|T jT<*t400 


IVEC 0,0,204,209) 


WRITE 
FORMAT 
CONTINUE 


(^H IlIn I) 
(9,4003$ 3 
(8H ICIR ,13) 


WRITE (9.9000) 
format (7H iLlN 


F) 


DO 900 



R 1 


3f'^b*BiVuJ (fl$EL) 


NUF 

■ 

p 

m 




|!SJIlT'?t3‘’?»c’'' 


,3(I3,1H,),I3) 

CONTINUE' 

•*** pRITE headlines 
WRITE (9,6000) 


nn 1 I C I 3 I OPnw I 

FORMAT M8H IwORIMAR S4|UP 28) 
WRITF (9,6001 ) A “ 

)Amat (/8H 6t 


FORMAT (/8H 
ill,' ?///8H 




WR 

FO 


//3H AT,)PEB 
6002) MAXIN 


Y/3H 0F/19H 8IRKELAND CURRENTS 


I 


WR 

FO 


BH 1 


IATT CS/10H MAXIN p i 1FE8,1,7H A/M**2/ 
ATT C6/10H MAXOUT p ,lPE8,f,7H A/M*»2) 
CP9 ,cM,CP3,CP2,CPI _ 


:sj 

TT O 

33 r~ 
O -n 

^3 


ITE (9,6003) CP9.CF4.CP3.CP2.CPI 

RMAT (14H lATT C0IOOW 6 , S ( // 3x , ! PE8 , 1 , 7H A/Mp*2),9H IMON) 


STOP 

END 




FORTRAN IV-PLU3 V02-5IF 
BUFuW.FTN /TR|ALL/WR 


I9I39IR1 15-0CT-B2 


C ***« PROVIOF9 A time UAO BUFFER FOR PLOTTING ON A 4027 
C 

I * 

FUNCTION BUFUNflOCL) 

BUFR ■ 3IN(IDFL/10,) 

BUFUN a TANHfBUFR) 

RtTURN 

FNO 


CT> 

•O 



OF POOR QUAUrt 


